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THIRD-ORDER RATE CONSTANTS OF ATMOSPHERIC IMPORTANCE 
Roger Patrick, and David M. Golden 
(International Journal of Chemical Kinetics, 1983, 15, 1189-1227) 


Third-Order Rate Constants of 
Atmospheric Importance 

ROGER PATRICK* and DAVID M. GOLDEN 

Department of Chemical Kinetics, SRI International, Menlo Park, California 94025 


Abstract 

Input data and results are presented for the calculation of a number of third-order rate 
constants of atmospheric interest using Troe's approximate method. A comparison with 
experimental data indicates that this approach provides a reliable method for predicting 
unknown rate constants and estimating temperature dependences. These calculations form 
the basis of the recommendations of the NASA review panel for third -order rate constants 
to be used in atmospheric modeling. 


Introduction 

Quantitative modeling of the atmosphere requires a comprehensive set 
of rate cons .ants for a number of chemical reactions. Many of the most 
important reactions are association reactions which are either at or close 
to their limiting third-order kinetic behavior under atmospheric conditions. 
The review of laboratory measured rate constants for these processes is 
currently undertaken by a number of panels, and their recommendations 
appear at regular intervals. Unfortunately few measurements have been 
made at stratospheric temperatures close to 200 K. and so it is necessary 
to extrapolate from data obtained at room temperature and above in order 
to obtain the relevant rate constants. In the past this has been done by 
assuming that the experimental rate constant data followed either Arrhe- 
nius, that is, k = .4 exp (-E/RT), or T n type behavior. Neither of these 
forms has any real theoretical basis, and although each may fit experimental 
data equally well, extrapolation at 200 K using the two equations may give 
differing results. 

In principle uiumolecular rate theory offers a sounder theoretical basis 
on which to fit data and formulate extrapolations. However, rigorous use 
of detailed theory presents a formidable computational problem, especially 
when the evaluation of large numbers of rate constants is required. 

* Postdoctoral Research Associate. 
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Troe [4,5] nas recently suggested simplified equations, based on RRKM 
theory, allowing the quantitative estimation of second-order rate constants 
for unimolecular decomposition which may be related through the equi- 
librium constant to the rate constants for the reverse third-order association 
reactions. In a further paper [6] he discusses treatment of rate constants 
in the fall-off region. This method has been adopted by the NASA and 
CODATA review panels [1,2] to estimate theoretically the temperature and 
pressure dependence of key atmospheric reactions, and on this basis, to 
evaluate and rationalize laboratory measurements. Th ' details of these 
calculations, along with the required thermodynamic and spectroscopic 
input data, are given in this paper. 

In addition to providing an improved method for extrapolating experi- 
mental data, these calculations allow a thorough testing of Troe’s formalism 
for a wide variety of molecules, many of which are larger than those that 
have been treated by this approach before. 

The formalism also provides a useful method of highlighting uncer- 
tainties and ambiguities in the experimental data. 


Theoretical Formalism 

The rate of dissociation in the low-pressure second-order regime is given 
by the rate of energization of molecules above the dissociation threshold 
Eq. If the molecular energies are separated into their vibrational and ro- 
tational contributions, then the strong collision rate constant may be ex- 
pressed as 


(1) *dUo = Z u C m f(E,J)dEdJ 

%J 0 Eq 

where Zu is the Lennard—Jones collision frequency, and f(EJ) is the 
distribution'of molecules among the available rovibrational states. Troe’s 
[4,5] formalism starts from eq. (1) and by a series of approximations fac- 
tors out contribution to f(E,J) by differing degrees of freedom: 

<o\ l sc _ 7 p(Ep)RT IzMc c v c 

(-) ^disa.O = ^LJ q GXp j* £** anh^fot** corr 

where p(£ 0 ) is the density of states at the critical energy £ 0 ; Fe,F. ,„h. and 
F rot are the correction terms for the energy dependence of the density of 
states, anharmonicity, and rotation, respectively; and Q vi b is the vibrational 
partition function. 
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Following Troe, Zu is given by 

(3) Zu = 8.09 X 10~ 10 /cm 3 molecule -1 s -1 

X >J (T/lOOO K)(20 g mol -1 /MAA/) <^«/5 A) 2 Q*& 


where the collision integral is approximated by 

(4) a* [0.697 + 0.5185 \og(kT/t AM )]-' 


In this expression the reduced mass v-am is given by m A mm/(m A + m.vf). 
the Lennard-Jones collision diameter a AAt is gi ven by (a A + n M )/2, and 
the Lennard-Jones well depth ( A m is given by V ( AA (mm 
p(Eo) is calculated according to the VVhitten-Rabinovitch approximation, 
and Fe and F an h are given by 


(5) 


*-> (s - 1)! / RT W 

u-o (s - 1 - o)! Uo + aEz) 


with a the Whitten-Rnbinovitch parameter and E z the zero-point en- 
ergy, 



with s the total number of oscillators and m equal to the difference in the 
number of oscillators in the reactant and product molecules. 

The barrier to association for the reactions discussed here is likely to be 
negligible, and so the critical energy E 0 may be identified with A/i o, the 
enthalpy change for the reaction at 0 K. Under these circumstances the 
relevant expression for F rot is 


(7) F rot = 


(■s - D! I E o + q£, \a/2 
(s + V 2 )! RT , 


2.15(g„/flD 1 ' 3 

2.15(£ 0 /FT) l/3 - 1 + (£„ + aE : )/(s -P '/ 2 )RT 


The treatment of internal rotors presents some difficulties, since for many 
of the molecules considered here, barriers to rotation are not known. 
Consequently all internal rotors are treated as low-frequency torsional 
vibrations. At the low temperatures relevant to these calculations, the 
difference in the heat capacity of these types of motion is not signifi- 
cant. 

The term F corr is not calculable and represents a correction term to ac- 
count for coupling between the different types of degress of freedom. 
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Following Troe, this was set equal to unity. This step may be justified a 
posteriori by the good agreement with experimental measurements found 
in most cases. 

The association rate constant may be easily obtained from /jJum.o 
via the equilibrium constant K^: 

(8) *Sc,0 = *&..0 Keq 


fc£c,o represents an upper limit for the rate constant and must be multiplied 
by a collisional deactivation efficiency d c (0 < /3 C < D to obtain a value that 
can be compared to experiment. The temperature dependence of /i c is 
given by 


(9) 


Pc <A£) 
1 - 0\ n ~ F e RT 


where < A E) is the average energy removed per collision. 

In the fall-off region of transition from high-pressure second-order to 
low-pressure third-order kinetics the pressure dependence of association 
reactions can show considerable deviation from classical Lindemann be- 
havior. Luther and Troe [6] have developed a simplified equation which 
accommodates this deviation in a parameter F, the broadening param- 
eter, 


( 10 ) 


k(M) = 


kv, -I- k 0 [M] 


F exp 


1 + 


^log 


[Ao[A/] 

2 

ao 

. 


where kc, and k„ are the low- and high-pressure limiting rate constants, 
respectively. Where F = 1, Lindemann behavior operates. F is made up 
of the product of two terms: a strong collision, F K and a weak collision F wc 
broadening parameter. F sc is dependent on the molecular structure and 
in general becomes smaller as molecular complexity increases. It may be 
calculated according to rules derived by Luther and Troe [5]. F wc is given 
approximately by 


(11) F wc 

For most of the molecules encountered here, a value of F ~ 0.6 is found, 
although for some of the larger species it may be considerably less. 


Table Format and Input Data 

The following units are used throughout the tables: i'/cm -1 ; 7/amu-A 2 ; 
77 K, A/7//kJ/mol; S°/J/mol-K; K„,/molecule/cm -3 : cru/A; p/cm; 
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Zu/cm 3 /molecule*9; fcdi»«o/cm 3 /molecule-s; fc r , c o/cm 6 /molecule 2 -9; 

(A£)/kJ/mol. 

Thermochemical data are generally taken from JANAF tables [7], except 
where footnoted. Data in parentheses are estimated either on the basis 
of known thermochemistry and spectroscopic parameters or from additivity 
rules [8]. Spectroscopic parameters v and / are also generally taken from 
JANAF. If frequencies are estimated, they are chosen so as to give the 
correct entropy for the molecule. Unknown moments of inertia are cal- 
culated on the basis of reasonable or known molecular structures using the 
method of Thompson [9]. 

Equilibrium constants are calculated from both the thermodynamic data 
K™ and the spectroscopic data in order to check the internal consis- 
tency of the input data. Agreement between these values is generally very 
good. Where there is disagreement, this usually results from the neglect 
of states other than the ground state in the calculation of the electronic 
partition function. For this reason the thermodynamic equilibrium con- 
stant is always preferred. 

The Lennard-Jones parameters, a and dk are generally not available 
for the molecules of interest. These are estimated from the relations 
[10| 

(12) a = 1.18 Vl* 

(13) dk = 1.2m 

where Tb is the normal boiling temperature (K) and Vb is the LeBas value 
obtained from tables of atomic volumes. Where Tb is unknown, it is esti- 
mated on the basis of analogous molecules. In all the calculations, the 
third-body M is taken to be nitrogen for which the values a =■ 3.61 A and 
dk = 91.5 K are us^d [9]. 

The value for k^ M at 300 K is generally the one r ^commended by the 
NASA review panel, and dr is derived by comparing this with the calculated 
value of the strong collision rate constant k ( A £ ) is then obtained from 

eq. (7). Following Troe [4] and assuming (A E) >s temperature indepen- 
dent allows dr at 200 K to be calculated and combined with k^ 0 to give the 
predicted value of k ?£ ,o at 200 K. ( In those instances where no experimental 
data exist for a reaction, a typical value of 0.3 is assumed for dr at 300 
K.) 

The temperature dependence of may be cast in the form 

(14) *r*.o-*^.o(7V300)-" 

The appropriate vaiues of n are given in the tables. 
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Sensitivity and Accuracy of Calculation 

The sensitivity calculations refer to the HO 2 NO 2 = HO 2 + NO 2 system 
which is fairly representative of the moderately complex species considered 
here. 


Variation of A Hq 

Since the exponentia' terms in eq. (2) and the equilibrium constant cancel 
when &£c.o is calculated, its dependence on the dissociation energy is rela- 
tively weak. An increase in AH JJ of 10% is reflected in a factor of 2 increase 
in o This arises almost entirely from an increase in the density of states 
at the dissociation limit p(E 0 ). Since HO 2 NO 2 has a number of low 
frequencies (125, 200 cm -1 ), the dependence of p(Eo) on A H” is more 
pronounced than it would be if all the frequencies were high. 


Variation of Vibrational Frequencies 

p(E 0 ) is the term which is most affected by changes in the vibrational 
frequencies. The extent of the variation depends greatly on the size of the 
frequency under consideration. Variation of frequencies close to 1000 cm -1 
by a factor of 2 gives rise to a variation in of ~20%, while a similar 
uncertainty in frequencies ~200 cm -1 gives rise to a change by a factor of 
2. In most cases the vibrational frequencies of the molecules are well 
known. Where they have to be estimated, the assigned values can be 
checked for consistency by calculating the overall entropy of the molecule 
and comparing this, via group additivity, to other species. 


Variation of Lennard-Jones Parameters 

These are the most uncertain of the input parameters. However, the 
values estimated from eqs. ( 12) and ( 13) are unlikely to be in error by more 
than ~20%. An error of 20% in (/k only leads to a 3% error in k J£ >0 for HO j 
+ NO_>, while a similar error in a is propagated more directly as a 23% error 
in the rate constant, assuming no errors in the corresponding values for 

N a . 

Thus it can lie seen that is comparatively insensitive to most of the 
input parameters and so, even when data are estimated rather than known, 
little error should be introduced into the calculation, provided reasonable 
choices are made. 
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Accuracy of Eq. (2) 

In addition to considering the uncertainties introduced into the calcu- 
lations via the input data, it is useful to try to assess the magnitude of the 
uncertainty associated with the approximations in eq. (2) itself. It can be 
shown [II] that eq. (I) may be written as 

• 15) ^-..-^J £ '7(£)p(£)exp(=D<i£ 


where Q r ' and Q r are the rotational partition functions for the transition 
state and the molecule, respectively, Q v is the vibrational partition function 
for the molecule, p(E) is the harmonic density of states at energy E, and 
y(E) is the appropriate correction for anharmonicity. The factor Q r f /Q r 
represents a correction for rotational effects and depends sensitively on 
the form chosen to represent the potential surface. The corresponding 
factor Frot in the Troe formalism is derived for a Lennard-Jones interaction 
potential with a centrifugal barrier, and following this assumption, it can 
be shown [8] that for the diatomic approximation 


(16) 


Q / (6Do| 1/3 

Qr “ Ur, 


where Do is the potential energy of interaction. 

Equation (15) was integrated directly using numerical integration 
techniques and input data for the dissociation reactions of HO-jNO-j and 
CINOo. The strong collision rate constants thus calculated were 5.8 X 10 -21 
and 1.6 X 10 -29 cm :i /molecule-s, and these may be compared with the values 
of 1.2 X lO -20 and 3.2 X lO" 29 cm a /molecule-s given by eq. (2). Thus the 
rate constants calculated by the different techniques agree to within a factor 
of 2, although it should be borne in mind that neither eq. (2) nor eq. ( 15) 
correctly addresses the equation of rovibrational coupling. 


Notes on the Tables 

These notes referring to Tables I— XIII. address specific points raised by 
the calculations for some of the reactions. 

Table II: .V,0 = N, + 0‘D 

The rate constant 1 12] for this association of N '2 and O'D is extremely 
low compared to the prediction for k^. This is a consequence of inter- 
system crossing to the lower energy N 2 + ( VP curve, providing an additional 


1196 


PATRICK AND OOLDEN 


dissociation channel. This channel cannot be accommodated by the simple 
Troe model, and so the dr value derived here does not represent a simple 
collisional deactivation efficiency as in the o«.her examples. It was assumed 
to be temperature independent in deriving the suggested value for at 
200 K. 


Table V: NO, = NO + 0 

Third-order calculations by both Troe [4] and Michael and Lee [13) also 
give rise to values of k^ 0 that are smaller than the recommended experi- 
mental value. The reason for the implied dr value > 1 is not clear since 
both the spectroscopy and the thermodynamics of this system are well 
defined. 


Table VII: H0N0 ^ OH + NO 

The recommended experimental value for the rate constant for the as- 
sociation reaction is 6.7 X 10 - - 1 cm 6 /molecule--s. In order to reproduce 
this value correctly using the Troe formalism, it is necessary to take account 
of the formation of both the cis and the tran s isomers of HONO. The value 
of dr was assumed to be the same for each of the channels. 


Table IX: HO , NO , = HO, + NO, 

The vibrational frequencies assigned to HO 2 NO 2 are those recommended 
by Baldwin and Golden (14). The first five frequencies are obtained from 
work of Niki et al. [15], while the remaining ones are based on those for 
FNO:i, loO-i, and HN0 ;1 , with the low-frequency torsion adjusted to give 
the correct entropy for the molecule. The entropy of HOjNOo <s based on 
a measured value for MeONO? (16), corrected for the change of internal 
rotor from a methyl to an OH group. The bond dissociation energy at 0 
K is that recommended by Baldwin and Golden (12|. 


Table XII: CINO, ^ Cl + NO;, C’ONU = Cl + NO, 

The calculated strong association rate constant for the formation of 
CINO.. is 7.71 X I0“ :u cm'Vmolecule J *s. which when compared to the ex- 
perimental rate constant for the disappearance of Cl. 1.8 X 10“ 10 cm"/ 
molecule J -s [1|, would lead to dr * 2.35 if this were the only product. This 
value for the deactivation efficiency is some eight times larger thun might 
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be expected and well outside the uncertainties of the calculation. However, 
Niki et al. [17] observed that CION D is in fact the major product and that 
the two isomers are formed in the approximate ratio 4-1 C10N0-C1N0 2 . 
The strong collision rate constant for CIONC formation is calculated to be 
2.83 X 10* 30 cm 6 / molecule 2 *, based on vibrational frequencies reported 
by Janowski et al. [18] and a value for A///(CIONO) given by Rohlack et 
al. [ 19). But as in the case of HONO, cir and trans forms of CIONO would 
be expected, and so, to a good approximation, the rata of Cl consumption 
by this channel would be doubled, that is, ■ 5.66 X 10" cm 6 /mole- 
cule 2 -s. The experimental rate constant thus represents the sum of the 
rate constants for the individual channels. The quoted value of ■ 0.28 
is assumed to be the same for each of the reactions. 


Table XIII: ClOO = Cl + 0 2 

The thermochemistry of ClOO obtained from the CODATA review of 
Baulch et al. [2], 


Table XIV: C10N0 2 = CIO + NO*- CIOONO = CIO + N0 2 

As indicated in the NASA review [1], experimental data for this associ- 
ation reaction fail into two sets that are in substantial disagreement. 
Several low-pressure measurements (20—231 of the rate of CIO disappear- 
ance in the presence of N0 2 are in good agreement and yield an average rate 
constant ho ■ 1.8 X 10* 31 cm 6 / molecule 2 *. No product identification was 
carried out and it was assumed that the reaction gave C10N0 2 . However, 
direct measurements [24,25] of the rate of thermal decomposition of 
ClONOj (/idiM * 6.2 X 10~ 23 cm 3 /molecule-s at 300 K) combined with the 
measured equilibrium constant gave k 0 ■ 4.5 X 10~ 32 cm 6 /molecule 2 -s. In 
Table XIV (a) the vibrational frequencies used for ClONO> are those given 
by Miller et al. [26], while the dissociation enthalpy AH!!(C10N0 2 > was 
taken from the study of Schonle et al. (25). If C10N0 2 is the sole product 
of the reaction, then to match the experimentally determined rate of CIO 
disappearance would require a large value of d r ■ 0.86. Chang et al. [3] 
have pointed out that this apparent disagreement can be resolved if it is 
assumed that an additional product. CIOONO. is formed in the reaction 
of CIO with N0 2 . The difference between the measured rate of CIO dis- 
appearance and the value of Schonle et al. for the rate constant for ClONOj 
production, that is. l.C X 10“ 3 ' cm'Vmolocule 2 *. is now identified as the 
rate of formation of CIOONO. The frequencies and thermochemistry of 
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C100N0 given in Table XIV(b) are those estimated by Chang et al. [3|. 
The Troe calculations are consistent with the suggestion that two isomers 
are formed since the necessary rate constants are predicted with reasonable 
values of 0 C . 

However, recent studies of the rate of formation of CIONO; from the CIO 
+ NO; reaction using diode laser spectroscopy [27) and of the Cl atom yield 
from photolysis of the products of this reaction (28) indicate that there is 
only one isomer of chlorine nitrate. This suggests that the most likely 
reason for the apparent discrepancy in the forward and reverse rate con- 
stants is the uncertainty in the thermochemistry of ClONOj. On the basis 
of the rate constant measurements, a “kinetic" equilibrium constant of 3.4 
X 10* molecules/cm 3 may be calculated, and this leads to A/f/lClONOj) 
■ 22.8 kJ/mol (cf. 26.4 kJ/mol as measured by Knauth [24]). In Table 
XIV(c) a calculation is presented based on this thermochemistry. The 
strong collision rate constant k £ is calculated to be 2.24 X 10“ 31 cm 6 / 
molecule 3 *, and so a ii r of 0.8 is required to match the experimental rate 
constant for loss of CIO. 


Table XV: BrGN0 2 = BrO + NO, 

The vibrational frequencies given for BrONOo Are those estimated by 
Sander et al. [29]. The thermochemistry of BrONOj is not known, but if 
it is assumed that AHlJlBrO-NOo) ~ A/f"(C10-N0i), the strong collision 
rate constant is then calculated to be 2.76 X 10 -31 cm 0 / molecule 1 ’* at 300 
K. Thus in order to reproduce the experimentally determined value of 5.0 
X 10 -31 cmVmolecule-* for BrO disappearance, a & t of 1.8 is required. As 
in the case of chlorine nitrate, this large value of d r would be explained in 
terms of a less stable isomer BrOONO, but there is no evidence for this 
species. Alternatively, the estimated thermochemistry of BrONO^ may 
be in error. It should be pointed out that Baulch et al. [2] recommend an 
estimated AH/tBrONO-j) ■ 20 i 30 kJ/mol * A/f/lClONO-..), but additivity 
arguments require that BrONO-. be substantially less stable than 

ClONOj. 


TabU XVI 

It should be noted that there is a large discrepancy between the values 
of A/Y/.:wm for FOO quoted by JANAF (12.55 k.J/mol) and the CODATA 
report (50 kJ/mol). For the purposes of this calculation, the JANAF value 
was used. If the heat of format ion were indeed as high as 50 kJ 'mol. then 
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AH(F - O 2 ) would be 28.91 kJ/mol and would be reduced to 1.14 X 
10 -33 cm 6 /molecule 2 -s, requiring a dc value of 14 to match the experimental 
rate constant. The results of the calculation presented here agree with 
those given by Jubert [30] who also used the JANAF thermochemistry. 


Table XVIII: FN0 2 F 4 VO*' FONO F + JV0 2 

Since there is some evidence suggest that less stable isomers may be a 
common occurrence among the halogen nitrites- calculations were per- 
formed for the corresponding fluorine analogs. The thermochemistry and 
spectroscopic parameters were estimated on the basis of similar molecules. 
The results indicate that the formation of the less stable species can easily 
be 2-3 times faster than that of the other isomer, largely as a result of its 
increased entropy. In the case of FONO, cis and trans forms need also be 
included. The possibility of isomers of fluorine compounds should 
therefore be noted by modelers and suggests a need for further experimental 
work in this area. 


Table XX: CH 3 0 2 ~ CH, + 0, 

The frequencies assigned to CH;|Oo were fixed so as to give an entropy 
consistent w*th that derived from group additivity The small value of d t . 
required to match the recommendation may be a consequence of incorrectly 
assuming F c orr = I. However, recent data [31] for this reaction suggest that 
it may be more complex at low pressures than believed at the time the ex- 
perimental data were evaluated. 


Conclusion 

It may be noted that for most of the reactions considered here, the value 
of dc required to give agreement between the calculated strong collision 
rate constant and the recommended experimental value lies in the range 
0.1-0.5. This is the range that might be considered reasonable on the basis 
of studies of energy transfer [32]. Some calculations require a dr value 
<0.1, which may indicate that the implicit neglect of coupling between 
rotation vibration involved in setting F corr may not be totally satisfactory. 
However, the general conclusion is that in the majority of cases this ap- 
proximate method of calculating third-order rate constants performs re- 
markably well. 
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ABSTRACT 

Experimental data concerning the title reactions are reviewed. An 
attempt was made to model the unusual pressure and temperature dependences of 
the rate constants within the framework of conventional RRKM theory. It was 
possible to fit all the experimental data with a single model but in order to 
do so, it was necessary to postulate a very looe^, highly entropic, H 2 04 
intermediate. Projections based on this model are made for the variation of 
the rate constant over wide ranges of temperature and pressure. 


*Thi8 work was supported, in part, by Contract DNA001-81-E-0146 (JPL 954815) 
with the National Aeronautics and Space Administration. 

March 1983; rev, May 1983 


Submitted to J. Phys. Chem. 


INTRODUCTION 


As a result of Its Importance to atmospheric and combustion chemistry the 
disproportionation reaction of gas phase HO 2 radicals has recently been 
studied by a number of groups. 

ho 2 + ho 2 -.i 2 o 2 + o 2 (1) 

Ibis work has indicated that the reaction rate shows marked pressure 1 J and 
temperature dependences* which are inconsistent with the original view 
that the reaction involved a simple metathesis in which one H0 2 radical simply 
abstracts a H atom from the other to give the products H 2 0 2 and 0 2 * Instead a 
mechanism involving the formation of an intermediate H 2 0 4 species is now 
favored: 


H0 2 + H0 2 t H 2 0 4 -*> H 2 0 2 + 0 2 (2) 

To date little attempt has been made to analyze this type of mechanism in 
detail to see if it can account for the observed kinetics or to make predic- 
tions of the rate constant at high pressures or temperatures. In addition to 
accurate rate data, it is necessary to have information on the spectroscopic 
and thermodynamic parameters of the H 2 0 4 intermediate in order to carry out 


2 


the required calculations. The available experimental data relevant to the 
problem are briefly reviewed In the following sections. Throughout this 
discussion, the rate constant for self-reactions, such as (1), Is defined by 



(a) Kinetic Data 

1 2 

Cox and Burrows and Thrush and Wilkinson'* first demonstrated that the 
rate constant for reaction (1) was pressure dependent. More recently Sander 
et al^ used flash photolysls/UV detection methods to show that the rate con- 
stant Is linearly dependent on pressure In a variety of diluent gases and 
varies from (1.91 ± 0.29) x 10" 12 CB 3 molecule” 1 s” 1 at 100 Torr to (2.97 
± 0.45) x 10“ 12 cm^ molecule”* s” 1 at 700 Torr In nitrogen. Extrapolation of 
their data to zero pressure suggests a pressure Independent contribution to 
the reaction with a rate constant of 1.75 x 10“ 12 cm^ molecule” 1 s” 1 . The 
form of the pressure dependence has been subsequently confirmed by Slmonaltis 
and Heicklen,^ who also used flash photolysis/UV-detection techniques. (There 
Is a small systematic difference In the rate constants reported by the two 
groups but this can be accounted for by the differences In the extinction 
coefficients assumed for the HO 2 radical.) The extrapolation leading to a 
zero pressure contribution to the rate constant has been supported by the low 
pressure studies of Thrush ond Tyndall®*® who reported a rate constant of 
(1.6 ± 0.1) x 10“ 12 cm^ molecule” 1 s” 1 between 7 and 20 Torr, using diode 
laser spectroscopy. 

The pressure dependence of the rate constant for 

D0 2 + D0 2 ->*D 2 0 2 + 0 2 (3) 


3 




t 
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c 


r. 


r 


r. 




has also been Investigated by Sander et al. They found the rate constant to 

be less dependent on pressure than the corresponding HO2 reaction, varying 

—13 3 —1 —1 

from an extrapolated 5 x 10 cur molecule L s at zero pressure to 
9 x 10 -13 cm -3 molecule -1 s -1 at 700 Torr of nitrogen. 

The temperature dependence of reaction 1 at pressures of 700 Torr or 
greater has been Investigated by several groups, and these studies all report 
strong negative temperature dependences. Cox and Burrows 1 generated HO2 
by the modulated photolysis of H2/O2/CI2 mixtures In a variety of diluents 
at a total pressure of 760 Torr, and they observed a temperature dependence of 
k^, between 275 and 339 K, which could be expressed as (3.8 ± 1.4) x 10~ 14 
expf(1250 ± 250) /T] cm 3 molecule -1 s -1 . Lll et al 4 used pulse radiolysis 
of H2/O2 mixtures at a total pressure of 2 atm and found a temperature 

—1 3 

dependence between 276 and 400 K which they expressed as (1.14 ± 0.16) x 10 
exp[(1057 ± 45) /T] cm -3 molecule -1 s -1 . Patrick and Pilling 3 used flash 
photolysis of CI2/O2/CH3OH mixtures at a total pressure of 700 ± 10 Torr to 
generate HO2 and found k to be given by (4.14 ± 1.15) x 10 -13 exp[ (630 ± 

115) /T] cm -3 molecule -1 s -1 between 298 and 510 K. Even at low pressure, the 

I 

temperature dependence was found to be strongly negative by Thrush and Tyndall** 
who reported a rate constant of 2.4 x 10 -13 exp(560/T) cm 3 molecule -1 s -1 
between 298 and 358 K at pressures of between 7 and 20 Torr. 

Q 

Very recently, Sander has also reported the temperature dependence of 
the HO2 + HO2 and DO2 + DO2 reactions, both determined by flash photolysl3/UV 
spectroscopy over the range 290-450 K. At zero pressure, the rate constants 
are given by: 

k t - 2.3 x 10 -13 exp(600/T) , k 3 - 2.2 x 10 -14 exp(900/T> 

4 



t 


while at 700 Torr, the blmolecular rate constants may be expressed as: 


k x - 1.8 x 10 " 13 exp(824/T) , k 3 - 1.7 x 10"* 4 exp(1201/T) 

The temperature dependence of k^ Is in quite good agreement with the previous 
data of Patrick and Pilling 3 and Thrush and Tyndall.® 

(b) Structure of H 2 O 4 

The molecular parameters of are well established and the values 

used In the subsequent calculation are given In Table 1. By contrast even the 
structure of the postulated H 2 O 4 Intermediate is open to speculation. At low 
temperatures, In matrix Isolation studies, a number of workers ** - * 4 have Iden- 
tified Infrared and Raman spectra of the products of electrical discharge pas- 
sed through water vapor, which they associate with a linear H0000H molecule. 
Only the skeletal motions of the 0^ chain have been identified In this way and 
there is substantial disagreement over the actual frequencies: 0-0 stretches 

are reported to be typically in the range 700-1000 cm - *, while 0-0-0 bends are 
in the range 400-600 cm - *. The actual values reported by Giguere and coworkers 
are given in Table 1 and these were used in some of the following calculations. 

Supportive evidence for the formation of tetroxides with 0^ chain struc- 
tures in the self reactions of peroxy radicals is provided by the observed 
products of the reaction of two CH 3 O 2 radicals :* 3 


2CH 3 0 2 

♦ HCH0 + CH 3 0H + 0 2 

(4a) 


♦ 2CH 3 0 + 0 2 

(4b) 


♦ ch 3 ooch 3 + o 2 

(4c) 






Table 1 MOLECULAR PARAMETERS FOR H 2 0 4 , D 2 0 4 , H0 2 , and D0 2 
® Vibrational Frequencies 


H 2°4 

w/cm * Type 

Ref. 

D 2 0 4 w/ cm 1 

v. 

3599 0-H stretch 

38 

2416 


• 

v 2 

3608 

0-H stretch 

38 


2029 


v 3 

855 

0-0 stretch 

11 


855 


V 4 

764 

0-0 stretch 

11 


764 


V, 

450 

000 bend 

11 


450 

* 

5 

V 6 

175 

H000 torsion 

(assumed) 


124 


V 7 

175 

H000 torsion 

(assumed) 


124 


v 8 

1402 

0-H bend 

38 


1015 


v ft 

1266 

0-H bend 

38 


966 


9 

v io 

823 

0-0 stretch 

11 


823 


V H 

430 

000 bend 

11 


430 


V 12 

100 

0000 torsion 

(assumed) 


100 

c. 


, -i 






ho 2 

w/cm 

Type 

Ref. 

do 2 



V 1 

3410 

0-H stretch 

10 

2411 



v 2 

1390 

bend 

10 

947 


r 

V 3 

1095 

0-0 stretch 

10 

1095 



Moments of inertia I/amu A“^ ^ 




c 

a 2°4 

ho 2 

°2°4 

CNI 

o 

Q 




16.6 

0.82 

18.3 

1.22 




113.2 

14.97 

120.0 

15.88 



t 

114.5 

15.84 

122.3 

17.04 




Lennard-Jones 

40 

Parameters 






e/k /K 

a /A 




t 

h 2°4 

577.0 

4.55 





N, 

71.4 

3.68 





t 
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Each of the sets of products Is consistent with the initial formation of a 
CH 3 O 4 CH 2 Intermediate. 

The HO 2 radical has a large dipole moment (2.1 Debye) , and it has been 
suggested 1 ^ that the intermediate in the disproportion reaction is in fact 
hydrogen bonded with a structure OOH-OOH or even a doubly hydrogen bonded 
cyclic species: 






o- rt 


(I) 


Lee and coworkers^ have used FTIR spectroscopy to investigate the products of 
the UV photolysis of a dilute matrix sample of t^CO or glyoxal in O 2 at 
12-20 K. They observed the IR peaks associated with RO 2 radicals and also 
other peaks at 3323, 1416 and 1108 cm - *, somewhat displaced from those of HO 2 , 
which they ascribe to a hydrogen-bondod HO 2 dimer. Since Lee and coworkers 
did not study the spectrum at longer wavelengths, it is not clear whether they 
were observing the same species as Glguere; their reported frequency shifts 
are consistent with those expected upon forming hydrogen bonds but, in addi- 
tion, they are not inconsistent with some of those frequencies expected for a 
linear HOOOH species. 

An experiment by Niki and coworkers has at least placed some restric- 
tions on the nature of the transition state Involved in the conversion of H 2 O 4 
to the final reaction products. Using FTIR detection they observed the iso- 
tope distribution in product H 2 O 2 formed from a mixture of and H^O^O 

radicals. Since a mixed product of the form H AO 0 lo 0H was never observed they 
were able to rule out (II) as a possible transition state. 




yO' '"C 

o' § l (III) 

ftt'n 0 


This experiment also excludes the possibility of a cyclic O 4 structure for 
H 2 O 4 but unfortunately, it cannot distinguish between a linear HOOOOH molecule 
leading to a transition state such as (III) or a hydrogen bonded intermediate. 

An additional unresolved question concerns the electronic state of the 
H 2 04 intermediate. Two H0 2 radicals in their ground state ( 2 A") can combine 
on either a triplet or a singlet surface but, if a linear HOOOOH Intermediate 
is involved, the triplet surface is likely to be repulsive. However, the 
assumed products H 2 0 2 and 0 2 ( 2 £g) can only be reached on the triplet surface. 

This raises the question of whether H 2 04 is a bound triplet state species, 
which may be possible for a hydrogen bonded species where the unpaired elec- 
trons would not be required for bonding, or whether the first formed products 
are the thermochemically accessible H 2 0 2 and 0 2 ( *A ) (Figure 1). 

O 

(c) Thermochemical Data 

1 9 

The heat formation of H0 2 has recently been redetermined by Howard who 
obtained it from a study of the H0 2 + NO t HO + N0 2 equilibrium. The reported 
value was AH^® ■ 10.5 ± 2.5 kJ mol“* . 

The heat of formation for H 2 O 4 depends on its structure and the nature of 
the bonding between the two H0 2 fragments. Two Independent investigations by 
Adamic et al . 20 and Bennet et al. 2 * - have been made into the low temperature 
(140 K) equilibrium in solution between t-butyl peroxy radicals and an assimed 


gun" o (II) 

• ^ 

Omtu'C 

A / 
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tetroxide species using ESR techniques. On the basis of these results, Nangla 
and Benson^ obtained a value for AH for the gas-phase reaction at 300 K: 

2H0 2 ♦ HOOOOH 

This was AH ■ -34.3 ± 4.2 kJ mol - *. Benson and Shaw^ estimated the corres- 
ponding entropy change on the basis of model structures as As ■ -143.1 J mol -1 
K -1 . The correction of the measured enthalpy change at 140 K in the liquid 
phase to the gas phase at 300 K was recalculated here (see appendix) and the 
corresponding entropy change was also derived. The values found were AH ■ 
-58.6 kJ mol”* and AS ■ - 237.7 J mol”* K”*. The origin of the difference 
between this value for AH and that reported by Nangla and Benson lies in the 
different assumptions made in estimating heats of vaporization. This shows 
that there is a large uncertainty in deriving the thermochemlstrv of gas phase 
HOOOOH. 

It is unlikely that a singly hydrogen bonded intermediate would be more 
strongly bound than Benson's estimate for HOOOOH, Indeed a value of about 
21 kJ mol”* is often suggested. This may presumably be doubled in the case of 
the cyclic species. 
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REACTION MECHANISM 


The form of Che pressure and temperature dependences for reaction (1) 
seem to be well established. 

The potential energy surface over which the reaction occurs may be con- 
sidered In terms of two limiting case.. Illustrated In Figure 1, both of which 
could, In principle, account for the observed pressure dependence of the rate 
constant, k^. 

Type (a) : Here two essentially Independent channels are operative (5,6). 

The zero pressure rate constant of 1.75 x 10“^ cm^ molecule s - ^ represents an 
abstraction reaction leading directly to products. The pressure dependence of 
the experimentally determined rate constant, k Qbg , is represented by a pressure- 
dependent approach to equilibrium between HO 2 and H 2 O 4 , the only loss process for 
H 2 0 ^ being redlsaoclatlon: 

H0 2 + H0 2 ♦ H 2 0 2 + 0 2 (5) 

[M] 

H0 2 + H0 2 t H 2 0 4 ( 6 , - 6 ) 

If H 2 0^ redissociates rapidly, then an equilibrium Is set up and the observe 
rate constant for loss of H0 2 Is simply kj. On the other hand, If k_g Is very 
slow, the H 2 O 4 is essentially stable and k obg is given by kj + kg[M]. In the 
Intermediate region, at high pressures, substantial amounts of the H 2 0 ^ complex 
will be formed and the observed rate constant will be given by 


Under euch conditions, k Qba Is not a true rate constant since It will be time 
dependent, but If k_g Is sufficiently small, then this would not be apparent on 
the timescale of experiments, such as those of Sander et al. , l.e., ~ 20 ms. 

Similarly, deviations from second-order kinetics In the observed loss of HO 2 
radicals would be small. Reaction (6) may be assumed to be at Its low-pressure 
limit In view of the weak bonding In H 2 0g, and Its relatively low density of 
states If the frequencies In Table 1 are used. The rate constants kg and k..g can 
then be calculated using approximate methods developed by Troe. Taking kg 

to be given by the experimentally determined 1.75 x 10 cm molecule a , it 
Is then possible to simulate the time dependence of H0 2 over a timescale of 20 ms 
using numerical Integration techniques, and hence extract the apparent second 
order rate conatant for loss of H0 2 as a function of pressure. Adopting this 
procedure. It was found that if Benson's suggested bond strength of 34.3 kJ la 
assumed for a HOOOOH species with the parameters given In Table 1, then the life- 
time of the HOOOOH molecule la too short to affect the observed H0 2 kinetics and 
k Q | }g 8 Imply equals kg, which la Independent of pressure. In order to Increase 
the lifetime of HOOOH sufficiently to account for the experimental pressure 
dependence In k Qbg , It la necessary to raise the bond strength to around 77 kJ 
mol - *'. 

An additional problem for this model, even If a stronger bond Is assumed, is 
presented by the temperature dependence. Ar the temperature Is raised, the life- 
time of such a HOOOOH Intermediate decreases rapidly until by 4C0 K the only loss 
of H0 2 la once again by reaction (5). Thus at temperatures greater than 300 K, 
the temperature dependence of k 0 ^ fl would be essentially that of kg, and It Is 
unlikely that the rate constant for a disproportionation reaction would have such 
a strong negative temperature dependence as that observed experimentally. 
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The laoe behavior would be expected If hydrogen bonded HOOHOO Intermediate 
were assumed, and In this caao It would probably be more difficult to Justify the 
increase In bond strength. 

Type (b): In this model It Is assumed that there Is no direct dispropor- 

tionation reaction at all. Instead, the reaction products H2O2 end O2 are formed 
by a molecular elimination from an HOOOOH Intermediate. The mechanism may be 
represented by: 

k. 


ho 2 + ho 2 


h 2 o 4 


* k r 


H 2 02 + O 2 


-a 


k.[M] *+ k_ f [M] 


h 2 o 4 


( 8 ) 


From this, k Qbg may be derived as 


k r k_ s [M][H 2 0 4 l/[H0 2 r 


( 9 ) 


obs k [M] + k +“k 

s J -a r 

where the rate constants kj, k r , k_ , k 8 , and k_ g are defined above. The zero 
pressure contribution arises from those vlbratlonally excited H 2 0 4 A molecules 
which rearrange to give H 2 02 and 0 2 directly, while the pressure dependence Is 
a result of the stabilization of some H 2 0 4 * Into the well. The stablized H 2 0 4 
molecules may or may not subsequently decompose thermally depending upon the 
depth that Is assumed for the well; the thermal decomposition products can be 
H0 2 radicals or R 2 02 + 0 2 , depending upon the relative barrier heights and 
A-factors. 

This kind of model Is more difficult to test numerically than type (a) 
and a two-channel RRIQf model was developed In order to make the necessary 
simulations. The calculation method In described below. 
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CALCULATIONAL METHOD 


H2O4* Is formed by a chemical activation process with a rate constant 

m 

In highly vlbratlonally excited states from which It may be stabilized 
by collision, decompose back to reactants with a rate constant k_ a (E) or 
eliminate H2O2 and 0 2 with a rate constant k r (E). If the thermallzed 
does not decompose on the experimental timescale, then the rate constant for 
the loss of HO2 radicals Is given by 

- k(E)+PZ[M] 

k obs \ •^e 1 k_ a (E) + k r (E) + PZ[M] f(E)dE (10) 

This expression Is derived In the same way as In the usual strong-collider 
RRKM theory, where a steady-state treatment for ^O^CE) Is used. E^ is the 

<9 

barrier to redissoclatlon, k^ is the high-pressure limiting rate constant for 

reaction (1), Z Is the Lennard-Jones collision frequency, P Is the colllsional 

2 9 

deactivation efficiency, and f(E) Is the chemical activation distribution 

k (E)N(E) exp(-E/kT) 

f(E) - -s (ID 

Jg k_ a (E)N(E) exp(-E/kT) dE 

where N(E) is the density of states of tkO^ at energy E. 

However, If the bonding in HjO^ is sufficiently weak, molecules deacti- 
vated Into the well may decompose either back to HO2 radicals or to products 
on the experimental timescale. Monte Carlo calculations 0 showed that the 
chemical activation process In the subsequent models at pressures of 100 Torr 
and greater leads to the rapid (< 5 ns) thermalizatlon of those H2O4 molecules 
that survive until the first collision occurs . The secondary (thermal) 
decomposition of stabilized H2O4 may therefore be treated by a conventional 


two-channel thermal decomposition mechanism. The rate constants for the 

27 

decomposition of thermallzed molecules are given by: 


c -a " SI 


k (E) PZ[M] 


E x k_ a (E) + k r (E) +PZ[M] 


f ' (E)dE 


( 12 ) 


and 

k f (E) 8 Z[M] 

k r " k (E) + k (E) + 0 Z [M] f ' (E)dE 

i —a r 


where f'(E) is the thermal distribution function: 


( 13 ) 


fl(E) _ N(E) exp(-E/kT) (14) 

Q vlb 

and E2 Is the activation barrier for the molecular elimination channel. 

The probability, p, of thermal decomposition leading to H2O2 + O2 rather 
than back to HO2 radicals is then given by: 


P * 



( 15 ) 


-a r 

If the thermal reaction is fast compared to the experimental timescale, then 
the overall rate constant for loss of HO2 nay be expressed in a modified 
version of ( 8 ): 


„ k(E)+PZ[M]p 

k obs " ^ E. k_ (E) + k (E) + PZ[M] f(E)dE (16) 

1 r 

The microscopic rate constants k__ (E) and k_ (E) were calculated using 
RRKM theory and assuming fixed transition states. Exact state counting^® was 
used to determine the sums of states for the activated complexes and the 
Whitten-Rabinovitch ^ 7 approximation was used to obtain the density of states 
of the molecule. Anharmonic corrections to the densities of states were made 
according to the Haarhof : s procedure. 
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The value of the high pressure rate constant k^, for this reaction is not 

known experimentally, and for the purposes of this calculation it was assumed 

equal to 1.4 x 10 -11 cm 3 molecule - * s -1 by analogy with the value for the 

31 

recombination of the isoelectronic ethyl radical. This may be an underesti- 

mate since using the Isoelectronic argument predicts that the high pressure 

—1 1 3 

rate for OH + HO2 would be the same as CH3 + C2H5, i.e., ~ 2.5 x 10 cm 
molecule -1 s -1 , (obtained by the cross-combination rule based on rate con- 
stants for C2H5 + C2H5 31 and CH3 + C^ 3 ^), whereas at 1 atm DeMore 33 reports 
the rate constant to be (1.2 ± 0.4) x 10 -1 ^ cm 3 molecule -1 s -1 and pressure 
dependent. The high-pressure rate constant for the DO2 + DO2 reaction was 
also taken as 1.4 x 10 -11 cm 3 molecule -1 s -1 , and both high-pressure rate 

constants were assumed to be Independent of temperature. 

2 S 

Following Troe the collision rate constant Z was set equal to the 
Lennard-Jones collision rate constant given by: 

Z - N a a^ M (8 RTMi a _ m ) 1/2 Q 2 » 2 (A-M) (17) 


where N. is the Avogadro number, p . M is the reduced molar mass, and o is 
A A*"M A"M 

the average Lennard-Jones collision diameter given by (o + o ) /2. The 

A M 

2 2 

reduced collision integral Q * (A-M) may be expressed as an approximate func- 
tion of kT / c a.-M* w ^ ere e A-M is tlie Lennard-Jones well-depth defined by 


/ xl/2 

£ A-M “ (£ AA e MM ) 


2 2 

Q * (A-M) 


[0.697 + 0.5185 log(kT/e A u )] 

A - n 


(18) 


The relevant Lennard-Jones parameters are given in Table 1. The temperature 

dependence of P, the deactivation efficiency, was described by the approximate 

25 

equation of Troe: 
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p 

~ 71 


<£E> 


(19) 


F E kT 

where < AE >, the average energy removed by collision, was taken to be 1.88 kJ 
mol - ^ and Independent of temperature, and Fg Is the energy dependence of the 
density of states. 

This treatment of "weak” collision effects Is perhaps not as realistic as 

the Monte Carlo calculations Incorporating an exponential model for energy 

2 & 

transfer, as described by Barker. However, a comparison was made, for a 
simple chemical activation case, between the results of a stochastic calcula- 
tion with starting energies selected from the chemical activation distribution 
function (equation 11) by Monte Carlo techniques and the results obtained by 
numerical integration of equation (10). These results were In complete agree- 
ment, which justifies the use of the simplified treatment of the energy trans- 
fer process (equation 10). The direct integration of equations (10) and (16) 
was found to be more efficient computationally than the equivalent Monte Carlo 
technique, which required large numbers of trajectories and long running times 
for systems In which the subsequent decomposition of thermalized molecules was 
important. 

The transition state parameters used In the first calculations, assuming 
a linear H0000H structure for the adduct are given in Table 2. The frequen- 
cies and structure of the transition state leading back to HO 2 fragments, 

(1)» were chosen so as to be consistent with the high pressure limit 

00 

rate constant. In order to ensure that was temperature independent some 
of the modes in this transition state were systematically tightened as the 
temperature Increased. The structure of the second transition state, 

(2), was assumed to be like (III) and the frequencies chosen were based on 
those of hydrogen peroxide with additional modes assigned according to 
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Benson's methods.^ The structure and frequencies of the corresponding deut- 

erated species were based on the original protonated ones and chosen so as to 

29 

be consistent with the Teller-Redlich product rule. 

Calculations were performed for three kinds of linear HOOOOH molecule in 
an attempt to account for the observed pressure dependences of the HO 2 + H0 2 
and D0 2 + D0 2 reaction and the temperature dependence of the H0 2 + H0 2 
reaction. The first of the species had frequencies as reported by Glguere 
and a dissociation energy AH(H0 2 - 0 2 H) as given by Benson, the second had 
the same frequencies but was more strongly bound and the third had looser fre- 
quencies and a weak bond. 


Table 2. TRANSITION STATE FREQUENCIES 


H 2 0 4 t (l) 
w/cm * 

Type 

H 2 0 4 f (2) 

, -1 
Gj/cm 

3500(2) 

0-H stretch 

3774 

1390(2) 

0^0 bend 

3788 

1000(2) 

0-0 stretch 

1266 

60(2) pseudo vibrations 

1380 

45(2) 

correlating with 

880 

25 

translation 

1580 


and rotation 

300 



300 



225 



150 



150 


D 2 0 4 t (l) D 2 04 f (2) 


Type <i)/cm * oj/cm 


0-H stretch 

2475(2) 

2800 

0-H stretch 

980(2) 

2800 

H*0 - 1 0 bend 

1000(2) 

1025 

00H bend 

60(2) 

1025 

0-0 stretch 

23(2) 

880 

0 — — 0 stretch 

32 

1580 

0*0 stretch 


300 

0.0 — — 0 bend 


300 

torsion 


175 

0*H*0 bend 


113 

0*0 *H bend 


113 
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A potential surface intermediate between (a) and (b) was also considered 
where some, but not all, of the experimentally observed zero-pressure rate 
constant may be ascribed to a direct disproportionation reaction. 

The additional uncertainties associated with possible hydrogen-bonded 
intermediates precluded any detailed calculations, but some conclusions can 
be drawn on the basis of the results of the calculations involving HOOOOH 
intermediates. 
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RESULTS 


(i) Tight molecule, weak bond: Using the frequencies and bond strength 

In Table 1 for the HOOOOH molecule and the transition state parameters given 
In Table 2, It was possible to obtain the observed zero pressure rate constant 
If the barrier to the decomposition of HOOOOH to products, E 2 , was set at 10.0 
kJ mol - *. However, no pressure dependence Is predicted when the pressure Is 
raised to 700 Torr. Within this model the microscopic rate constants for back 
dissociation, k_ a (E), are too fast for competition by collislonal stabili- 
zation to occur. It was found that the zero pressure rate constant is 
determined primarily by the difference In energy between the two barriers, E^ 
and E 2 , for the decomposition of the intermediate molecule, while the pressure 
dependence Is determined by the density of states of the molecule at energy 
Ej_. In order to obtain the observed experimental pressure dependence, it Is 
necessary to either increase the well depth over the value suggested by Benson 
or increase the entropy of the molecule by reducing some of its vibrational 
frequencies from the values reported by Glguere. These options were investi- 
gated in turn. 

(11) Tight molecule, strong bond: If E^ is increased to 140.0 kJ mol - * 
and E 2 is set at 116.0 kJ mol - *, then it Is possible to obtain the experi- 
mental values of k^ at zero and higher pressures, as shown in Figure 2, using 
equation (10). With this deep a well, the lifetime of stabilized HOOOOH is 
~ 2 x 10^ ^ for unimolecular decomposition at 700 Torr, which is much longer 


than the experimental timescale. However, It Is known that H 2 O 3 has a very 
similar absorption spectrum to and so it may be expected that the 

spectrum of H 2 O 4 would also resemble that of H 2 O 2 . Thus the residual absorp- 
tion observed after the disappearance of the HO 2 radicals In the flash photol- 
ysls experiments might have been due to H 2 O 4 rather than to 

The predicted temperature dependence for this model Is Illustrated In 
Figure 3 along with available experimental results. The prediction Is based 
only on the room temperature data of Sander et al ., 3 and there are small 
systematic differences between that data set and the other temperature- 
dependent data. Variation of the frequencies for ^0^(2) Indicated that any 
loosening of the transition state would change the temperature dependence from 
negative to positive which would then be In disagreement with the experimental 
data. 

The remaining data which have to be taken into account within the model 
are those for the pressure and temperature dependence of the DO 2 + DO 2 
reaction. By making a suitable choice of frequencies, consistent with the 
Teller- tedlich relation, for the transition state for decomposition to pro- 
ducts it Is possible to arrive at values of E^ and E 2 , based on the protonated 
case, but corrected for zero point energy differences, which lead to the 
experimental value for kj at zero pressure. E^ Is then 140.8 kJ mol -1 ' and E 2 
is 128.8 kJ mol” 1 . Thus although Ej^ is not affected greatly by deuteration, 

E 2 Is substantially raised when deuterium is substituted for hydrogen. How- 
ever, it can be seen from Figure 2 that the p-essure-dependence predicted by 

—12 3 —1 

this model is too steep, with k^ predicted to be 2.9 x 10 cur molecule 
s ” 1 at 700 Torr. The enhanced pressure dependence arises because the density 
of states of D0000D at E^ is greater than that of HOOOOH, and this reduces the 
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redissociation rate of the first formed DOOOOD* and enables colli9ional sta- 
bilization to compete more effectively.. Within this model, the only way in 
which the pressure dependence of the DO 2 + DO 2 reaction could be less than 
that for HO 2 + HO 2 would be if the well were sufficiently shallow that the 
stabilized molecules also decompose on the experimental timescale. The 
product distribution of this later decomposition would be determined largely 
by the relative barrier heights and since E 2 is higher for the deuterated 
case, back decomposition might be more favored than in the corresponding 
protonated case. This "loose molecule-weak bond" model is Investigated in the 
next section. 

(ill) Loose molecule, weak bond: If the value of is taken to be that 

suggested by Benson, then it is necessary to reduce some of the frequencies of 
the HOOOOH molecule from the values reported by Giguere in order to match the 
experimental pressure dependence (Table 3). In addition, as E^ is reduced, it 
is necessary to take account of the thermal decomposition of stabilized HOOOOH 
by using equation (16) to calculate k^. The lifetime of thermallzed "loose" 
HOOOOH is ~ 2 x 10”^ s at 700 Torr and 300 K. With this model it is possible 
to reproduce the experimental pressure dependence for the HO 2 + HO 2 reaction 
with E^ ■ 34.3 kJ mol” 1 and E 2 ■ 10.0 kJ mol” 1 using the molecular parameters 
in Table 3 and the previous transition state parameters given in Table 2. It 
can be seen that it was necessary to Increase the entropy of the molecule 
quite drastically with the 0-0 stretching frequencies reduced to 100 cm” 1 and 
the 0-0-0 bending frequencies reduced to 95 cm” 1 . The realism of this step is 
open to question. The predicted temperature dependence for this model at zero 
pressure and 700 Torr is similar to that calculated in (il) and presented in 
Figure 3. It can be seen thac the temperature dependence of the rate constant 


is somewhat underestimated by these calculations. The agreement can be 

CO 

Improved If it Is assumed that k^ has some negative temperature dependence and 
the result for the case when 1^1 T ‘ Is also illustrated In Figure 3. The 

CO 

true temperature dependence of k^ Is, of course, unknown but by analogy with 
alkyl radical recombination^ 1 Is might be expected to be slightly negative. 

As In (11), using the same frequencies f or the transition state leading 
to products, D 20 ^( 2 ), It Is possible to fit the zero-pressure rate constant 
for the DO 2 + DO 2 reaction (Table 3). The corresponding values of E^ and 


Table 3. MOLECULAR FREQUENCIES FOR -LOOSE" HOOOOH MODEL 
H2O4 

. -1 , -1 
w/cm u/cm 


3600(2) 

2546(2) 

1500(2) 

1060(2) 

100(3) 

100(3) 

95(3) 

95(2) 

80(2) 

72(2) 

60 

42 


of E^ and E 2 are 35.1 and 22.8 kJ mol -1 , respectively. When the pressure is 
raised to 700 Torr, kj Increases to 1.2 x 10“*^ cnr* molecule”* s”*, in good 
agreement with the experimental value. The pressure dependence of k 3 is 


moderated by the thermal reaction and so, in this way, It Is possible for the 
pressure dependence of k^ to be less than that of k^. 

The predicted temperature dependence of kj at zero pressure and 700 Torr 

Q 

Is given In Figure 4 along with the recent data of Sander. This prediction 

Is based only on the room temperature pressure dependence for the HO2 + HO2 

reaction and the assumption that k^ Independent of temperature, and not 

optimized for the DO2 + DO2 reaction, but even so the agreement Is quite 

■ _o 

good. The results of a calculation In which k^ Is arbitrarily assigned a T 
dependence are also Illustrated. 

It Is possible to obtain the same agreement with experiment for Increased 
values of E^, up to ~ 60 kJ mol” 1 If the frequencies of the HOOOOH molecule 
are raised accordingly. For well depths greater than 60 kJ mol”*, the life- 
time of H2O4 becomes too long for thermal decomposition of the stabilized 
species to play a role, and the predicted pressure dependence for reaction (3) 
Is then too great to give agreement with experiment, as In (1). Molecular 
parameters which give satisfactory fits to all the experimental data are given 
in Table 4, along with those for a model based on Glguere's frequencies which 
fits only the HO2 data. 


TABLE 4. MOLECULAR PARAMETERS FOR HOOOOH USED IN FITTING KINETIC DATA 


/kJ mol” 1 

E2/kJ mol” 1 


“nnn /c »~ l 

S/J mol” 1 

34.3 

10.0 

100 

95 

379.5 

40.0 

15.7 

145 

100 

369.7 

60.0 

35.7 

300 

150 

346.5 

140.0 

116.0 

Glguere 

Model 

298.9 


Figure 5 shows an extended promotion of the temperature dependence of 
reaction (1) from stratospherlcally Important temperatures to those more rele- 
vant to combustion processes. Above ~ 1000 K, it can be seen that the reac- 
tion rate is predicted to be Independent of pressure and begins to Increase 
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with temperature, with a rate constant of ~ 3 x 10 cnr molecule s at 
3000 K. The same temperature dependence Is obtained for all values of up 
to 60 kJ mol - *. Figure 6 examines the pressure dependence of k^ at 300 K up 
to pressures near the high pressure limit. In the loose molecule-weak bond 

m 

case, It can be seen that the high pressure limit for Is not k^ but 

Ol • GO 

p k^ , where p Is given by equation (IS), when k r and k_ a have their limiting 
high pressure values. However, In this Instance, p is ~ 0.94, and so the 
results for models (11) and (111) are essentially the same. Below 1 atm the 
pressure dependence is linear, but curvature becomes apparent at only a 
few atmospheres, and this pressure regime should be easily accessible 
experimentally. 

In addition to type (a) and (b) potential surfcces, It Is possible to 

envisage an Intermediate case where some, but not all, of the zero-pressure 

rate constant may be ascribed to a direct disproportionation reaction. This 

possibility is more difficult to assess since the magnitude of this direct 

contribution to the zero-pressure rate constant, k d , not known. However, if 

- 1 3 

kj Is assumed to be temperature-independent, It must be less than 8 x 10 
cnr* molecule - * s - *, the slowest experimentally observed value for the zero- 
pressure rate constant at 450 K. Taking k d ■ 7 x 10 cnr molecule s , 
and assuming a 60 kJ mol - * bond for H0000H, it was possible to fit the experi- 
mental pressure dependence of the HO2 + H0 2 reaction with 0-0 strecches of 300 
cm - *, 0 0 0 bends of 180 cm - * and E 2 ■ 40.0 kJ mol - *. 


The temperature 


dependence of Che rate constant at 700 Torr was calculated to be the sane aa 
the case where kj ■ 0, however, the temperature dependence at zero pressure is 
somewhat less marked than In the latter case. Thus, it appears that it is not 
possible to rule out the possibility that some contribution to the total rate 
constant Is due to a direct disproportionation reaction. 

(lv) Hydrogen-bonded Intermediates: Due to the additional uncertainty 

associated with this kind of molecule, extensive detailed calculations were 
not attempted. However, on the basis of the previous results and some limited 
calculations on this model, It is possible to draw some useful conclusions. 

For a linear hydrogen bonded species, H00H00, the binding between the 
radical fragments Is unlikely to be stronger than the 34.3 kJ mol -1 assumed by 
Benson for HOOOOH, Indeed It may be substantially less. Thus, In order to 
achieve a high enough density of states to account for the pressure dependence 
of the reaction, It too would be required to have a number of very low fre- 
quencies. In addition, the transition state for product formation will be 
much looser than that for the previous model and will resemble more closely 
that for back decomposition to reactants. Consequently, the values of k r (E) 
will be much faster, and It will be necessary to raise E 2 relative to E^ In 
order to obtain the correct zero pressure rate constant. Depending on the 
exact nature of this transition state. It may even be necessary for E 2 Co he 
greater than E^ , In which case the temperature dependence of the zero pressure 
rate constant will no longer be negative, and hence will he at variance with 
experimental results. 

It has been suggested 16 that a doubly hydrogen bonded Intermediate may 
be formed. However, these two bonds would have to be formed consecutively 
and, given the weakness of the first hydrogen bond and the lower A-factor for 
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cycllzatlon compared Co redlssoclatlon, this process is not likely Co be 
favored. Even If cycllzatlon should occur, Che concomitant reduction In the 
density of states at an albeit greater value of Ej would make It difficult to 
accommodate the observed pressure dependence. 
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CONCLUSION 


The calculations presented here indicate that any model based on a type 
(b) potential surface is quite heavily constrained in terms of the molecular 
and transition state properties, required to reproduce the experimental 
pressure and temperature dependences of the HO2 + HO2 and DO2 + DO2 reactions. 
The observation that the pressure dependence of the DO2 + DO2 reaction is less 
than that of the HOj + HO2 reaction, requires that the interfragment bonding 
be weak ( 30-60 kJ mol“^), but in this case, the H2O4 molecule must be allowed 
a higher entropy chan Giguere's frequency assignment permits in order for any 
pressure dependence to be predicted at all. In addition, in order to get a 
negative temperature dependence for either reaction, the transition state 
leading to products must be made very tight. Even so, the temperature depen- 

GO 

dence is underestimated unless k^ is also assumed to have a negative temper- 
ature dependence. This problem cannot be overcome by further tightening the 
transition state since it is then becomes necessary to reduce E2 in order to 
reproduce the zero pressure rate constant, and under these conditions the bar- 
rier leading to products is so low that the rate constant is no longer pres- 
sure dependent. 

It is not clear from the currenc literature whether a very loo^e struc- 
ture for H2O4 is possible or not. The frequencies assigned by Giguere and 
coworkers to H 0000 H observed in a matrix are certainly consistent with those 
that might be expected on the basis of analogous molecules, such as H2O2 or 
H2S4, assuming similar bonding. However, the bond strengths in these stable 
molecules are substantially greater than that in any of the H2O4 species 
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listed in Table 4 (e.g., AH(HO-OH) » 215 kJ mol”*), and arguments may be made 
based on Badger's rule to support the suggestion that the 0-0 frequencies 
should be less than those in ^ 02 * Also, in addition to possible problems 
Involved in assigning frequencies in a complex system such as that analyzed by 
Giguere, it is not clear that the intermediates they observed at low tempera- 
tures in a matrix are necessarily the same as that formed in reaction ( 1 ) in 
the gas phase. 

It is also difficult to find any potential problems with the experimental 

data on the kinetics of reaction (1). All the rate data which display a 

pressure dependence were obtained using UV spectrometry to detect the HO 2 

radical, but the reported gas phase and solution spectra for HO 2 are very 

similar, and so it is unlikely that any unexpected pressure dependence in 

the spectrum is possible. However, confirmation of the pressure dependence of 

the rate constant by some other technique would be valuable. It would also be 

useful if more kin-tic data were available at very high pressures or in non-polar 

solvents to validate the pressure dependence and provide a measure of the high- 

pressure limit rate constant. In aqueous media, the HC >2 + HO 2 reaction is very 

slow, 1.26 x 10- J cm molecule s , presumably due to solvent effects. By 

37 

contrast, studies by Howard and Ingold suggest that in tetrachloromethane the 
rate constant may be as fast as 9 x 10 cm molecule s , although this value 
was determined in an indirect experiment. It is always necessary to use fairly 
complex chemical systems to generate HO 2 radicals, but sufficient cross-checks 
seem to have been made to confirm the expected chemistry. Even so, it is some- 
what unexpected that the ratio k^/k^ is independent of pressure. 

In view of the expected weak bonding in H 2 O 4 , and because when it is first 
formed, energy is already in the reaction coordinate for redlssoclation, it might 


23 



be suggested that a statistical RRKh calculation of k_ fl (E) may be Inappropri- 
ate. But if the statistical assumption does bre&ic down, k_ fi (E) would be faster 
than calculated by RRKM techniques, and it would be correspondingly more diffi- 
cult to account for the experimental pressure dependence. 

Further experimental work is required in order to clarify the mechanism of 
this reaction, although it is unlikely that any single experiment will yield 

3D 

sufficient information. In addition to offering an opportunity to measure k^, 
studies in solution may lead to Improved spectroscopic characterization of H2O4. 
It may also prove fruitful to investigate the electronic state of the oxygen 
molecule produced in this reaction. If ©2 ( ^ ) is formed, then following the time 
dependence of its production may lead to information on the lifetime of H2O4. 

The main conclusion of this study is that either H2O4 indeed has a very 
loose structure or that some other mechanism, rather than one involving associa- 
tive comple'-: formation, is required to account for the experimental data. 
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Appendix 


BOND DISSOCIATION ENTHALPY FOR tBu0 2 -0 2 tBu 

Following Nangia and Benson, ^ the experimentally determined thermoche istry 
for the equilibrium: 


tBu0 2 ~0 2 tBu * 2 tBu0 2 

20 21 

in solution at 140 K is taken as a starting point: ' 

AH 1 140 " 36,8 ± 4,2 ^ mol_1 

AS* 140 * 142,3 * 25,1 J 1,101-1 K_1 

The superscript (+) refers to measurements in liquid-phase systems (hydrocarbon 

solvents) at 140 K with the standard state of 1 Molar. 

In order to translate these data to thermodynamic quantities for the gas 

phase at 300 K, it is first necessary to "correct" the liquid phase data to 300 

K, and then to transpose them to the gas phase. 

The appropriate quantity required to correct the enthalpy and entropy to 300 

K is the liquid-phase heat capacity as a function of temperature. To approximate 

this quantity estimates of the heat capacities of 2,2,7,7-tetramethyloctane and 

2,2-dimethyl butane were used as surrogates for ditertiary butyl tetroxide and 

the tertiary butyl peroxy radical, respectively. Using the group additivity 

methods of Luria and Benson, 141 AC is found to vary *rom 25.5 J mol -1 K -1 at 140 

P 

30 






K to 7.8 J mol - * K - * at 300 K. In view of the other approximations involved, it 
is reasonable to use the average value, AC p ~ 18.8 J mol - * K - *, at all temper- 
atures. Thus the liquid-phase parameters at 300 K may be calculated as: 


AhJ 300 - AhJ 140 + AC p (300-140) 10 -3 - 39.7 kJ mol -1 

iS I,300 - AS I,140 +iC p ln(30 ° /l40> 


156.5 J mol - * K - * 


In order to convert these values to the corresponding gas-phase quantities, 
the methods of Mendenhall and Benson^ were used which lead to: 


AH 


° _ * AH. _ + 2 ah (tBu0 o ) - AH (0. )+ AC [T - T - 2T_] 
g,T Jt,T vap 2 vap 4' vap 1 o B J 


AS 


g,T 


- Ast _ + Rln(V ) + TC + AC In ( 


T T 


*,T 


vap 


) 


where the superscript (o) indicates the 1 atm standard state, Tg is the boiling 

T B 

point of tBu02» T Q is the boiling point of tBuO^tBu, AH v (tBu02) and 

T Vap 

AH ° (0,) are the heats of vaporization at the boiling points, V 0 is the molar 
van 4 “ 


vap 

volume of the solvent, and TC is Trouton' s constant. 


It is observed empirically^**^ 3 that AC ~ 50.2 J mol - * K - * and so it was 
T vap 

estimated that AH ° (0, ) ■ 41.8 kJ mol - * with T rt ■ 467 K (based on an average of 
vap 4 o 

values for dodecane and di-t-amyl peroxide) and AH vap (tBu02) ■ 27.2 kJ mol -i with 
Tg ■ 323 K (based on the values for 2,2-dimethyl butane). Rln(V g ) ~ 20.9 J mol * 
K - * for many solvents, and so 


AH g,300 


A %,300 


+ 18.8 ■ 58.6 kJ mol - * 


AS' 


ASl 


+ 83.7 - 238.5 J mol“* K -1 


The value for AH 

.22 


g,300 


g,300 A ,300 
is almost a factor of 2 larger than that deduced by Naiigia 


and Benson^ 4 from the same experimental data. This discrepancy arises from the 
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use of different hydrocarbon species as models for tBuO- and tB^O^ and different 
methods of correcting ^ H vap fro® the boiling point of the liquids to 300 K. This 
Indicates that the real uncertainty In the bond dissociation enthalpy at 300 K In 
the gas phase Is at least ± 24 kJ mol -1 . 

The value for AS° ■ 238.5 J mol -1 K -1 Is substantially greater than the 
value of 150.6 J mol -1 K -1 obtained by additivity methods where group values were 
estimated assuming a "normal” structure and frequencies. On the basis of this 
value of AS° the entropy of ^0^ would be only 217.6 J mol -1 K -1 which 

implies a tighter molecule than even that suggested by Giguere (Table 4). 

If, on the other hand, the additivity value of AS° ^qq ^ 8 assumed correct 
and converted to AS^ ^q, then a value of 220.0 J mol -1 is obtained. This may be 
combined with the experimental value of AG^ ■ 16.7 kJ mol -1 to yield 


AH 


■47.7 kJ mol -1 . However, these parameters are not consistent with the 


i ,140 

reported temperature dependence of the equilibrium constant. It thus appears 
that while a looser ^0^ species than expected from additivity methods is 
required by the kinetic data, a tighter one is needed to fit the equilibrium 
data. However, it should be borne in mind that the uncertainties in this kind of 
analysis are large and so the reliability of the gas-phase thermodynamic 


parameters is difficult to assess 
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1. INTRODUCTION 

A complex chemical system Is descrlbable In tens of discrete elementary 
steps. The rate constants for these steps are the Input coefx’lclonts to the 
ensuing differential equations which may be written to describe the time- 
dependent change of concentrations. Solutions of the mathematical problem 
represented by these equations have physical meaning onxy to the extent 
that the appropriate set of steps has been taken into account (mechanism) and are 
quantitatively correct. Furthermore, the uncertainty In the solution depends on 
the combination of uncertalrty and sensitivity of these rate constants. 

We will discuss herein the current best framework of theory for evaluating 
and estimating gas-phase rate constants. We also will discuss some experimental 
methods by which rate constants are measured. This combination of measurement 
and the framework provided by theory Is the basis for our very pragmatic approach 
to real complex chemical problems. A scenario for this approach to a given 
problem Is as follows: 

(1) Write a mechanism based on best available chemical Information. 

(2) Estimate rato constants and uncertainties for those steps that 
have not been previously studied. 

(3) Calculate the desired time-dependent concentrations, their 
uncertainties, and the sensitivity of these to the Input parameters. 
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(4) Measure those quantities where the combination of sensitivity 
and uncertainty are greater than acceptable. 

(7) Continually perform experiments designed to test the limits 
of accepted theory. 

It should be noted that rate constants for elementary steps may be functions 
of temperature, pressure, and the nature of other molecules present. In addition, 
some complex processes (such as smog) involve photochemical steps as well. 

For this symposium, I *111 discuss the status of Items (2), (4), and (5) 
for gas-phase homogeneous processes. 



2. THERMOCHEMICAL KINETICS 


Background 


The development of an empirically based framework for 
extrapolation and estimation is intimately tied to the avail- 
able methods for measuring rate constants. These are discussed 
separately for organizational convenience. This framework has 


! 
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* 
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been given [l] the name "Thermochemical Kinetics" by S. W. 
Benson. I have discussed "Thermochemical Kinetics" in the 
following manner with reference to tropospheric measure- 

ments , "combustion-problems [2], and modelling [3] recently, but I vrtl-1 
repeat it here for completeness. 

Therraochemi s try 

It is impossible to begin a discussion of the theoretical 
basis for critical evaluation and extrapolation of thermal 
rate data without first discussing methods for estimating 
thermochemical quantities, such as AHj t , AS° , and c£ T for 
molecules. 

Group Additivity 

When a sufficient data base exists, we have found [41 the 
method of group additivity to best fit the need for accuracy 
and ease of operation. The basic concept and assumptions 
involved in the group additivity method are as follows: 

For the disproportionation reaction 
rnn'r + snn's S rnn's + '•nn'r 

any additivity approximation assumes that A5 * A$ a * where i 
is any molecular property, and A$ a is the contribution to that 
property due to symmetry changes and optical isomerism. For 
the molecular properties of interest here, ABy -* o, Ac p>T "* °» 
and As*. - s . R In K , where K - <j<RRNn' R) cj(SNN # S) / CT (RNN 'S) 
o (SNN R) , c (X) being the symmetry number including both inter- 
nal and external symmetry. An additional term for entropy of 
mixing, due to the existence of optical isomers, must also be 
included. 

If the molecular framework NN' is two atoms or greater, 
these relationships imply the additivity of group properties, 
which include all nearest-neighbor interactions, since a group 
is defined as an atom together with its ligands (e.g., in the 
group C-(H) 3 (C), the central C atom is bonded to three H atoms 
and one C atom) . Thus the equation 

ch 3 oh + ch 3 ch 2 och 3 « ch 3 ch 2 oh + ch 3 och 3 

implies the additivity of the properties of the groups 
C-CH) 3 (C), C-(H) 3 (0), 0-(C)(H), C-(H) 2 (C) (0) , and 0-(C) 2 , if 
the appropriate A0 . 
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We have developed group additivity methods that permit the 
estimation, for many organic chemicals in the gas phase, of 
heats of formation to ± 1 kcal/mole, and of entropies and heat 
capacities to ± 1 cal/ (mole-K) , from which free energies of 
formation can be derived to better than ± 2 kcal/mole. 

It should be noted that entropy And heat capacity are 
molecular properties that can be accurately estimated under 
much less stringent conditions* than energy (or enthalpy) . 

Thus the method of bond additivity seems tc work quite well 
(± 1 cal/ (mole-K) ) for estimating the former properties, but 
not at all well (± 4 kcal/mole) for the latter. 

Structural Considerations and Model Compounds 

If sufficient thermo chemical data is lacking for the 
estimation of group properties, entropy and heat capacity can 
often be adequately estimated from structural parameters of 
the molecule. (Enthalpy estimates are more difficult, requir- 
ing. a better knowledge of potential functions than are usually 

available.) The methods of statistical thermodynamics may be 

0 0 

used to calculate Cp and S directly for chose molecules where 
a complete vibrational assignment can be made or estimated. 

Also, •" reasonable" structural and vibrational frequency 
"corrections” to the corresponding established thermodynamic 
properties of "reference” compounds may be mads. A suits? le 
choice of reference compound, i.e., one similar in mass size 
and structure to the unknown, assures that the external rota- 
tional and translational entropies and heat capacities of the 
reference and unknown compounds will be the same and that many 

of the vibrational frequencies will be similar. The basic 

o o 

assumption is that S and Cp difference can be closely esti- 
mated by considering only low-frequency motions thought to be 
significantly changed in the unknown. Fortunately, entropies 
and heat capacities are not excessively sensitive to the exact 
choice of these vibrational frequencies, and estimates of 
moderate accuracy may be made with relative ease. 

Kinetics 

The extension of thermochemical estimation techniques to 
thi evaluation of kinetic data rests largely on the validity 
of transition state theory. 
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The transition state theory expression for a thermal rate 
constant is, 

k - (kT/h) expf — flG^/RT] 

(the units are sec“l for the first order and atm"l for the 
second order) and, 

, tty* = AH^qq - TAS§Jo + < AC p * > [ (T-300) - Tln(T/300) ] 

(In the ideal gas approximation we can drop the standard state 
notation on AH* and AC*.) if the empirical temperature depen- 
dence is represented by 

k » AT B exp(^C/T) 

A - [k/h(300) <Ac p /R> ] exp[( dS3 00 “ <AC*>)/r] 

B - (<AC^> + R)/R 

C - (ah 3 * 00 " < AC p > (300)) /R 


k ■ Boltzmann's constant 




Planck's constant 

entropy of activation at 
300 K, standard state of 
1 atm. 

enthalpy of activation 
at 300 K 

average value of the heat 
capacity at constant pressure 
of activation over the 
temperature range 300-T °K. 


If we wish to express second-order rate constants in con- 
centration units instead of pressure units, we must multiply 
by RT in the appropriate units. This has the effect of 
writing: 

o' 

k » a't exp (— C/T) 

where A ' » AR and b' » B + 1 - (< ACp > + 2R) /R 

Thus, simple "Arrhenius behavior" which will be sufficient 
for lower tropospheric temperatures is characterized for 
first-order reactions by Ac* ■ — R; (AC* » AC^ ■ AC*), and for 
second-order reactions using concentration units by Ac p - 
- 2R (or ACy* » - R) . 


In the case of simple Arrhenius behavior: 
k - A exp( — B/T) 
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log A - log tek < T>/h> + AS*/R; B - <AH* + R<T>)/R 

Thus, the quanties A H * , AS*, and ACp are of interest. He 
apply similar methods to those already discussed with respect 
to thermochemistry to view rate data in a rational framework. 
These techniques are discussed in some detail by Benson[l], but 
certain points are worthy of re-emphasis here. 

We begin by classifying reactions as unimolecular or bi- 
molecular. (The only termolecular processes of interest to us 
will be energy- transfer controlled bimolecular processes. 1 

Unimolecular Processes 

Simple Fission: AB " A + B 

Complex Fission: Molecule *• Molecule + Molecule 

Cor radicall 

Isomerization: Intramolecular atom rearrangement 

Fimo lecular Proces ses 

Li;:ect Metathesis: A + BX - AX + B 

Addicion: A + Molecule Stable Adduct Greverse 

of complex fission) 

Association: A + B - A — B (reverse of simple 

f issionl 

The first thing to notice is that of all these reactions, only 
direct metathesis reactions are not subject to becoming energy 
transfer limited at high temperatures and low pressures Ci.e., 
in the "fall-off" region!). This means that not only does the 
so-called high pressure rate constant need to be estimated or 
known, but the extent of fall-off, as well. Methods are 
available for making fall-off corrections [s]. 

In hydrocarbon reactions in the troposphere, we may expect 
that most direct metathesis reactions will involve the exchange 
of a hydrogen atom between larger groups. A simple, semi- 
empirical prescription exists for estimating the value of AS* 
for these types of reactions. First, one realizes that these 
values are limited between the "loosest” possible model 
(A-factor equals gas kinetic collision frequency! and the 

"tightest" possible model in which R»*.H*»*R is represented 

/ 

by the molecule R-R . Experience using data in the range 


300<T/K < 700 has taught us that generally the AS* value 
corresponds to a transition state only slightly looser than 
the tightest possible value. 

Since the other two classes of bimolecular porcesses are 
the reverse of unimolecular reactions, we may consider them in 
that direction. (The equilibrium constant is either known or 
estimable.) Once again, using experimental results as our 
guide, we note that model transition states which correspond 
to the values of AS* are generally "tjht". That is, we may 
visualize them as minor modifications of the reactant molecule, 
usually involving some increase in rotational entropy due to 
slight enlargement of certain bonds. The dominant entropic 
feature is usually the stiffening of internal rotations as a 
result of multiple bond formation or ring formation [l]. 

Bond scission reactions present a particular problem, since 
it is particularly difficult to locate a transition state. 
Recent work [e] both experimental and theoretical, indicates 
that these reactions can be modeled with a transition state 
which becomes tighter as the temperature rises. Lover limits 
of A-factors can be estimated fairly accurately for H-atom 
metathesis reactions by making use of model transition states 
Cl]. We illustrate the method for the atmospherically important 
reaction of OH with HNO3 which is supposed to proceed via: 

HO + H0N0 2 “• H 2 0 + NO3. Reported valuesC?] are k~10 8 M" 1 s” 1 
independent of temperature. Our analysis suggests that the 
reported rate constant does not represent the reaction as 
written above. 













3 . MEASUREMENT TECHNIQUES 


Bimolecular 

Bimolecular reactions of interest are usually radical- 
molecule or radical-radical interactions , although occasionally 
a molecule-molecule reaction will be an elementary process. 
(Gas-phase ion-molecule reactions are of interest sometimes , 
and they can be included in the same general framework with the 
realization that attractive forces are longer ranged.) Since 
bimolecular rate constants can be measured only if concentra- 
tions are known, the standard "trick” is to arrange for "psuedo 
first-order" conditions by running the reaction in a large 
excess of one of the reactants. This allows the relative 
measurement of the other reactant to determine a first-order 
rate constant, k 1 , which, when divided by the concentration of 
the overwhelming component, will yield k 11 , the second-order 
rate constant of interest. Most experimental techniques take 
advantage of the psuedo-f irst-order trick over the ranges of T 
and P accessible. 

Several methods exist for concentration (or relative con- 
centration) measurements. Most prominent are mass spectroscopy 
and all the wavelength ranges of optical spectroscopy, both 
absorption and fluorescence. 

The experimental environment can also vary widely* Flow 
tubes that monitor time as distance after an interaction 
region, real-time measurements that are rapid enough to 
determine actual concentration changes, shock tubes that create 
high- temperature conditions rapidly, and some recently exploited 
infrared laser-heating techniques (vida infra) . 

Unimolecular 

Unimolecular reactions are measured using many of the same 
general techniques as for bimolecular reactions. In general, a 
substantial energy barrier must be surmounted in order for these 
processes to proceed, 30 that moderate- to-hign temperatures are 
required. In addition, since unimolecular processes are depen- 
dent on the total pressure and the specific nature of other 
constituents of the gaseous mixture, all these properties need 
to be under experimental control. 
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Vary Low-Pressure Pyrolysis (VLPP) 

A particular technique that can be used is the Very Low- 
Pressure Pyrolysis (VLPP) technique developed at SRI Inter- 
national C*3 . VLPP is ideally suited for the measurement of the 
rate of initial bond-breaking reactions in the pyrolysis of 
organic molecules where secondary reactions often interfere 
with the characterization of the initial step. The technique 
has been described in detail previously [8,9]. 

The procedure consists of allowing the reactant to flow 
through the inlet system at a controlled rate and switching 
valves in the inlet lines so that the reactant alternately flows 
through a Knudsen cell reactor or through the bypass directly 
into the mass spectrometer. Since the flow rate is held con- 
stant, the difference between the mass spectrometer signals 
for the reactant in the bypass position and the reactor position 
corresponds to the amount of reactant that is decomposed as it 
flows through the reactor. 

The data are interpreted with the aid of various steady- 
state expressions derived as shown below. At low flow rates, 
the treatment for a simple irreversible unimolecular decompos- 
ition is appropriate, and the extraction of rate parameters is 
straightforward. At higher flow rates (higher reaction pres- 
sures) and smaller escape apertures, rapid bimolecular reactions 
can compete with unimolecular decomposition and with escape from 
the reactor. These interactions must be included in the anal- 
ysis. The reason for using VLPP at pressure high enough for 
secondary reactions to occur is shown by the description given 
below of our recent study for the pyrolys*«» -of l-ethyl naphthalene sad i-ethyl 
awthrscon* [ io ] .. Briefly, observation of competition between 

unimolecular bond scission and radical recombination amounts to 
measurement of an equilibrium constant: 



Reliable measurement of an equilibrium constant can provide very 
good third-law values for &H°, since these values are not s .b- 
ject to the systematic errors that can markedly affect the slope 
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of Arrhenius plots. Thus, by providing for measurement of an 
•quilibrium constant that is otherwise not raadily measured at 
higher pressures or in static systems, the VT,PP system allows 
two largely independent measurements of bond strength and, 
therefore, a valuable internal consistency check. 


For irreversible, unimolecular decomposition of a substance 
A, only three rate processes are considered: 



A 

products 

mass spectrometer 


( 2 ) 


where R A is the rate at which the reactant is allowed to flow 
into the reactor and k^ is the first-order rate constant 
describing escape from the reactor. 


Steady-state analysis provides the expression: 


A escaped 

^reacted 


(A) 


* ' ss 

(A) - (A) 

o , ss 


ss 



( 3 ) 


where (A) _ is the steady-state concentration in the reactor 

o , ss 

when there is no decomposition. 

When recombination of the radical fragments produced by 
unimolecular decomposition competes with escape from tha reac- 
tor, the following reactions must be considered: 



A 

B* ♦ C’ 


} escape to mass spectrometer 


( 4 ) 


The controlled flow rate into the reactor is given by R A , 
and the first-order escape rate constants for the various frag- 
ments are related by square roots of their masses so that the 
escape constant for the 1-naphthylmethyl radical (B*) is related 
to that for 1-ethylnaphthalene by 

k .B ‘ * ".A' 1 ' 05 ' 


( 5 ) 
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Steady-state 
lane provide* an 


analysis of this sequence for 1-ethyl naphtha- 
expression similar to aquation (3) t 


^escaped m 

^reacted 


(A) 


ss 


1X1 .. 



k r F CA) 

* k d (3V'33J 


( 6 ) 


Equation (6) indicates that decomposition and competitive re- 
combination will provide a straight line of slope k ^/k^ and 
intercept k^/k^. In the limiting case where recombination is 
unimportant, equation (6) reduces to equation (3) . 


Fhotolytic Radical Production 

We have recently [ill combined the phenomenon of infrared 
multiphoton dissociation of organic molecules with the VLPP 
technique to produce a method for determining the rate constant 
for radical-mclecule reactions at temperatures determined by 
reactor wall temperatures and completely independent of any need 
to heat the radical precursor. (Current powerful dye lasers 
wj.ll allow this technique to be useful for UV-vis photochemical 
radical production as well.) The current application has been 
ti. reactions of CF3 radicals from CF3I, but extension to 
aliphatic and aromatic systems requires only time and funds. 

In this experiment, as previously, the effusive molecular 
beam was mechanically chopped in the second (differentially 
pumped) chamber before it reached the ionizer of the quadrupole 
mass filter (Finnigan 400) . The signal was demodulated by a 
lock-in amplifier (PAR 128 A) whose output was now stored in a 
signal averager (PAR 4202), which also served to trigger the 
laser. The two-aperture Knudsen cell was fitted with Kd win- 
dows and had an optical pathlength of 20.5 cm and a volume of 
approximately 105 cm 3 . The cell was coated with Teflon by 
rinsing with a finely dispersed Teflon slurry in a water/aro- 
matic solvent mixture (Fenton Fluorocarbon, Inc.) and curing 
at 360 C. 

The Lumonics TEA-laser (Model K-103) was operated at the 
R(16) line of the 9.6* transition at a pulse repetition fre- 
quency of .25 Hz, slow enough to permit > 99% of the reaction 
products to escape from the cell before the next laser shot. 

The multimode output of the laser consisted of a pulse of 
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approximately 5.0 J, directed through the photolysis cell after 
being weakly focused by a concave mirror (fl * 10 m) , .'hich 
gave a beam cross section of 2.67 cm2 a t the KCl-entrance 
window of the cell. 

The typical experiment consisted of averaging the time- 
dependent mass spectroscopic signal intensity of the products 
for a number of laser shots (10-100) as a function of the flow 
rate of the reactant gases at constant CF3I flow rate and con- 
stant energy per pulse. Although experiments could be performed 
on a single-shot basis, the signal/noise ratio was improved by 
averaging the results of a number of laser shots. 

The total yield of product formed in the reaction of 
interest was then determined by integration of the accumulated 
time-dependent signal of the signal averager on a strip-chart 
recorder fitted with an electronic integrator (Linear Instru- 
ments, Inc.). 

Using the apparatus described above, reaction products could 
be observed as well as the transient depletion of CF3I. An 
advantage of the low-pressure technique is that the effects of 
secondary reactions are minimized, although they must be 
considered in the data analysis. 

The following chemical reaction mechanism is appropriate for 
studying the reaction CF3 + Br2 "* CF 3 Br + Br: 

CF 3 + I 

CF^Br ♦ Br 

(wall loss not yielding CF^Br) 

No heterogeneous first-order reaction of CF 3 to produce CF^Br is 
included, since the data interpretation does not suggest it. As 
in the usual data treatment for VLPP studies, each molecular and 
radical specie escapes from the reactor with a characteristic 
first-order escape rate constant; for CF^I, Br 2 , CF^, and CF^Br, 
the escape rate constants are kj, k^, k^, and k^. 


CF 3 I + nhv “* 
CF, + Br, i 

cf 3 -• 
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Analysis of the reaction mechanism and solution of the 
appropriate differential equations give an expression for the 
time-dependent mass spectrometer signal due to CF^Br. The 
total yield of CF^Br, (Y) , is related to the rate constants as 
follows: 



where pseudo-first-order conditions are assumed to hold. In 
this expression, a is a mass spectrometric sensitivity factor, 

0 is the fraction of the initial CF^I that is dissociated by 
the laser pulse, and V is the volume of the cell. The initial 
L CF 3 l3o ■ f cf 3 I^ V * *2^ ’ where FCF3I is the flow rate of CFjI 
into the reactor, and V is the reactor volume; similarly, [Br 2 ] 

■ F Br ^/(V * k 3^* For eac ^ F Br2' tWO esca P® apertures can be 
used, giving two different values for the escape rate constants, 
corresponding to the "big" and "small" apertures. Plots of Y“1 
versus F”^ give two straight lines with intercepts c g and c b 
(small and 2 big apertures) given by i = s or b) : 


‘2i 


O0VF, 


cf 3 i 


The slopes of the straight lines are given by 


m. 


c l V(k 4i * V k 31 


Laser-Powered Homogeneous Pyrolysis (LPHP) 

In VLPP experiments aimed at obtaining unimolecular rate 
information, we rely on the walls of the VLPP reactor to be the 
source of heat through gas-wall collisions, while at the same 
time being non-catalytic for the destruction of the substrate. 
Sometimes this latter condition is not met. Therefore, we may, 
in addition to VLPP, use the technique of laser-powered homo- 
geneous pyrolysis [12] . This technique is a valuable complement 
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to the VLPP procedure, since it essentially provides a "wall- 
less" reactor. The total pressure in the reactor is on the 
order of 100 torr, consisting mostly of bath gas and SF g . An 
IR laser is used to heat the strongly absorbing sensitizer (SF g ) , 
using a wavelength at which the substrate does not absorb. The 
SF g transfers its thermal energy by collision to the substrate 
molecules, and decomposition takes place. As described by Shaub 
and Bauer [ 12 ] , the technique worked well for compounds of 
relatively high vapor pressure, using a static reactor system 
and a CW laser; we are currently adapting the technique for 
application to poly-nitro aromatics and other low-vapor-pressure 
substrates, using a pulsed C0 2 laser and a flow system with GC 
detection. The reasons for these modifications are: (1) When 
operation is with a pulsed laser, reaction times are short 
because of rapid cooling by contact with the off-axis cell 
contents, and secondary reactions are either unimportant or can 
be minimized by suitable choice of a scavenger, and C2) when 
low-vapor-pressure substrates are being studied, quantitative 
recovery and measurement of products and unreacted starting 
material is simpler with a flow system. 

We have preliminary evidence from this technique [” 13 ] which 
indicates that nitrotoluenes decompose by NC^-aromatic bond 
scission. This seems to be true of those with ortho-methyl 
substitution as well, in contrast to earlier reports. 

4. DISCUSSION 

m here are various limitations on the usefulness of thermo- 
chemical kinetics. These range from quantitative to qualitative 
uncertainties. We need to test some of the preceding ideas with 
experiments conceived for just that purpose. Since the ideas 
are based on the transition state theory formalism, it is impor- 
tant to address the question of limits of validity of transition 
state theory. In general, these testing reactions should be 
measured under conditions where isolated reactions can be 
observed, as the extraction of individual rate constants from 
complex reacting systems is fraught with difficulty. 

There are many examples of reactions for which rate con- 
stants have been studied, but product studies are lacking. Thus, 


in reactions of OH with olefins current smog and combustion 
models must arbitrarily decide on branching ratios. This is 
equally true in aromatic systems. 

In all of the above discussion of estimation of rate data, 
the importance of thermochemical values for all species has 
been emphasized. It is particularly important to have a good 
set of values for the entropy and heat of formation of organic 
fre -2 radicals. 

Very few spectroscopic assignments exist for modest-to- 
large size organic free radicals . Entropies (and heat capa- 
cities) have generally been estimated by methods discussed 
earlier. Uncertainties arise from changes in hindered rotation 
barriers and changes in skeletal bending frequencies. 

Many bimolecular reactions are not direct, rather involv- 
ing bound complex as an intermediate. These give rise to what 
seem to be unusual parameters for bimolecular processes [14-]. 

I 

A simple treatment for the pressure dependence of uni- 
molecular processes, which would lend itself to easy use in 
large models, is seemingly close at hand, but care must be 
exercised. 

In summary, we have a framework for the codification and 
extrapolation of rate data, but much testing and modification 
will be necessary. 
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Abstract: Rate constants for i on-molecule association reactions have been calculated according to a simple model. The model, 
which mimics RRKM results, is constrained to reach limiting values corresponding to ADO (or Langevin) collision frequencies 
at the high-pressure limit and the strong-collison energy-transfer rate constant at the low-pressure limit. For those reactions 
at the low-pressure limit, the required information is the density of states (entropy) of the association complex. Low-pressure-limit 
rate constants and those in the "fall-off" regime are compared with data. 


Introduction 

Ion-molecule clustering or association reactions are of interest 
in physical, organic, and atmospheric chemistry. The overall 
process, typically written 

A* + B + M === AB* + M (1) 

can more descriptively be illustrated as 
A Aim 

A* + B = (AB*)* * AB* (21 

*» 

The collision rate constants. k< and k„ are given by Langevin' 


1 Postdoctoral research associate. 


and ADO theory^ -4 as appropriate. 

It has already been pointed out 3 6 that this type of reaction is 
totally analogous to neutral radical-combination reactions. Thus, 
since it is really the microscopic reverse of a unimolecular 
bond-scission process, all the considerations associated with the 
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Ser A. 3*9. 1 (1978); (e) E. Herbal. J Chem. Phys , 70, 2201 (1979); (d) 
D R. Bates. ■ bid.. 71. 2318 (1979) 


Kinetics for Ion- Molecule Association Reactions 
Table I. Low-Pressure-Limit Rate Constants* 
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reaction 

k/cm* 

ealed 

molecule' 1 s'* 

exptl 

ref 

(1) 

H,cr + H,0 + Ax - H,0-H,0* ♦ Ax 

5(-28) 6 

6 (-28) 

18 




2 (-28) 

19 


H,cr ♦ H,0 r He - H,OH,0* + He 

4 (-28) 

6.7 (-28) 

11 




7.2 (-28) 

20 




1.2 (-27) 

12 

(2) 

NCT + CO, + N, - NO*CO, ♦ N, 

2 (-29)' 

15 (-29) c 

21 

(3) 

O,* + N,0 + N,0 - 0,’N,0 + N,0 

3<-29) c 

5.2 (— 29) c 

22 

(4) 

0,* + 0, + 0, -• 0,* + 0, 

1.9 (-30) 

15 (-30) 

23a 




2.8 (-30) 

23b 




14 (-30) 

23c 

(5) 

0,* + N, + He-O.’N, + He 

1.5 (-29)“* 

1.9 (-29) d 

22 



7 (-31) 

8 (-31) 

15 

(6) 

o,* + h,o + o, - o,*h;o + o. 

1 (-28) 

1 (-28) 

24 




13 (-28) 

25 


* Calculations and experiments at 300 K unless otherwise cited. b Number in parentheses is the power of ten. e 200 K. d 80 K. 


origin and nature of pressure dependence, as well as temperature 
dependence, arise when quantitatively discussing the rates of such 
processes. Ion-muecule association reactions are therefore treated 
by RRKM theory 3 ** in order to Minify and extrapolate the data 
»o various conditions of temperature and pressure. 

RRKM treatments involve complicated computer codes and 
'wailed molecular information, and Troe and co-workers’-* have 
presented a model that allows a simpler calculation We have 
adapted Troe's model to ion-molecule reactions, by taking into 
account *iie differences in the collision process. 

Back{r<<ad 

The ,.todel for the calculation of the association rate constant 
can be treated in three parts: high-pressure limit (second order), 
low-pressure limit (third order), and intermediate- pressure regime 
(fall-off region). 

FOth-Premre Lint. For the high-pressure-limit rate constant, 
we use Langevin and ADO theory to calculate k^ the collision 
frequency between A* and B. 

Low-Pressure Limit We First calculate (£,)<>, the low-pres- 
sure-limit rate constant in the unimolecular direction, for the 
decomposition of AB* 'see eq 1). Wc then calculate the low- 
pressure-limit rate constant of interest from the overall equilibrium 
constant ( K ■ k t /k,). 

In complete analogy with Troe 7 we calculate the collision 
frequency between AB* and M by using ADO theory and a simple 
harmonic oscillator model followed by a few correction factors 
which account for (1) anharmonicity, (2) energy dependence of 
the density of stales, (3) overall rotation, and (4) internal rotation, 
i.e. 

)kT 

K * 0Zado 0— exp(-E 0 /kT)F lhk F E F rai F m „ (3) 

where d * collision efficiency (vida infra), Z\do m ADO collision 
frequency, = harmonic oscillator density of states. * 
vibrational partition function, and £ 0 = critical energy (bond 
dissociation energy at 0/K). (Explicit formulas 7 * for the various 
factors used in this calculation are listed in Appendix I.) Since 
one of the oscillators is represented by an /T* potential, there may 
be a difference in anharmonicity correction factors between 
ion-molecule and neutral-molecule reactions, but the difference 
due to this one oscillator should be small, especially when the 
associated molecular ion has many atoms, d is assumed to be 
unity. 3 Equation 3 arises from 7 

(*,)o _ [AB»*1 ^ g A> -» _ X/ (£V " />T d£ _ (*,) o 

*«[M] [AB*’] (?as* Gas* &Z ado 

(4) 


(7) J. Troe. J Chtm Phv i . 6*. 4758 1 1977) 

(8) K. Luther and J. Troe. "Proceedings of ihe 17th International Sym- 

posium on Combusuon". Leeds. England. 1978. p 535. 


Table tL Collision Efriciencea-Small Complex 

O,* + H,0 ♦ M — ;0,H,0)* + M 


M 

*axpti a 

ref 

He 

8.7 X KT” 

24 


9 x 10-” 

25 

Ax 

1.7 X 10-" 

24 


10 X io-“ 

25 

0, 

13 X 1(T“ 

24 


1 X ltr*' 

25 

N, 

15 x lO”* 

24 


18 x 10-” 

25 

a cm* molecule' 1 s'* 

; 300 K. k ctled = 1 X 

10"“ cm* molecule'* 


f 1 for all cases. 


Sole that this is independent of transition-state or RRKM theory. 

This low-pressure-limit decompsition rate constant is then 
multiplied by the equilibrium constant to calculate the low- 
pressure-limit association rate constant. Unfortunately, ther- 
mochemical data are not always available. In these cases, we 
estimate some of the thermochemical data from geometry and 
vibrational frequencies to calculate the equilibrium constant. 

FaJl-OfT Rrgiou. After we have obtained both the high- pres- 
sure-limit rate constant and the low-pressure-limit rate constant, 
the intermediate-pressure rate constant is obtained from a for- 
mula, 1 empirically designed to mimic RRKM curves, which has 
been successfully used’ 10 in neutral molecule associauon reactions. 
For each temperature, this formula is 

MM] rli*|loftMMl/*.)lir' 

1 + kolM]/*/ c 

where k 0 and k m represent low-pressure-limit and high-pres- 
sure-limit rate constants, respectively. 

F q is a constant that can be evaluated* for different conditions. 
We fir.d that for our limited purposes, a value of F c =* 0.6 is 
sufficient to reproduce values computed with a full RRKM 
calculation to within a factor of less than 2. ( More accurate 
reproduction of RRKM can be obtained by usuig the more detailed 
formulas in ref 8.) 

Results and Discussion 

(A) Low- Pressure- Limit Association Rate Constants. The re- 
sults of low-pressure-limit rate constant calculations for six as- 
sociation reactions are listed in Table I. The molecular parameters 
and tbermochemical data used in the calculations are listed in 
Table IV As shown in Table I, the agreement between the 
experimental data and the calculated rate constants is excellent. 
(This is also true for reactions shown in Tables 11 and III.) This 
not only means that this model works extremely well in relating 
the thermochemical and the kinetic data but also allows us to 


(91 R. Zellner. Btr Bunstnyts. Phyt. CStm.. 12. 1172 (1978) 

(10) D M. Golden. J S. Cheng, end A. C Baldwin. NASA Workshop 
on Stratospheric Chemiatry JPL Publication No. 79-27. Evaluation No. 2. 
April 15. 1979 
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Tabic I1L Collision Efficiencies- Luge Complex 


C.H,* ♦ C.H, ♦ M - (C.H.),- » M 


M 

10 ,, A: c *icd° 


ref 

He 

3.8 

4 X 10-” 

26 

Ne 

1.2 

<5 x 10-** 

26 

Ar 

1.3 

7 X 10-” 

26 

Kr 

1 

1.1 x 10-” 

28 

Xe 

1 

1.6 x lC** 

26 

C.H. 

2 

1.2 X lO”* 

26 


• cm* molecule'* f 300 K. 


understand and predict many gas- phase ion-molecule association 
rate constants. 

Of particular note is the inference that f or some reactions that 
are difficult to study, we may predict rau constants if we have 


the requisite thermochemical data. In fact, we argue 6 * that pursuit 
of the thermochemical information would be the more useful 
endeavor. 

(B) Fall- Off Phenomena. Ion-molecule association rate con- 
stants will change from third order in the low-pressure regime 
to second order at the high-pressure limit. At the medium 
pressures where the transition takes place, "fall-off phenomena 
are expected to be observed. RRK.M treatments of this are 
adequate, as previously demonstrated. In this paper, we use 
eq 5 to predict the fall off. The result of this prediction is very 
close to an RRKM calculation which employs the same values 
of fco and k m . 

A few experimental techniques have been widely used in 
studying the kinetics of ion-molecule association reactions. Most 
of them employ relatively low-pressure conditions such as the 


Tabts IV. VibraUonal Frequencies and Thermochemical Data 

complex* freq A H. b AS* ref complex* freq A H, b AS* ref 


(1-1) H,OH,CT 
410(1) 

3100(4)* 
1800(1) 
1600(2) 
500 (4) 
410(1) 
142(2) 
50(1) 

AW = -31.6, 
AS = -24.6 

6s 

(1-2) NO- CO, 

1343 

667 

2349 

2377 

1600 

500 

100 

100 

500 

A/fw — 11.3, 
AS = -25.0 

21 

(1-3) 0,*-N,0 

1580 

AW = - 12.9, 

22 


1275 

589(2) 

2225 

1275 

400 

AS = -29.0 

estimated 


1483(2) 

1326(1) 

1310(1) 

1178(2) 

1146(1) 

1089(10) 

1037 (2) 

1010 ( 1 ) 

992(1) 

985 (1) 

970(2) 

847(3) 

800(2) 

703(1) 

600(1) 

400(7) 

300(2) 

70(5) 

50(1) 

(F) CH,NH, CH,NH,* 3300(4) Af/ = -21.7, 

2970(6) AS = -23.6 3 

1600(1) 


100 

100 

(1-4) O.* 1110 

1043 
800 

400 

400 

100 

(1-5) 0,*N, 1100 

1100 
1200 

400 

400 

75 

(1-6.11)0,- 11,0 1100 
1043 
3600 
3600 
400 

1000 

1000 

1400 

100 


A/f = -9.6, 27,23b 

AS = -20.6 27,23b 


aH = -7, 22,27 

aS = -20 


tT)C,H,aCH,- 


aH - - 16. 
AS = -17 9 


24. 25 


1450(6) 

1175(4) 

1000(3) 

900(2) 

820(2) 

740(2) 

600(1) 

400(1) 

300(2) 

200 ( 2 ) 

175(2) 

100 ( 1 ) 

3002 AH = -30.7, 

580 AS = -30.7 

3184 12) 

1383(2) 

1020 ( 2 ) 

1460(1) 

1020 ( 2 ) 

1000 ( 1 ) 

1350(1) 

3000(5) 

1150(2) 

400(1) 


dll) (C.H.),* 


3082(6) a// =-15.1, . 

3080(2) AS » -23.2 61 

3062(1) 

3060(1) 

3047(2) 

1597 (6) 

1596(2) 


125(2) 

100 

100 

550(5) 


“ 1-1 through 1-6 refer to Table I. II tefers to Table II; 111 refers to Table 111; F refers to the figure; T refers to the text b AM for the reaction. 
A’ ♦ B — AB" in kcal inoP (300 1C). e AS for the reaction: A* * B - AB* m cal mol'* K* 1 (300 K). 
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Flgar* I. Pressure dependence of k, at 300 K for the reaction CH,NHj 
4- CH,NH,* + M ^ CH,NH r H 4 NH J CH, 4- M (O. refl 3; 0, ref 14): 
( — ) a poon beat estimates from methods outlined heron: (---) reduction 
of t, by s factor of 5 to accommodate data of ref 14; (•••) comp omiie 
value to beat accommodate high- and low-pressure data. 

flowing-afterglow technique 11,12 and ion-cyclotron resonance 
spectroscopy at mili torr pressures, leading in general to near 
low-pressure-limit rate constants. However, some of the reactions 
may still show deviation from the low-pressure limit, particularly 
when mass spectrometric techniques in the torr range are em- 
ployed. Extrapolation out of the experimental pressure range will 
be unreasonable 12 unless this is taken into account. 

Prediction ts compared to experiments using, as an example, 
the reaction 

CH)NHj* 4- CHjNHj 4- M — CHjNH, 4 CH,NHj + M 

Two groups have studied this reaction by two different techniques. 
Netlson et al. 12 studied this reaction over the pressure range of 
1 x i (T* to 3 x 1(T 2 torr of CH,NH 2 by using the ioo-cycloiron 
resonance mass spectroscopic technique and concluded that the 
low-presaure-limit and high-pressure-limit rate constants were 1.8 
X 10" 23 cm* molecule' 2 and 2.1 x HT 11 cm 2 molecule ! s' 1 , re- 
spectively. 

Meot-Ner and Field 14 used a pulsed high-pressure mass spec 
trometer to study this reaction over the pressure range of 0.5-2 
torr isobutane and concluded that the low-pressure-limit and 
high-pressure-limit rate constants were 15 X 10“ 24 cm* molecule 2 
s" 1 and 4.0 x 10* 1# cm 2 molecule' 1 s' 1 . 

The disagreement between these two reports results, in part, 
from systematic differences in the experiments and, in part, from 
an incorrect extrapolation of the data. 12 Figure 1 is illustrative. 
The higher pressure data 14 can be fit by assuming that It. ■ Ic^jo 
■ 1.75 x 1 0~* ctr 3 molecule' 1 s' 1 and k^ ■ 3 x 1 O' 2 * cm* mole- 
cule' 2 s' 1 . We compute a value of *o * 1.1 x 10" 25 cm* molecule' 2 
s' 1 (in isobutane) by using the methods outlined herein. This is 
not out of the bounds of possible error. However, the lower 
pressure data 12 fit well the computed vaiue (in methylamine) of 
Jlto ■ 1.5 x 10 22 cm* molecule 2 s' 1 . Taking this latter at face 
value suggests that the higher pressure data are in error by a factor 
of ca. 2. Of course, a value of lc 0 ■ 6 x 10" 2 * cm* molecule' 2 s' 1 
splits the difference accommodating reasonable opposing errors 
in the two data sets, as well as an acceptable deviation from the 
calculation. (This latter line gives values within a factor of 2 of 
the RRKM calculations of Olmstead et al., 2 who had only the 
high-pressure results at their disposal.) 


( 1 1 ) V M Bierbeum, M F UoKle. and F Kaufman. J Chem Phyi . *5. 

2715 « 1976). 

(12) C. J. Howard. V M. Bicrbaum. H W Rundle. and F Kaufman. J. 
Chem Phyi . 57. .’491 (1972) 

(13) P V Neilson. M T Bowen. M Chau. W R. Davidson, and D H. 
Attt, J Am Chem Sot . IM. 1649 (1975) 

(14) M Meo« -Ner and F H Field. J Am Chem Soc .. 97. 5339 ( 1975). 

(15) L. Baaa. W J. Cheanavich. and M T Bowers. J Am Chem Soc.. 
79. 1501 (1979) 


Another example of the usefulness of this approach is illustrated 
for the reaction 

CjHj* + CHjCl + CH« - (C 2 H,ClCH,r + CH 4 

The rate constant for this process has been measured 1 * at ~4 torr. 
The authors, having supposed the process to be in the low-pressure 
limit, interpreted the data to yield k m 1.7 x 10~* cm* molecule' 2 
s' 1 at 300 K. However, this rate constant multiplied by the 
pressure (i.e., 1.7 x 10 -2 * cm* molecule' 2 1 ' 1 x 1.3 X 10 17 molecule 
cm 2 ) yields kf ■ 12 X 1(T* an 2 molecule' 1 s' 1 . The limiting values 
(/co(cm‘ molecule' 2 s' 1 ) ■ 3.5 x 10* 22 , calculated by methods of 
ref 7; It. (cm 2 molecule' 1 s' 1 ) ■ 102 X KT*. calculated from ADO 
theory along with eq 4) predict a value of k(cm 2 molecule' 1 s' 1 ) 
■ 1.8 X KT*. well within all the error limits. Notice that rather 
than following the author's supposition, at this pressure the process 
is very close to the high-pressure limit This is to be expected for 
this large complex which is strongly bound (30.7 kcal mol' 1 ). 

(C) M Effect (N stare of the Bath Cm). The assumption that 
3 is unity (the strong collision assumption) has been discussed 
previously. There is no reason to expect that the efficiency for 
the rare gases in not somewhat less than unity. We report here 
that for small ion-molecule complexes (AB 4 )*, 8 “ 1 (this as- 
sumption is good within experimental scatter), but that for large 
(AB 4 )* (see Table III), the efficiency of small bath gases seems 
to be low. These effects can be demonstrated with the following 
two reactions. 

For small AB 4 molecule:, the reaction 

O, 4 + HjO + M — Oj 4 -H 2 0 + M 

is used as an example. The M effect in the calculations resulted 
from the difference in Langevin collision frequencies which is 
proportional to (aM,MM,/ a M/*M,) l/J for different M, where a is 
polarizability and a is the reduced mass of AB 4 and M. Usually, 
the calculated values of the rate constant for different M's are 
fairly dose. The calculated and the experimental results are shown 
in Table II. All the experimental rate constants for different M's 
are within a factor of 3 compared with the calculated valuta, which 
lends credence to the notion that all bath gases are ‘’good 4 third 
bodies. 

For large AB 4 molecules, the reaction 

C*H « 4 + C»H« + M — (C»H 4 )j 4 + M 

is used as an example. Table III shows all the experimental and 
calculated results. These results suggest that all rare gases are 
inefficient energy receivers or earners, while C 4 H* itself works 
well. This suggestion of inefficient energy transfer from a large 
ion-molecule complex to “small" bath-gas molecule warra .ts 
further investigation. 

(D) Otter Coasidmtiowa. The temperature dependence of the 
rate constant can also be calculated according to this model. The 
results arc similar to the results from RRKM calculations, which 
have also been shown to yield good agreement with some ex- 
penmen ts. 5 -* 

The uncertainties of around a factor of 2 or 3 in the low-pressure 
rate constant calculation are mainly due to current insufficient 
knowledge of the tbermochemical data and of the geometry 
(structure) of the associated ion-molecule complex. For these 
association reactions, the influence on the rate constant from 
uncertainties in enthalpies is not as significant as those from 
entropies. This has been discussed previously with respect to both 
neutral and charged reactants.**- 1 ’ Therefore, it seems clear that 
any studies designed to yield entropic information on the ion- 
molecule complexes created by association reactions such as those 
dealt with herein will enhance the ability to predict rate constants 
for these processes. 
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CoaciwMM 

Thermal ion- molecule association reactions may be treated via 
RRKM theory 3 * 15 and thus also by the simplified computational 
schemes of Tree and co-workers •' if neutrrl-neutral collision 
frequencies arc rrp laced by ion- neutral frequencies (Langevut and 
ADO). Crucial inputs are the thermodynamic parameters of the 
association complex, especially the entropy. 
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Appendix L List of Formulas 

(a) ADO collision frequencies 

Z ^° m2 ^[ a ' li + C ^{^f)'] (A1) 

where q is the charge of the ion. m is reduced mass between 
colliders, a is polarizability of the molecule, is the dipole 
moment of the molecule, C is the correction factor as derived in 
ref 2, and k is the Boltzmann constant. 

(b) Harmonic oscillator density of states 
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where £ is the energy, £ z is '/iZia/Ari, r is number of oscillators, 
a(£) is the empincal Whitten-Rabinovitch corrector, and e, are 
the vibrational frequencies. 

(c) Anharmonicity correction factor 
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where m is the number of cscillaton that disappear during the 
dissociation reaction. 

(d) Energy dependence of the density of states corrector is given 
as 
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Abstract: The significance of entroplc barriers to dissociation 

and recombination in the prototype ionic system CH^ + CH^ + has 
been investigated. Ion-molecule systems are shown to react through 
an entirely different dynamics than neutral systems, due to intrinsic 
differences in the shapes of the relevant potential surfaces. 
Consequence* with regard to the interpretation of experimental rate 
parameters in tne ion-molecule area are discussed. 
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Introduction 


To a first approximation, ion-molecule dissociation aid 
recombination reactions might appear to be analogous to the cor- 
responding neutral processes of nni molecular dissociation and the 
reverse association. In this paper, however, we show that while 
much of the analogy is valid, ion-molecule systems react through 
an entirely different dynamics than neutral systems. Specifically, 
ion-molecule reactions proceed through "orbiting" transition 
states, even though radical-radical processes occur via "tight" 
transition states representative of an entropic barrier. 

Uninwlecular reaction rates of many neutral reactions have 
been shown to be well modeled through the use of variational transi- 
tion state theory (VTST).^ ^ Problems arise in the application of 
traditional transition state theory for reactions in which there 
is no potential energy maximum between reactants and products, 
such as simple bond homolysis; an example is the methyl radical/ 
ethane system, whose potential is outlined in Figure 1. Depicted 
are curves for J ■ 0 and some average J ■ (j); the latter shows the 
centrifugal barrier (greatly exaggerated) that results from angular 
momentum conservation requirements. Neglecting minor centrifugal 
additions onto the carbon-carbon (CC) bond stretch reaction 
coordinate, a potential maximum occurs only at infinite methyl 
radical separation. One is faced with the problem of where to 
place the transition state. In his analysis of this same system, 
Hase located the transition state at the site of minimum number of 
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bound quantum states. which of necessity lies on the side of 
a potential "hill". In doing so he was successful at modeling the 
rate parameters both of the dissociation of ethane and the recombine* 
tion of two methyl radicals. 

In terms of a statistical theory of unimolecular reactions, 

the above treatment takes on the following form: The micro- 
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canonical rate constant for dissociation is written 


k(E) 


olssh 

hN(E) 


where G^(S^) * G^(E-E ) is the number of states in the transition 

c 

state with non-fixed energy E - E £ ■ E*. E £ is the "critical 

energy" or barrier height, and N(E) is the density of states in 

the decomposing molecule at total energy E. The purpose of the 

VTST algorithm (often stated in terms of flux through phase space^) 

is to minimize k(E) in the course of minimizing the number of states 

G + (E*). For a canonical ensemble this sum of states minimization 

is equivalent to locating the point of highest free energy along 

the reaction coordinate; thus the procedure accounts for possible 

entropic as well as enthalpic barriers in the course of a reaction.®*’ 
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This free energy criterion is arguably superior to the more 
traditional potential energy criterion for transition state place-- 
ment. In any event, numerous Investigators have found this extended 
definition valuable for rate parameter modeling in the neutral 

domain.* ^ 

In contrast to neutrals, ion-molecule association processes 
have been modeled with an orbiting transition state located at the 


maximum of the centrifugal potential (see Figure 1; the actual 

•4 

potential will look slightly different owing to the long-range r 

-6 

ion-induced dipole Interaction as opposed to the r attractive 
potential in the neutral case) . A sum of states "bottleneck" on 
the side of a potential hill, however, could possibly control the 
rate parameters of such association processes. ^ Indeed, the 
significance of the minimum quantum sum transition state in ionic 
systems finds support in a number of recent reports. For instance, 
entropic barriers have been invoked by Chesnavlch e_t al . to explain 
the rate parameters of C^Hg + decomposition,* 0 and by Kebarle to 
model the temperature dependence of rates of hydride transfer 
between alkanes and alkyl carbonium ions.** 

Nevertheless, in general this model fails to explain the 
behavior of experimentally determined rates of ion-molecule reac- 
tions. We note two important facts: 

1) Many ion-molecule reactions proceed at the Langevin (or 
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encounter) limit; 

2) If a transition state of the kind described by Hase^ were 
always operative, ion-molecule reactions would never proceed at 
the Langevin limit. 

We conclude that the reaction coordinate leading directly to 

covalent bond formation is not the operative reaction coordinate 

/ 

for ion-molecule association. Ion-molecule reactions differ from 
radical-radical combinations due to the strength of the long-range 
interactions. These ion-dipole and ion-induced dipole interactions 
lead to reactive intermediates formed via the orbiting transition 


state; these relatively long-lived complexes can either decompose 
Co products or break apart to regenerate reactants. Consequently 
ion-molecule association processes with no potential energy barrier 
are encounter controlled. Those reactions that proceed more slowly 
than the encounter limit, do so as tne result of potential energy 
barriers rather than as a result of the entropic barrier introduced 
by minimizing G + (E + ). 

Calculations 

For our computations we have chosen the reaction 

ch 3 + + ch 3 > ch 3 ch 3 + 

i.e. the association reaction of a methyl radical with a methyl 

cation to yield ethane radical cation. This system is the ionic 

5 a . c 

analogue of the well-studied neutral reaction 
CH 3 + > C^CEj 

In each of the above reactions 'there is no potential maximum along 
the (incipient covalent) carbon-carbon bond distance reaction 
coordinate (see Figure 1) . In fact, our calculations show (see 
below) that the variation of numbers of states along the reaction 
coordinate is remarkably similar for both the neutral and ionic 
systems . 

In order to compute numbers of states, one needs to specify 
three quantities that vary with the reaction coordinate r: the 


ncn-fixed energy E(r), the vibrational frequencies v^Cr), where 
the index i runs over the 3N - 7 internal degrees of freedom (the 
CC stretch is not included since it is the reaction coordinate), 
and the shape of the potential curve V(r) along the reaction 
coordinate. The non- fixed energy is approximated as 

E(r) - <E> - V eff (r) 

The average energy (E) is estimated through the Boltzmann-weighted 

6 

sum of states expression 

E(E - E.)G(E - E )exp[-(E - E ) /kT) 

/r\ . 2 u 0 

W IG(E-E 0 )exp[-(E-E 0 )/VtT] 


where G(E-E^) is the sum of molecular quantum states between 

energies Eg and E, Eg is taken to be the potential energy at 

infinite fragment separation, and the sunsnations are taken over a 

suitably fine grain of energies E. The expression is evaluated at 
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the orbiting configuration (located at the Langevin fragment 
separation of 4.95 A), and (e) is assumed to be a constant indepen- 
dent of fragment separation. For T ■ 300 °K and 600 °K one finds 
(E) ■ 2.30 kcal/mole and 5.85 kcal/mole, respectively. 

An approximate effective potential function V e ^(r) is con- 
structed by summing the two components V. , ._(r) and V (r) . 

4 . —o - 1 z roc 

14 

^4-6-12 ^ t * le conmon ^ orm an ion-molecule interaction 


-6-12 ' 2 ' ‘XT/ - 3 < 1 - Y >(f)‘] 


( 1 ) 


where D g is the dissociation energy, r g is the equilibrium bond 


length, end 0 < V < 1 is a parameter that measures the relative 

_4 

importance of the r ion-induced dipole interaction term compared 

to the r”^ term typical of neutral-neutral interactions. The 
-12 

r term corresponds to a repulsive core, y is determined given 
D g , r^, and the polarizability a of the neutral fragment via 


V - 1 


2 

e a 


3 D r 
e e 


where e is the charge on a proton. 

The dissociation energy of CH^CH^"*" is estimated to be 78.3 
kcal/mole, given the known heats of formation of CH^»^ 
and CH^CH^,^ and the adiabatic ionization potential of ethane. 17 
In their theoretical study of the radical cation, 1 ^ Pople and 
co-workers have determined the equilibrium CC bond distance to be 

2.0 I, where the ionizing electron has been removed from the CC 

.3 

sigma bond. With oKCH^) ■ 2.2 A one can thus solve for y ■ 0.82 
-4 

(indicating seme r character in the potential) and generate the 
4-6-12 potential function (1). 

In addition, one must account for conservation of the orbit- 
ing angular momentum of the two fragments about their common center 
of mass. This is done by adding onto the potential ^. 5 .^ 2 ^ 
term 


V „(r) 
rot ' 


m (J(J + l))n 

2 2 
8tt nr z 


( 2 ) 


where J is the quantum number for rotation, and 11 is the reduced 
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mass. The average value (j) is estimated to be about 30 at 

room temperature. V t (r) constitutes a minor correction to 

V. c. n( r ) at reasonably low temperatures and (j) values. The 

resultant curve V e ^(r) 1* plotted in Figure 2. 

Loosening in the critical configuration along this curve was 

accomplished by lowering uie frequencies in the four rocking 

motions of the cation (analogous to the two doubly degenerate 
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frequencies and of neutral ethane ) according to a decay- 
ing exponential in r,^ 1 l.e. 

Vg * 822 exp[-0.73(r - 1.54 k) } cm * 

" H90 exp{-0.73(r - 1.54 A) ) cm 1 

The choice of the factor 0.73 will be explained below. 1.54 A is 
the equilibrium bond length of neutral ethane; thus note that this 
bending mode loosening gives lower frequencies in the radical 
cation equilibrium configuration, in accord with theoretical pre- 
dictions.^ The other degrees of freedom in the critical conflgura- 
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tlon are taken to be the remaining normal vibrations of ethane, 
except that the CC stretching mode becomes the reaction coordinate, 
and the methyl groups are allowed to rotate freely.^* Sums of 
states were calculated by a direct count of vibrational states and 
a semi-classical count of rotational states.^ 

For comparison we also calculated the sum of states at the 
orbiting transition state. In our computations we treated the 
..ysCem as ivo wholly disengaged fragments, with geometries and 
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vibrational frequencies corresponding to those reported lor the 

21 22 
methyl radical, and with free tumbling of the fragments. This 

model for the orbiting transition state has been referred to as 
23 

the Gorin model for neutral molecule decompositions. The factor 
of 0.73 used above in the exponentials describing bending mode 
loosening was chosen such that both models — the loosened vibrational 
and the fragmented— yield the same number of states at the orbiting 
CC separation of 4,95 A. 

See Table 1 for the input parameters for all of the sums of 
states computations. 

Results and Discussion 

Results of CH^/CH^" 1 " Model Calculations 

Using the above parameters for the CH^/CH^" 1 " system, we find 
the sums of states minimum to occur at a CC separation of 4.3 A, 
where the potential attraction is of the order of 1-2 kcal/'mole. 

In Figure 3, we plot the potential and the sums of states against 
the reaction coordinate. According to statistical unlmolecular 
reaction rate theory, ** the microscopic reaction rate constant k(E) 
is proportional to the sum of quantum states G^(E^) available to 
the transition state, where E and E are the non-fixed energies of 
the reactant and transition state, respectively. Thus the effect 
of the minimum is such that if the unlmolecular dissociation of 
ethane cation were to occur along this reaction coordinate (bond 
stretch), the thermal rate constant for dissociation at 300 °K 
would be about half of that predicted if the orbiting complex were 
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association 


taken to be the transition state. The rate of CH^/CH^* 
would be likewise depressed. 

Note that one can interpret the minimum sum transition state 
in terms of a "locking" of the external rotations of the fragments 
as they recombine. The locking process converts rotations into 
vibrations, thus lowering the number of states; the play between 
this locking, and an increased amount of non-fixed energy as the 
recombination reaction proceeds, determines the location of the 
transition state. Thus we shall henceforth refer to the minimum 
sum transition state as the "locked-rotor" transition state. 

In Figure 4 we plot the variation of sums of states with 
temperature for each of the models outlined above. The positions 
of both the locked-rotor and orbiting transition states along the 
reaction coordinate are relatively weakly dependent on non-fixed 
energy, and are here assumed to be Independent of temperature. 

Hase was able to parameterize loosening in the ethane homolysis 
process so as to locate the locked-rotor transition state and fit 
experimental rate parameters.^* Since analogous data for ethane 
radical cation dissociation do not exist, the placement of the 
locked-rotor transition state on the ion-induced dipole attractive 
coordinate is admittedly somewhat arbitrary. Nevertheless, we feel 
that despite the quantitative uncertainties of our model, several 
important observations can be drawn cnceming the CH^/CH^* system, 
and, we suggest, for unlmolecular ionic dissociation and association 


reactions in general: 
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First, the locked-rotor transition state exists for ionic as 
well as for neutral reactions. Second, the locked-rotor transi- 
tion state is located "Inside" the orbiting complex, that is, at 
a smaller fragment separation on the coordinate leading to CC 
covalent bond formation. This makes good physical sense: the 

free rotors can lock up only after they have passed the orbiting 
complex in the recombination process. Third, the locked-rotor 
configuration is located in a potential regime where only a few 
kilocalories of non-flxed energy are available to it. Fourth, the 
locked-rotor transition state becomes more Important relative to 
the orbiting complex as the temperature is raised. Finally, the 
low number of quantum states associated vlth the locked-rotor 
transition state must result in significantly decreased micro- 
scopic rate constants for association and dissociation reactions 
relative to those expected with an orbiting transition state. 

Many incoming reactive configurations would reflect off such a 
barrier. 

Proposed Mechanism for Ion-Molecule Reactions 

We propose that ionic systems take advantage of their Inherent 

ion-dipole and ion-induced dipole attractions, which can allow 

24 

quite substantial energies (10-15 kcal/mole ) of unstructured 

25 

association between charged and neutral fragments. In Figure 5 
we employ an energy contour description of the system 

dynamics. One of the methyl fragments is fixed at the origin, and 
the other is allowed to approach it. The radial coordinate 
r (0 £ r < ») corresponds to intercarbon separation. The angular 
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coordinate. 6(-n/2 < 6 < + rr/2) corresponds to the acute angle 
becveen the line joining the carbon atoms, and the normal to the 
plane described by the atoms of the non-fixed fragment: 
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The approximate potential energy contours were generated by 
using the construction 

V(r,9) - V dlr (r)cos 2 9 + 

where v j ir ( r ) i* Che effective potential function V # ^(r) 
direct CC bond formation, developed in the calculations section. 
V^ d (r) is the appropriate 4-6-12 function for the formation of a 
loose association complex, and is generated in the following way: 

The binding energy of the association complex can be estimated as 

5 kcal/mole, given the polarizability of methyl radical. If we 
model the CH^ fragments as spheres with radii of 1.5 k--CR bond 
length plus the radius of H— then r^ ■ 3.0 A. From these parameters 
a reasonable potential curve is constructed using equation (1). 

Note that Vfr,B) reduces to V d ^(r) for 0*0, and to V^ nd (r) for 

6 - +tt/ 2, -r/2. 
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The ouCer circle et zero energy corresponds Co the orbicing 
transition scace ac a CC separaClon of 5 1; lCs locaCion is 
IndependenC of relative orientation. The molecular cation equi- 
librium CC separation at r ■ 2.0 A shows up as a very deep well 
(V(2.0 A, 9 ■ 0) * -78 kcal/mcle), while the potential minimum 
along a coordinate of constant 6 decreases in magnitude as 9 varies 
toward +n/2 or -rr/2. At 0 ■ -m/2, -rr/2 the minimum of -5 
kcal/molc arises solely through ion-induced dipole attraction, 
with no covalent component to the potential. 

The dashed line corresponds to the approximate location of 
the set of loeked-rotor transition states on the two-dimensional 
surface. The inclusion of this line is based on the assumption 
that the two methyl groups lock up only if (1) the approaching 
fragment is aligned to within a specified angle 9 (9 is estimated 
as rr/4 in the Figure), and (2) the fragment separation is within 
the distance r * 4.3 A, which is the location of the loeked-rotor 
transition state on the CC covalent bond potential profile (9 * 0) 
that we determined above. This assumption is in accord with the 
arguments presented above, namely, that the formation of a loose 
association complex does not lock up the free rotations of the 
fragments. The entropic barrier lies inside a large portion of 
the surface that is stabilized by ion-induced dipole interactions 
to the extent of 5-10 kcal/mole. 

So long as the methyl fragments do not lock up, and are able 
to freely tumble about each other, one can make the reasonable 
assumption that the frequencies and moments of Inertia of the 


fragments remain constant as r is varied, i.e. these parameters 
can be considered to be independent of fragment separation in the 
process of forming a loose association complex. Thus the number 
of states for a loose complex Increases monotonically , and dras- 
tically, vith decreasing CC separation, as the system gains non- 
fixed energy. The fragments can lock easily, vith a negligible 
sum of states bottleneck, if the locking occurs in a region of the 
potential surface where many lccal of non-fixed energy are available 
to the system. 

The trajectory in figure 5 near the 9*0 coordinate corre- 
sponds to direct CC bond formation. It crosses the entropic 
barrier early in the bonding process. On the other hand, the 
second trajectory provides an example of indirect association, in 
which the formation of an unstructured complex precedes covalent 
bonding. This trajectory crosses the barrier later (in the sense 
that the system is at lover potential energy) in the bonding 
process. 

Thi.s we believe that the ion-neutral system is able to access 
a well that does not exist in the neutral system; an ion and 
molecule can avoid the locking up of the free rotors that charac- 
terizes the entropic barrier to association. The neutral and 
ionic fragments approach each other and the system passes through 
the orbiting transition state (which is the operative transition 
state for the process). At this point, however, the system "feels" 
the entropic barrier that would be brought on by the formation of 
a chemical bond and concomitant locking of the free rotors, and 
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avoids that locking by channeling away from the locked-rotor into 
a well representing an unstructured intermediate. This well does 
not exist in the case of neutral -neutral interactions, and the 
locked-rotor transition state cannot be circumvented in the same 
manner. Of course, the ion-neutral system must still undergo a 
locking up of the rotors at some point, regardless of reaction 
path; however, because the formation of an unstructured complex 
has released some 5 kcal/mole of non-flxed energy into the system, 
the ion-neutral intermediate can pass with ease into the well 
representing the structured molecule. This process corresponds 
to a trajectory crossing the dashed line In Figure 5. The locking 
process chat occurs upon forming a chemical bond Is klnetlcally 
unimportant In a potential regime so far below the energy of the 
separated fragments. Moreover, we have treated a rather conserva- 
tive system in the CH^/CH^ case: most loose ion-molecule complexes 

have association energies significantly greater than 5 kcal/mole, 
since larger neutrals are more polarizable than CH^ . 

Thus the locked-rotor transition state has a chance of being 
rate -determining only if the system channels through it relatively 
early In the recombination process, away from any deep potential 
wells. If the system can somehow avoid locking until later lu the 
process, then the locked-rotor state is unimportant. Ion-molecule 
systems can channel away from an early confrontation with the 
locked-rotor entropy barrier through the dynamical mechanism 
described above; neutral systems do not have that same choice, but 
must channel through the locked-rotor state relatively early in a 
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recombination reaction, where the entTopic barrier could be 
significant. 

One can also consider this notion in terms of an ensemble 
of incoming reactive systems, constituting a distribution of 
oulecular orientations and impact parameters. An ion and molecule 
are mutually attracted long before they are close enough to discern 
how they are oriented with respect to each other; they will form 
an unstructured complex regardless of orientation. If properly 
oriented, they can do the reaction directly, but the probability 
of this is so low and the rate of the process so small that this 
pathway is unimportant. 

Note that this model implies certain unusual characteristics 
in the dynamics of the ethane radical caticn dissociation/recombina- 
tion system. The incoming fragments, given appropriate trajec- 
tories to form a chemical bond upon close approach, actually 
"prefer” to avoid forming that bond, and instead tumble about each 
other in a loose association complex prior to chemical structuring. 

By microscopic reversibility, the ethane cation dissociates not 
through simple CC bond stretching, but through some kind of bend- 
ing and subsequent collapse to the same unstructured intermediate. 
This process would be followed by diffusion of the fragments away 

from each other, with the Langevin orbiting complex serving as the 
13 18 

transition state. ' 

Actually, for the system the locked -rotor configura- 

tion does not constitute a severe entroplc barrier; in general, 
the barrier may be large, depending on the particular system treated. 


Both Hase's and our modeling suggest that severe bottlenecks (as 
evidenced by a sharp minimum in the sums of states) occur in the 
regime of a few kcal of non-fixed energy, since more energy than 
chat w:'ll Increase the sum of states greatly. 

Consequences with Regard tc Interpretation of Experiments 

The results of this work have important implications in 
interpreting some experiments in the ion-molecule area. Chesnavlch 
et al_. , in their analysis of C^Hg* decompositions, have pointed 
out that apparent inconsistencies exist between PIPECO (photoion- 
photoelectron coincidence spectroscopy) results that imply barriers, 
and the often observed orbiting transition state rate constant 
for thermal ion-molecule reactions.^ In addition, the product 
translational energy distribution appears to be statistical, even 
in cases where the decomposition reactions are slow. They suggest 
that a locked-rotor transition state can account for the above 
results . 

As we have shown above, the locked-rotor cannot be playing 

an important role in any of the dynamics. Thus, an alternative 

explanation is required. In spite of the statement of Chesnavlch 

et el. that no model has been postulated which accounts for the 

literature results, any double well potential surface will suffice. 

The double well potential has been suggested to account for slow 
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ion-molecule reactions, and it has been noted that its effects 
will be Increasingly important at higher temperatures or energies. 
There is good reason to believe that many reactions which are 
encounter controlled at room temperature will be slow at higher 
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temperature owing to barriers, and in at least some cases there is 
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experimental verification of this. Any potential surface of 
this type will produce statistical translation in the products, 
because the loose complex following the barrier can equlpartitlon 
its energy. The locked-rotor transition state will, of course, 
also produce statistical translation, even if no other barriers 
are present. 

The central barrier (or barriers) will, in general, be chemical 
in nature. For example, rearrangements in the C^Hg + system surely 
involve activation energy (as in solution). It is very unlikely 
that fragmentations which Involve such extensive rearrangements 
could occur on a monotonically increasing surface. An analysis of 
the details of the surxace is not in order here, although there 
may well be a conmon intermediate or intermediates present. In 
any case, the system resembles many others whose dynamics have been 
interpreted in terms of this general model; there is nothing about 
its characteristics which are at variance with the simple double 
well potential picture. 

Conclusion 

In summary, we have shown that the locked-rotor transition 
state occurs late, and has negligible kinetic importance in ion- 
molecule reactions, in contrast to neutral systems where it can 
play a crucial role in the dynamics. Quantitative differences 
between ion-neutral and neutral-neutral interactions produce a 
qualitative change in the chemical mechanism: namely, the 
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> 5 kcal/mole energy stabilization of unstructured ion-molecule 
Intermediates always allows the circumvention of early entropic 
barriers. Potential surfaces incorporating true chemical barriers 
can explain the reaction patterns of ions; in particular, the 
previously suggested double well potential accounts for the 
seemingly anomalous rates of decomposition and product ratios of 
energized ions. In reactions such as dissociation in which there 
is no major potential barrier, centrifugal additions onto the 
zero-angular momentum reaction coordinate determine the distance 
of the (orbiting) transition state. 
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Table 1. Parameters Used in the Calculation of Sums of States. 



Orbiting transition 
state 

. Locked -rotor transition 
State 

v i 

3162 (4) 

2975 

(4) 


3044 (2) 

2915 

(2) 


1396 (4) 

1472 

(2) 


617 (2) 

1460 

(2) 



1388 

(2) 



160 b 

(2) 


• 

lll b 

(2) 

I i* s i 

1.23, 2 (2) 

3.56, 3 

(1) 


1.23, 1 (2) 




3.56, 3 (1) 




moments of inertia in amu k 2 . 3 ^ is the symme try number 

rotation. 

^Values that minimize the number of states; see text. 




Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Potential energy functions for CH^CH^ homolysis, plotted 
against CC bond distance coordinate. Note that account- 
ing for nonzero angular momentum raises the potential 
at all values of r, and creates a maximum in V gff (r) at 
r “ r_ . In the case of ions r^ is the orbiting (or 
Langevln) radius. is the dissociation energy for 
J - 0. 


Potential energy function v e ff( r ) * V 4-6-12^ + 

V t (r) for CH^CH^ dissociation, plotted against CC 
bond distance coordinate. The potential energy zero 
(here and in Figure 3) is taken as that of Infinite 
fragment separation. 


Potential energy function v e ff( r ) a3 shown more fully 
in Figure 2, along with the sums of states function 
G t (r), both plotted against CC bond distance coordinate. 


Variation of the sums of states function G with tempera- 
ture, shown for the orbiting configuration (r “4.95 k, 
upper curve) and for the locked -rotor configuration 

(r “ 4.3 k, lower curve). 

/ 


Energy contour plot shoving the dynamics of the CH^CH^ 
system. See text for discussion. 







r / A 


I 



009 



Chapter VI 
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Reaction Rate Constant for OH + HC0N0 2 — ■ Products over the Temperature Range 
246 to 324 K 


Paula L. Trevor, Graham Black, and John R. Parker* 

Pfryaoa Scanom 0**tv * SW/ ImaraOonaL Marttc Pa*. CaStamW 94025 (rtmaamt Plenary 1Z 1901: 
In nat form Umamtar 2X 1 941) 


Absolute bimolecular reaction rata constants for the title reaction have been determined for temperatures ranging 
from 246 to 324 K. The laser flash- photolysis resonance- fluorescence iLFPRF) technique was used to generate 
0('D) which reacted with H, and/or H,0 to produce OH radicals. The bimoleculai rate constants for the title 
reaction showed no dependence on total (He) pressure over the range —3 to 15 torr, and they did not depend 
upon initial [OH] or upon its mode of formation The H^O? impurity was explicitly measured in ail experiments, 
and the rate constants were corrected for ita contribution A weighted least-squares analysis cf the data obtained 
at nine temperatures (226 data poiuia) gave the Arrhenius expression (h ± l<r) » (8.06 ± 5.69) X l(T i2 exp(-193 
± 194 /T) an 3 s'* with covariance 1.098 X 10*. A simple weighted average (temperature independent) fits the 
data just ae wed, and when the effects ol ryitemaUc errors ere taken into account, our recommended rate constant 
is ( k ± 2<r) - (4.0 ± 1.6) X l(T u cm 3 a* 1 . 


Introduction 

Both HO] and NO? play crucial roles in the chemistry 
of the upper atmosphere. For several years, these two 
species were thought to react via a radical disproportion- 
ation reaction, although it wao suggested 1 that a longer- 
lived complex could be formed according to reaction 1. 


HOj + NO, 4* M cs HOONO, + M (1) 

The importance of reaction 1 has been verified by direct 
observation of H0 2 N0 5 using Fourier transform infrared 
spectroscopy. 3 Moreover, the rate of reaction 1 has been 
measured, 1 * and its reverse reaction (-1) has been studied.* 


(2) K. Niki, P. D. Maker. C. M. Srvage, end L P. Broitenbach. Clam. 
Phyt. ieic. 44. 564 (1977). 


Ill R. Simonaitia end J. Heicklen. J. Phyt. Cham.. 89, t (1976). 
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ngwe 1. Sctiemeflc dtepam of laeer Path ohcuctytk. reeorv w-fir-irss cr cm expnnmer.iei apparatus For experiment aeriee d, a Krr exdmer 
laeer (246 rwn) wi ueed, and for experiment seriea C and 0. Vie quadrupled output (2M nm) of a NdtYAQ laaer was used. 


Also, RRKM theory has been applied to these data lor 
reaction -1 so that its rate can be reliably estimated pj> 
function of both temperature and pressure . 5 In another 
publication, the effect of H0-N0 2 on tropospher >ho- 
tochemical smog chemistry 4 was discussed. 

A concerted effort is underway to investigr e p* l’.tric 
acid (PNA) and determine the rates of various i .» 
me chan isms that affect its residence time in the strato- 
sphere. In addition to the experimental studies reported 
on its tnermal decomposition rate, results of investigations 
on optical absorption cross sections have been reported,''* 
permitting calculation of the photolysis rates appropriate 
for the atmosphere. 

In our laboratory, reactions of PNA with several at- 
mospheric species have been studied to determine bimo- 
Iccul*' rate constants for the destruction of PNA, including 
reactions with 0 atoms , 10 H atoms , 11 and 0 3 . 12 Each of 
these three reaction* is too slow to be important for at- 
mospheric chemistry. 

Of the atmospheric species, the OH radical is probably 
the most reactive, and its reaction with PNA can have a 
profound effect on the chemistry of the stratosphe-e. The 
reaction products are not known, but several sets are 
thermoc heroically possible. 


OH e HOONO, - H O + 0,NO, 

(2s) 

1 


0, a- NO, 


- H,0, + NO, 

(2b) 

- HO, * HNO, 

(2c) 


If the reaction products are given by eq 2a, model calcu- 
lations 11 show that the net result of reactions I and 2 is 


(3) C. J. Howard, J. C/um. Phyt.. «7. 5258 (19TO. 

(4) R. A. Graham. A. M Winar. utd J. N. Pitt*. Jr., Chem. Phyt. Lett.. 
St. 215 usm; J. Chem Phyt.. M. 4505 (1978). 

(5) A. C. Baldwin and D. M. Gold.n. J. Phyt. Chem.. 82. 644 (1978). 

(6) A. C. Baldwin. J. R. Barkar, 0. M. Gulden, and 0. C. Hendry, J. 
Phyt. Chem.. 81. 2433 (1977). 

(7) R. A. Graham. A M. Winer, and J. N. Pitta, Gtophyt. Rn. Lett., 
S. 909 (1978). 

(8) R. A Cox and K. Patrick, /nr. J Chem. Kinet., 11. 635 (1979) 

(9) L. T. Molina and M. J Molina. 14th Informal Conference on 
Photochemiatry, Mar 30- Apr 3. I960. Paper C2; M. J. Molina. Report No. 
F AA- EE-80-0'’. U S. Department of Tranapcrtation. Federal Aviation 
Admmiatratlon. 1980. 

(10) J. S. Chan*, P L Trevor, end J. R. Barker Tut. J. Chem. Kinet., 
13. 1151 (1981). 

(11) P. L. Trevor and J. R. Barker, lu. J. Chem. Kinet.. 13. 1163 
(1981). 

(12) P. L. Trevor end J. E. Dr- enporv. pnvate communication. 


the NOj-catalyred disproportionation of OH and H0 5 
radicals. The effect of this reaction is to reduce the cal- 
culated depletion of the stratospheric ozone coiumn by 
about a factor of 2. 13 Because of the sensitivity of strat- 
ospheric chemistry calculations to the rate and products 
of reaction 2 , we have made a great effort to obtain a 
reliable experimental measurement of the reaction rate 
constant . 14 Unfortunately, we are not yet able to identify 
the reaction products. 

Although the primary reason for studying PNA reactions 
is their potential importance in the atmosphere, PNA also 
is an interesting chemical that is part of the fabric of HNO 
species. In this sense, it is a member of the homologous 
series of nitrogen acids (e.g„ HONO, HONOj, HOONOj), 
as well as being the first member of a series of perm Crates 
(e.g., HOONOj, CH 3 OONO 5 , CH 3 (C0)0N0 2 ; in addition, 
it is a hydroperoxide (e.g., HOOH, HOOR, HOONOj). All 
these species are important for various reasons, and data 
are needed on their reactions. Sucn data can then be used 
as a basis for predicting reaction -ate constants for species 
that have not, as yet, been studied experimentally. 

Experimental Section 

The experiments were performed in four series, each 
differing from the other in various ways, mostly due to 
improvements in technique. Seriea A was performed using 
a flow-tube apparatus and was unsuccessful, but senes B, 
C. and D were performed using variations on the flash- 
photolysis resonance- fluorescence (FPRF) technique cou- 
pled with mass spectrometric measurement of [PNA], and 
they were fully successful. Although the flow-tuDe ex- 
periments were unsuccessful, much of the apparatus is 
common to the other series and the flow tube itself was 
used as a mixing mainfold in the other series, justifying 
a brief description of the apparatus. 

Seriea A: Flow-Tube Apparatus. 1 * The apparatus 16 
consists of two main parts, the flow system and the mass 
specTronetric detection system. The flow system consists 
of a flow tube, a gas-metering system, a mass spectrome- 
tric-sampling device, and a high -capacity pump. As shown 
in Figure 1, the main flow tube is a Pyrex tubo, 1 m long 
and 2.5 cm i.d., equipped with a manifold to prepare the 


(13) la) R. Cicerone and S. Wallen. 14th Informal Conference on 
Photochemistry, Mar 30-Apr 3. 1980. Paper C8; (b) W. H. Duewer and 
D. J. Wuanbtaa. ibid.. Papa;- J6: (c) J. Steed, pnvate communication; Id) 
S. C. Liu. private communication. 

(14) Preiimi.ery reauite for tha flow-tuba date - ere presented earlier 
P. Trevor, J. S. Chen*, end J. R. Barke-. 14th Informal Conference on 
Photochemistry. Mar 30-Apr 3, 1980, Papsr C3. 

(15) J. S. Chan* and J. R. Barkar, J Phyt. Chem.. 83, 3059 11979). 
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Reaction Rate Constant for OH + HOONO, — Product* 

main flow mixture before it if introduced into the flow 
tube. Typically, a small amount of hydrogen mixed with 
helium constitutes the main flow which ia then passed 
through a microwave discharge (2450 MHz) before being 
introduced into the flow tube. A Wood’s horn ia used to 
prevent light generated by the microwave discharge from 
entering the main flow tube. The He flow is controlled 
with a needle valve ard ia measured by a Fisher-Porter 
rotameter, reactant flows aro measured by the pressure 
changes in a calibrated volume fitted with a Validyne 
pressure transducer (±0.1 psi). 

Sampling of the reactants and the products from the 
flow tube to the mass spectrometer is accomplished by a 
Teflon cap that is located at the downstream end of the 
flow tube. In the center of the Teflon cap, a 0.010 in. 
diameter pinhole, situated close to the center of the flow 
tubs, ssrves as a molecular beam source; the beam ia 
modulated by a mechanical chopper and detected by the 
mass spectrometer. The pressure in the first pumping 
chamber after the pinhole is typically ~ 1-3 X 10"* torr, 
while the pressure in the second chamtwr (containing the 
Finnigan quadrupole) is ~1(T* torr. The same modulated 
molecular- beem mass spectrometer (MMBMS) has been 
used for detection in various other studies, such as very 
low-pressure pyrolysis and very low-pressure photolysis, 
and a description can be found elsewhere. 1 * 

Series P: 248-nm Flash- Fhotoly sis Retonance- 

F-uorescence ( FPRF ). The apparatus used in experiment 
series B, C, and D is shown schematically in Figure 1. For 
this setup, the flow tube described above was coupled to 
a resonance- fluorescence cell as shown in the figure. The 
cell is not the same one used for series A and is 70 cm on 
the long axir to reduce scattered Light from the CaFj 
windows that transmit the laser beam. The top and side 
windows are quartz to transmit the OH radical reso- 
nance-fluorescence light centered near 3100 A. Light from 
the OH resonance lamp (2450 MHz; 9 torr He passed over 
HjO) is collimated or weakly focussed with a quartz lens 
and directed through the side window of the cell. The flow 
tube is coupled to the cell directly opposite the side win- 
dow. Since tne long axis of the cell is oriented at ~35° 
from the axis of the flew tube, the light from the resonance 
lamp does not pass do.vn the length of the cube out is 
stopped by the bend of the coupling tubulation. This 
arrangement acts as a crude Wood’s horn to reduce scat- 
tered light from the amp. 

The photomultiplier tube (EMI Model 6255B) was wired 
for single photon counting and viewed the OH fluorescence 
through the top window of the ceii. The 3100-A inter- 
ference filter was used to reduce scattered light. Output 
from the photomultiplier was amplified and an amplifi- 
er/ discriminator was used to detect single-photon pulses. 
The output of the discriminator was stored in a multi- 
channel analyzer (Nuclear Data, Model ND100) that was 
operated in the multichannel scaling mode, where each 
channel corresponds to a time increment After several 
thousand laser shots, the memory contents were recorded 
on on z-y chart recorder and were then analyzed to de- 
termine the fluorescence decay rate. 

The gas flow entered the cell at the port opposite the 
side window and then turned 90° downward through the 
exit tube that is opposite the top window. About 8 cm 
from the optical center, the gas stream was sampled 
through a Teflon cap, perforated with a hole —0.008 in. 
in diameter that acted as the originating orifice of the 
modulated molecular- beam mass spectrometer described 


(16) J. S. Ch*n(. J. P. Barker. J. E. Davenport, and D. M. Golden. 
Chtm. Phyi. Lett.. SO. 386 ! 1979). 
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above. The gas flow pattern in the fluorescence cell is 
somewhat complex (sen Figure 1) but probably consists of 
a rapid, turbulent flow in the optical viewing region, but 
nearly stagnant behavior in the long side arms. The mass 
spectrometer sampling orifice ia only about 3 cm directly 
downstream from the optical center and appeared to give 
a reliable measure of gas composition at that center, as 
demonstrated with test reaction!. At the exit port 01 the 
cell (see Figure 1), a 10- mm glass or Teflon-bore stopcock 
fitted with Viton O-rings was used aa throttlo valve to 
control gas flow velocity and total pressure. From this 
valve, the gas flow was conducted directly to the main 
flow-tube pump. 

Temperature control of the cell was maintained by 
circulating fluid irom a constant- tamperture bath (Neelab 
LT-50) through coils of tubing wrapped around the cell. 
A sheet of closed-cell foam Insulation covered the coi*s and 
five chromel-conatantan thermocouple* attached at var- 
ious points around the call body showed the temperature 
to be constant and uniform witnin ±2 °C during the ex- 
periments. 

In these experiments, the flow tube was used primarily 
as a gas- mixing manifold; lines* flow velocities were usually 
of the order of 0.5-3 m/s. but faster valociti is were used 
occasionally. Concentration calibrations for NO* Hj, and 
HNO-) were performed by measuring the rata of pressure 
increase in calibrated volumes, as described previously. 17 
Dry HNO, was prepared by mixing NaNO, with concen- 
trated HjSO, in an evacuated bulb. Ozone concentration 
calibrations were performed daily by measuring the at- 
tenuation 1 * of light to a Hg 2537-A pe alight source 
viewed through an interference filter by a RCA 1P28 
photomultiplier down the long axis of the photolysis celL 

Ozone was prepared in a separate vacu'im line by an 
electric discharge in 0 } and wsc trapp-.d and stored in a 
silica gel-packed trap maintained at -78 °C. The ozone 
concentration was controlled by sweeping it cut of the trap 
with a stream of helium diverted from the main flow. Maas 
spectrometric analysis showed that HjO vapor waa released 
from the trap along with 0 5 and Oj. 

Pulses of OH radicals were generated according to the 


following scheme 

Oj + M248 cm) — 0(‘D) + 0 7 (3) 

0(‘D) + H, — OK + H (4) 

O('D) + PNA — • OH + 0, 8 NOj (5) 
H + Oj -* OH(vib) + Oj (6) 


The 2 13- am light was generated by a Lambda-Physik 
ezeimer laser (Model EMG 101) operated with a mixture 
of krypton, fluorine, and helium. Output energies were 
measured with a Scientech calorimeter-type power meter 
(Model 360203) and averaged ~80 mJ/pulse, of which 
energy about 25% actually entered the cell after collima- 
tion. Typically, pulse repetition frequencies used were 
-5-10 Hz. Since gas residence times in the cell body were 
— ‘/j s, the Oj was significantly depleted by the laser beam, 
but the PNA concentration showed less than ~5% de- 
pletion, demonstrating the validity of pseudo- first-order 
conditions. It was found by experimentation that the 
oboerved OH fluorescence decay rates were independent 
of Oj pressures (at this Laser pulse energy ' as long as they 
were less ihsn ~ 1 mtorr, corresponding to an initial OH 


(17) R. L Ksnlsv. P. L Trevor, ind B. Y. Ian. J. Am. Chem. So c.. 103. 
2207 (1981). 

(18) M. Griggs, J. CStm. Phy.. 4», 857 (1968). 
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concentration of 5 2 X 10 u radical* /cm 3 . Hydrogen con- 
centration* were high enough (2 X lu is molecules/ cm 3 ) so 
that all of the 0(*D' reacted by reaction 4 within — 10 *is 

after the laser pulse 

Ti e troublesome possibility existed that vibrationaily 
excited OH was present initially, and an observed spectrum 
showed our resonance lamp could excite OH(o * 1), if 
present. Thus, the fluorescence data for hmee shorter than 
15 ms after the laser pulse were ignored to allow some time 
for vibrational relaxation of OH(vib). Several spikes due 
to electronic pickup from the Laser appeared on each dray 
curve, providing a convenient time marker for the 15- ms 
delay. 

In this series of experiments, it wss found that PNA 
decomposition on the Pyrex surface could be somewhat 
inhibited by application of a thin coating of haiocarbon 
wax (the preparation and absolute concentration mea- 
surements of PNA are discussed below). This treatment 
permitted 'he incorporation of a cold trap in the PNA 
preparation train that helped to reduce unwanted HjOj. 
Moreover, higher total pressures were obtainable without 
excessive decomposition cf the PNA. The primary high- 
preasure limitation, however, waa the residual gas pressure 
tolerable for the mass spectrometer. 

Series C: 266- nm FPRF. 70 This series of experiment* 
was performed in a manner similar to series B, with a few 
differences. The light source for this series was a Quanta 
Rey Nd:YAG laser whose output was quadrupled in fre- 
quency to give 266- nm light (~20 mJ/pulse, 90% of which 
entered the cell). The scarce of OH was 

0, + hr) 266 nm) — 0( ; D) + 0 : (3) 

0( l D) + HjO -* 20H (7) 

This source of OH was “cleaner" than in series B, because 
the initisi vibrational excitation in OH produced by :e- 
action 7 is quits low, n and there is no subsequent reaction 
of H + Oj to produce highly vibrationaily excited radicals. 
Due to the absence of vibrationaily excited radicals, the 
data analysis could be carried out for all data from t = 0. 
The initial OH concentration is estimated to be £ 2 X lCP 
radicals/ cm 3 for 1 mtorr 0 3 . 

A new quadrupole mass spectrometer (Balzere Models 
311. 140) was substituted for the Finnigan instrument and 
a new series of relative mass spectrometer sensitivity 
measurements was performed. 

One additional change was made in the data analysis 
procedure. The trigger sequence was arranged so that 
background fluorescence was recorded for several tenths 
of a millisecond before the laser was pulsed, enabling us 
to digitally subtract the background light intensity. After 
subtraction of background, the memory contents of the 
MCA were passed through a logarithmic amplifier prior 
to the x-y recorder. These steps considerably reduced the 
time consumed in the data analysis process. 

Series D. 2t€-nm FPRF. This aerie* of experiments was 
performed in just the same way as series C, except that 
the lens elements of the mass spectrometer were cleaned 
and the electron multiplier waa rejuvenated (because of 
an accidental vacuum failure), necessitating a new series 
of relative sensitivity measurement*. 

PNA Preparation, Purity, and Concentration Deter- 
mination. The preparation, purity, and concentration 


119) J. A. D«v\d*oo. H. L Sctuff, G. E. Stmt. J. R. McA/m. A. L. 
Scimwluiopf. tad C. J. Howard. J. CKtm. Ph y«.. (7. 5021 (1977). 

1 20 1 Praiiminarv rtaulta ofacaioad in apartment Sanaa 0 wara reported 
aerliar 1. R. Berker. G. Black and P. L Trevor. ’980 American Geo- 
phyi'cai Unrna Pall Maeon*. San Fraacuco. CA. 8-12 December 1980. 
(2)) K. H. Gencka aad F. J. Comae. Diem. Phyt. Lett.. !i. 83 (1980). 
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determination of PNA require a separate section because 
this aspect of the experiments is the single largest source 
of error and uncertainty. 

Our experience is tnat PNA tends to decompose rela- 
tively rapidly (seconds) on Pyrex and metal surfaces. 
Haiocarbon wax and Teflon tend to be more inert, but 
some decomposition still takes place. Due to this behavior 
by PNA, we did not use static mixtures of PNA in storage 
bulb*, but generated the PNA batch wise and mixed it with 
the gaa stream in the flow tube. Small preparative batches 
of PNA are advised because, during this study (CAU- 
TION!!) two minor explosions took place, one of which 
destroyed a heated glass tube, even though the partial 
pressure of PNA was probably less than 5 tore and the total 
pressure was less than 20 tore. 

As described elsewhere, PNA was prepared 17 by cau- 
tiously and slowly adding 0.6 g of NOjBF 4 (Alfa-Ventron) 
to 90% HjO] maintained at 0 °C and rapidly stirred. This 
preparative solution waa kept at 0 °C until the smell vessel 
(—25 cm 3 ) was O-ring sealed to the preparative train and 
pumped down to the flow tube pressure. PNA was carried 
out of the preperative vessel by a stream of helium (di- 
verted from the main flow) that was passed over the liquid 
surface. Bubbling the carrier gao through the solution did 
not significantly improve the PNA yield. The flow-tube 
PNA concentration was controlled by varying the helium 
flow rate and by varying the bath temperature of the PNA 
preparative solution up to —15 °C; it was deemed too 
dangerous to allow the temperature to go higher. 

In experimental series B, C, and D, a U-trap followed 
the preparation vessel, and when heid at a temperature 
of - — 15 °C, it seemed to reduce the relative amount of 
impurity H 2 0 2 that was released. 

Potential impurities, such as N0 2 , H 2 0 ? , and HNOj, 
were specifically sought using the mass spectrometer, and 
it was found that HN0 3 either was not present or its level 
wss <5% of the PNA (see below). Hydrogen peroxide 
appeared as a highly variable impurity, ranging from un- 
detectable levels ( m/e ■ 34) up to concentrations several 
times as large os the PNA concentration. NO- waa often 
present as an impurity, but its level was usually £10%, 
and never more than —50% of the PNA concentration. 
As discussed later, the reaction rate of CH with N0 2 at 
the total pressures of He used here are slow, and thus the 
effect of N0 2 is small. H 2 0 and 0 2 were also impurities, 
but they were not monitored. 

The basic ides behind the relative mass spectral sen- 
sitivity measurements for PNA and N0 2 is that the 
thermal decomposition of PNA quantitatively yields N0 2 

dot 

HOONO, — ' HOO + NO, (8) 

2HO, •* H,0, + O, (9) 

H,0 + '/, O, 

The experimental procedure was as follows. For a 
constant flow of PNA: (1) Measure the mass spectral 
intensity at m/e ■ 46: / ( (millivolts). (2) Using the 
flow-tube microwave discharge, add 0 atoms and titrate 
away any N0 2 impurity using the N0 2 * afterglow tech- 
nique; measure the mass spectral intensity: /-. (3) With 
the 0 atom concentration from step 2, divert the PNA flow 
through the heated bypass [a ‘/<-in. X 50 cm stainless steel 
(because of explosion hazard' U-tube maintained at 300 
°C], and measure the mass spectral intensity: I y (4) Add 
more 0 atoms to titrate the N0 2 produced by the thermal 
decomposition of PNA, and measure the mass spectral 
intensity: /«. (Titration of N0 2 by 0 atoms is an effective 
technique, because the rate constant O 4- N0 2 is about 


Reaction Rat* Constant for OH + HOONO, — Products 

TABLE I: Relative Maas Spectrometer Seuaitivitias 

for PNA and NO, 

data mass 

exptl series S t Id points spectrometer 

B 7.48 t 1.05 21 Finnigan 

C 4.44 t 0.74 9 Salzera 

D 5.07 t 0.34 14 Balzeri 

three orders of magnitude faster than that of the reaction 
0 + PNA. As long as a modicum of care is taken. PNA 
is not depleted by reaction with 0 atoms under these 
conditional. 

The difference (7j - /*) is proportional to the amount 
of NO] produced by the thermal decomposition of PNA. 
The difference (7, - IJ is proportional to the amount of 
PNA decomposed. The ratio of these quantities is the ratio 
of the mass spectrometer sensitivities for the two com- 
pounds 

(7, - 7*)/(7, - IJ - S^ 1 (10) 

and [PNA] “ 5 -1 £r NOt -7 1) where a^o, is the calibrated NO, 
sensitivity (molecult* cm' J mV' 1 ). Since a small fraction 
of the PNA can survive the heated tube, I t is usually not 
identically zero; the portion of PNA that survived de- 
pended upon flow conditions and ranged from nearly zero 
up to ~2Q%. 

To determine whether HNOj is an impurity originally 
present in the PNA or is produced in its thermal decom- 
position and contributes to the obeerved mass spectral 
intensity of m/e ■ 45, the PNA ie diverted through the 
heated bypass and a great excess of 0 atoms (i.e., >10 u 
atom/cm 3 ) is added to remove the NO, and the PNA that 
survived the heated tube. Under theee conditions (reaction 
time ~ 1 s), it is obeerved that virtually all of the intensity 
at m/e " 4€ disappears, indicating the absence of HNO, 
since HNO, does not thermally decompoae under these 
conditions, and it reacts only very slowly with 0 atoms. 
On the other hand, PNA both thermally decomposes and 
reacts with O atoms at a moderate rate (k ~ 10" 14 cm 3 s' 1 
at 300 K). 10 All experiments of this type that were oer- 
formed showed negligible amounts of HNO,. 

It was concluded that the thermal decomposition of 
PNA under our conditions quantitatively gives NO, for the 
following reasons. (1) HNOj was shown not to be a sig- 
nificant product. (2) The ratio of intensities at m/e m 46 
and at m/e m 30 were coua-stent with product NO,, al- 
though the experimental scatter was large due to PNA 
interference and due to taking the differences of relatively 
noisy values. (3) Experiments in which mixtures of NO, 
and HjOj were passed through the heated bypass showed 
that the NO, was unaffected by the treatment. (4) In the 
highly oxidizing environment containing HOONO, and 
H,0,, it seems unlikely that the PNA decomposition 
product would be HONO or other reduced forms of ni- 
trogen. 

The results obtained for the relative sensitivities of NO, 
and PNA in experiment senes B, C, and D, are presented 
in Table L The uncertainties quoted include effects from 
aii sources. 

Results and Discussion 

The conventional flow-tube resonance fluorescence ex- 
periments (series A) were performed under conditions 
where pseudr-first-order conditions obtain. In the course 
of tbs experiments, it was found that an inordinate amount 
of time was necessary for the signal to stabilize after 
moving the probe position. Moreover, there was a large 
amount of scatter in the data, perhaps symptomatic of wail 
reactions involving OH, PNA, or both species. 



0 1 .0 2.0 


TIME I mil 

Ftgiw* 2. First-order plot of tog /, vs. time. Rats data giving k 1 - 
1155 s*' art for [PNA] ■ 1.06 X 10 u moiacutas/cm 2 an) [HyO,] » 
1.05 rWvolts. Data giving k> - 2531 a' 1 are for [PNA] - 1.51 X 
10 U m ofe ci ie s /cm 3 and [H 2 Ojj “ 21.9 mWvofla (too fxgfi for ua# In 
second-order analysis). Both seta of data are from experiment senes 
C. T - 324 K. 

After many efforts to eliminate the anomalous scatter 
in the data, we have concluded that this behavior may be 
due to a wall reaction of CH and PNA that ia dependent 
on the PNA concentration and/or H,0, concentration and 
the history of the flow tube. Halocarbon wax is among the 
least likely substances to cause adsorption of PNA and 
thus is one of the best for inhibiting wall reactions. Since 
it has turned out that this coating is not effective, we 
concluded that further experiments usinr the flow-tube 
technique on this reaction would not be fruitful, and we 
turned to the FPRF approach. 

For the FPRF technique, the OH concentration change 
with time is conventionally written 23 

d[OH]/df - -(If, [PNA] + V[° H 1 (ID 

and the inverse decay lifetime of the fluorescence intensity 

7, is 

r' 1 - A, [PNA] -f- lip (12) 

where A, is the b (molecular reaction rate constant and A„ 
is a parameter that describes the effective ‘rate constant ‘ 
for loss of OH due to physical processes. At high pressures. 
A., may be identified with the fint-ordcr diffusion process 
alone. For the particular conditions in the present ex- 
periments, the physical processes include the translational 
motion of gas parcels as they pass into and out of the 
intersection of the detector field of view and the region 
illuminated by the resonance lamp. Such motions are due 
to the gas-stream velocity as it is pumped through the cell 
and to mass diffusion of [OH], Neither of these processus 
is exactly represented by the first-order ‘rate constant" 
A p since the pumping velocity is clearly independent of the 
OH concentration, and spherical diffusion subsequent to 
excitation in a cylindrical profile coupled with crosaed- 
cptical-beam detection is not expected to be so simple. 
These effects are exhibited in that A„ appeared to be time 
dependant, resulting in OH fluorescence decay plots that 
deviate slightly from exponential decay (Figure 2). 
Nevertheless, eq 12 is a good, pragmatic representation as 
long as the expression is employed for the local slope at 
the same average elapsed time followed the laser pulse and 
the time dependence of A p is not too pronounced. 


1 22) For lumpli. om rtf 24b. 
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Flgwe 3. Seoond-ordar pfc* of v*. [HjOj from experiment hum 
C. 324 K (unoartalntfM shown are 1<r). 



Flgwc 4. Saoone-crttr plot of fr‘ va. [NO,]. The effoOfve MmofeaJw 
rate constant ri *“ - (9.23 ± 0.22) X 10"'* cm 3 a' 1 (uncertutUM 
shown are la). 


Under pseudo-first-order conditions, a plot of log I t vs. 
time will give a nearly straight Line, if k p is only slightly 
dependent on time. Two such plots are shown in Figure 
2 for experiment series C. For series C and D, the protocol 
was established that tho initial slopes on first-order plots 
be taken; this protocol is supported because second-order 
plots of r' 1 vs. reactant oncentrations give good straight 
Linns. Two such examples are presented in Figures 3 and 
4 for the inaction 

OH + H a Oj — H,0 + HO, (13) 

and for 

OH ♦ NOj + He — HONOj + He (14) 

Both second-order plot* show very good straight-line be- 
havior. The absolute rate constant for reaction 13 was not 
determined because the mass spectrometer could not easily 
be calibrated for sensitivity to H,0 2 . (This is due to the 
unknown extent of H?0? decomposition in the inlet lines. 
In subsequent analysis, we have assumed that the mass 
spectrometer signal at m/e m 34 is linearly proportional 
to the unknown H 2 0 } concentration. This assumption is 
consistent with our experience using mass spectrometers 
for detection.) The absolute rate constant for reaction 14 
determined from the data in Figure 5, is 1.6 X If® cm* 
molecule' 2 s' 1 , in excellent agreement with literature data.' 3 
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Figure S. Second-onW. plot of * va. [PNAj for e xp eriment senee 

C,r«2MK. 

TABLE □: Teat Reactions and Result* ot Present Study 



reaction 

T, K 

P, tore 

fcE obad. 

OH 

4. <C 4 H„. 

298 

9.2 

3.6 x 10'" 



324 

10.1 

(3.62 t 0.40) x lO"* 



267 

10.1 

(2.70 i 0.20) x 10"* 

OH 

+ NO, + Ha 

246 

16.1 

(9.23 t 0.22) X lO ' '• 


Similar teats were performed using isobutane as reactant, 
also leading to rate constants in excellent agreement with 
literature values. 24 The teat teaciion rates determined are 
summarized in Table II. These results show that con- 
centrations monitored S cm downstream from the optical 
center are accurate. 

For experiment series B, the protocol for data treatment 
was somewhat different for two reasons. First, because it 
was felt that OH(vib) could cause problems, the first 2.5 
ms of data from each run were ignored, and the first-order 
plots were constructed from data beginning at 2J> ma after 
the laser pulse. The second reason for this procedure waa 
that two or three ‘spikes' due to electronic pickup from 
the excimer Laser appeared in each run, causing some in- 
terference. The last spike occurred at 2.5 ms after the laser 
pulse, providing a convenient time murker. 

Because the data for series B correspond to lower OH 
concentrations (due to decay after the pulse), the 
fluorescence data were significantly noisier than those of 
series C and 0, but they, again, gave good straight sec- 
ond-order plots and absolute rate constants in good 
agreement with the literature. The excellent linearity 
exhibited in Figures 3 and 4 and the excellent results 
presented in Table □ support the use of the approximately 
accurate eq 12 for the second-order plots 

In all these series — B, C, D — the intercept of the sec- 
ond-order plots can be identified with k v It waa observed 
that kp is larger at lower pressures, as expected for diffusive 
loss, and it seemed to depend upon day-to-day variation 
in optical alignment and linear flow velocity through the 
optical intersection region, as expected. Quantitative es- 
timates of linear flow velocity through the ~ 1 cm optical 
center are difficult due to the geometry of the cell, but they 
are probably about the same as in the flow tube (50-300 
cm/s). At the low pressures of the present experiments, 
OH diffusive loss rates are probably of the order of 50-300 


(23) (a) For a summary, s m R. F Hampsoi.. Jr., and D. Garvin. SBS 
Sptt. Publ. (US.), No. 513 (19T7); (bl P R Win*, N. M. Kreuttar. and 
.V R. Ravuhankara. J Ph yt. Cktm., S3. 3191 119791. 

241 (ai N. R. Grtinar. J. Chtm. PH yi„ S3, 1070 119701; (b) K. R. 
Dareall. R. Atkmaon, and J. N Pitta, Jr., J. Ph y». Chtm.. 82, 1581 1 1978), 
and r*fcr*ncaa citad therein. 
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[PNAl x 10 - ’ 3 — molxcutrs'cm 3 

Flgura «. Second-ortw plot tf v*. [PNA| for experiment serine 

C. r - 2*8 K. 


s' 1 . These values are consistent with the observed values 
for k~ 

When H,0j and NO, are present as impurities, eq 12 
must be generalized to 

r- 1 - J^PNA] + * u [He][NO,] + *u[H,Ot] + (151 

As mentioned earlier, the NO, impurity is relatively small, 
and at the pressures and temperatures of the present 
study, its effect on the observed r is also small and can be 
neglected (see below). Since k p varies on a day-to-day 
basis, the beet way of combining all of the data for a given 
temperature and pressure is to consider only the net 
first-order rate constant resulting from the following ex- 
pression 

* kj[PNA] - r-‘ - - kj, (16) 

For experimental runs performed using just HjG* we have 

^*HfOi * fcnlHA] + k p (17) 

For each sot of conditions urns were carried out using only 
HjO, as reactant, and a linear least-squares analysis gave 
values for k l3 and k„. Since H,0, concentration was 
measured by the MMBMS as peak height (millivolts) at 
m/e » 34, has units of (mV)' 1 s' 1 ; is obtained in units 
of s' 1 . Once k p and k l3 are known. H,0j and PNA con- 
centrations are monitored mass spectrometncally and re- 
corded along with r' 1 . From these quantities, is de- 
termined, and second-order pkna of vs. [PNA] can be 
constructed as shown in Figure 5. 

The data for 2S8 K presented in Figure 5 were obtained 
on several different days for three different total pressures 
and several different Linear flow velocities. Clearly, the 
data are quite consistent with one another, showing no 
significant dependence on total pressure. For the purpose 
of illustration, the data for 246 K are presented in Figure 
6, and again, there ia no significant dependence on total 
pressure. 

From the 0, concentration, laser energy, and ozone 
absorption coefficient, 11 we can estimate that the initial 
OH concentration does not exceed —2 X 10 u radicals /cm 3 . 
Several tests weie performed to determine whether the 
observed fluorescence decay rates depended upon either 
the source of OH radicals or on their initial concentration. 
In experiment series B, the O, partial pressure was varied 
from 0.23 to 0.8 mtorr with no perceptable variation in 
either decay rate or in the shape of the fluorescence decay 
curve, although the absolute fluorescence intensity varied. 
When the 0, partial pressure was raised to 3 mtorr. a faster 
decay rate occurred at early times, but by the 2.5 ms time 
marker, there was little effect. 


TABLE OI: Rate Constant* for OH + PNA -* Product* 



daui 10‘ : * 

10‘* 

10"- 

10** 

T, K 

point*® cm 1 ;* 

»FPRF a NOy PNA 

a n«t 



Serie* B: 

248-nm 

F??F 


251 

11 

3.84 

0.28 

0.54 

0.61 

273 

17 

4.60 

0.41 

0.66 

0.77 

298 

23 

6.05 

0.52 

0.84 

0.99 

316 

12 

4.S5 

0.31 

0.69 

0.76 



Serie* C: 

266-nm 

FPRF 


246 

42 

4.80 

0.46 

0.80 

0.92 

267 

22 

4.83 

0.21 

0.81 

0.84 

298 

58 

3.61 

0.16 

0.60 

0.62 

324 

20 

3.79 

0.21 

0.63 

0.66 



Sene* D: 

266-nm 

FPRF 


254 

21 

3.24 

0.30 

0.22 

0.37 

« Total = 

226 data point*. 




In experiments series C, more extensive testa were 
performed. With 30 mtorr H,0, the 0 3 partial pressure 
was varied from 0.4 to 4.0 mtorr, resulting in a variation 
of ~10% in the OH + H,0, bimolecular rate constant 
The 4 mtorr 0 3 data for OH + PNA may have been ~15% 
lower than the usual, but this variation is not significant 
and the data were included in the overall data set 

In another aeries of test runs in series C, H } was sub- 
stituted for H^O, to reproduce the conditions of experiment 
series B. Subeitution of 85 mtorr H, for the usual H,0 
gave no significant difference in observed rate constant 
Indeed, runs with no H, or H 2 0 added to the gas stream 
also agreed with the rest of the data. This result can be 
understood, since reaction of CK’D) with PNA may give 
OH and since H 3 0 was observed evolving both from the 
PNA preparation train and from the 0 3 trap. Thus, it 
appears that H,0 was always present in significant 
quantities. Since HjO + 0(‘D) gives relatively ‘cool’ OH, 
and because H,0 ia an efficient auencher for vibration&lly 
excited OH, it appears that if any residual vibrational 
excitation remained, it played no significant role in the 
observed rate constants. 

In addition to varying the 0 3 , H 2 0, and H } partial 
pressures, the laser power was also varied in experiment 
series C. When the laser pulse energy was reduced by 
about a factor of 10 from the usual 20 mJ / pulse, there was 
a reduction in fluorescence intensity, but no significant 
variation in observed rate constant. Thus, we may con- 
clude that (1) secondary reactions of OH and reaction 
products are not significant, (2) depletion of PNA con- 
centration ia not significant, and (3) OH disproportionation 
ia not significant. Thus, the system appears to be kinet- 
icaily well behaved. 

The measured bimolecular rate constants for each tem- 
perature are summarized in Table HI and Figure 7, along 
with least-squares and propagation of errors estimates of 
the associated experimental uncertainties. The first un- 
certainty trrRT listed for each rate constant is just the 
statistical uncertainty associated with the second-order plot 
of k 1 ^ vs. [PNA]. The second uncertainty listed is that 
associated with the PNA/NO- relative mass spectrometric 
sensitivity. The third uncertainty listed is the propagation 
of errors combination of the? * two sources of error 

A major contributor to fff-PRr is the uncertainty associ- 
ated with the qusntity (k l3 lH 2 OJ + in eq 16. For small 
values of [Hj 0 2 ], this uncertainty is quite small, but for 
a few runs, was very large, posaibly resulting in 

large uncertainties in k l M . To avoid this problem, the few 
runs giving < 0.5k l3 [H 2 O 2 ] were neglected, although 
they would not have affected the final bimolecular rate 
constant. It should be pointed out that pseudo-first-order 
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conditions were assumed in the analysis, but some runs 
may have had a [PNA]/[0H] ratio of less than 5, intro- 
ducing some error into the analysis. B :auae the deter- 
mination of k :s and h v introduce correlation into the runs 
performed on a given day, it is possible that slight errors 
are introduced into the rate constants derived from the 
least-squares analysis. We feel that these errors must be 
very small, however, since rate constants derived from 
experiments performed on different days usually she w 
good consistency. A potential source of systematic error 
in the derived rete constants is the behavior of the 'static* 
gas in the long side arms of the cell. The magnitude of 
error from this source is not known but we feel it to be 
small, since variation of linear flow velocity and pressure 
did not significantly affect the derived bimolecular rate 
constants. 

The effect of any NOi impurity can be shown to be 
small. In the wont case, for our lowest temperature (246 
K) and highest pressure (15.1 torr), the effective bimole- 
cular rate constant for reaction of N0 } with OH is *i«[He| 
■ (9.23 ± 0.22) x 10" 13 cm 3 s' 1 (F’gure 3). If. for our wont 
case estimate, the N0 S impurity is 50% of the PNA con- 
centration, the relative contribution of N0 2 to the observed 
rate constant is less than 10%. Since this wont case 
scenario is intended to be pessimistic and the effect is 
smaller than the other sources of error, effects due to the 
small N0 5 impurity have been neglected. 

The data from Table III show very little variation with 
temperature, but a weighted least-squares calculation 35 of 


(25) R. J. Cvetnnovic. R P Ovirrad. and G. P aroekavopouloua. Int. 
J. Chtm. Kirn 1, 249 (19751; R J. Cvetanovi:. D. L. Stnflaton. and G. 
Paraahevopuuioa. J. Phyt. Chtm , S3. 50 (19791. 


TABLE IV : Rate Constant Companions 

reaction rate constant, cm* s' 1 ref 


OH ♦ H.O, 

(2.51 t 0 6) x 10-“ 

29 


exp(— 126 t 76/T) 



(2.96 t 0.50) x 10'“ 

30 


exp(— 164 t 52/T) 


OH + HOONO, 

S3 x 10'“ 

4 


(4.0 t 1.6) X 10'“ 

present work 


independent of T 



(2TJ) x 10'“ 

26 


4 X 10"’ 

31 

OH * HONO, 

(1.52 t 0.38) X 10' 1 *- 

27 


exp<649 t 69 IT) 


OH HONO 

7 X 10'“ 

32 


the Arrhenius parameters was performed, giving (k ± la) 
- (8.05 ± 5.69) X 1(T U exp(-193 ± 194/ T) cm 3 s' 1 , co- 
variance “ 1.098 X 1CT*; at the center of our temperature 
range, these parameters give a one standard deviation 
uncertainty of about 8%. 

Because the temperature dependence of the rate con- 
stant is so small, a simple weighted average of the nine rate 
constants in Table III is of equal usefulness. Such an 
ai'erage gives {k x Id) » (3-97 ± 027) x HT 13 cm 3 s' 1 . Since 
the uncertainties quoted refer only to statistical fluctuation 
and do not include systematic errors, our recommended 
rate constant and 95% limit uncertainty is (k ± 2a) m (4.0 
± 1.6) x 10" 11 cm 3 s' 1 , where the uncertainty is assigned 
by a subjective estimate of systematic error. 

Our recommended value of (4 0 ± 16) x 10~ 13 can be 
compared to the “upper limit" rate constant of ~3 X 10" 12 
estimated by Graham, Winer, and Pitts 4 on the basis of 
a competitive rate measurement in a smog chamber. More 
recently, an excellent molecular-modulation spectroscopy 
study was reported by Littlejohn and Johnston* that gave 
a rate constant fc 2 * (2 **) x 10* 13 . This last result is 
noteworthy since it depends upon the accuracy and com- 
pleteness of a reaction mechanism consisting of more than 
a dozen reactions, and yet the result is in good agreement 
with the present study. Moreover, the PNA in that study 
was prepared by an entirely different technique and the 
reaction mixture was analyzed spectroscopically, rather 
than by mass spectrometry. Thus, the agreement of these 
two very different studies is quite satisfactory. Further 
experim nts using the molecular modulation technique 
may give information on the reaction products as well. 

In Table IV tne present results are compared to rate 
constants measured for H 2 0 2 , HONO, and H0N0 2 . The 
most perplexing of these rate constants is the recent re- 
sult 27 for OH + HNOj, since it exhibits a "negative acti- 
vation energy." (The Arrhenius form of data expression 
is probably net appropriate for such reaction rate con- 
stants, but it is convenient in the present context.) Al- 
though the reaction products are not beyond question, 
indications art that it is a simple abstraction reaction 35 
making interpretation of the reaction dynamics even more 
obscure (a displacement reaction could be explained on 
the basis of long-range chemical attraction). The OH + 


(26) D Littlejohn and H. S. Johnaton. 1980 American Geophysical 
Union Fall Meet inf, San Franriacc. CA. Dec 8-12. 1980. 

(27) P H. Wine. A. R Raviehankere, N. M. KreutUr. R C. Shah. J 
M. N i cum c h. R L Tbompeon. and D. J. Wuebblea. J. (jeophyt flee., in 

(28) W. J Mannelll. R H. Neleon. and R S. Johnaton. 1980 American 
Geophysical L'nica Fall Meetinj. San Francieco. CA. Dec 8-12. 1980 

(29) L. F Keyeer. J. Phyt. Chtm.. M. 1659 (19801. 

(30) U. C Sndharon. B Keianonn. and F Kaufmen. J. Chtm. Phyt., 
71. 1288 (19801 

(31) F Zabel. private communication. 1981. 

1 32) R A Cox. R G. Derwent, and p M. Holt. J. Chtm. Soc., Faraday 
Trent. 1, 72. 2031 (1976). 



PNA reaction rate cone tent seems to be unexceptional in 
this group, since no clear trends ere apparent. The A 
factor is appropriate either for abstraction or for dis- 
placement or complex formation. Further work should be 
done on this series of reactions to elucidate the dynamical 
pathways involved. 
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Low pressure photolysis of tert -butyl nitrite at 248 nm with an excimer 
laser includes secondary photolysis of products under conditions where these 
same products do not photolyze. The addition of HI to the photolyzing mix- 
ture suppresses secondary product formation. 

1. Introduction 

An experimental apparatus designed to investigate photolysis pathways 
has been built in our laboratory. A schematic diagram is shown in Fig. 1. The 
very low pressure photolysis (VLP<1>) technique [11 employed allows gas 
molecules to flow through a Knudsen cell at low pressures (10 1 * cm" 2 3 or 
less) where the average residence time is controlled by the rate of escape 
through an aperture [2). The modulated molecular beam mass spectrometer 
analyzes the composition of the effluent gas; the Knudsen cell (or reactor) is 
fitted with high quality Suprasil windows to permit laser irradiation of the 
sample gas. The high sensitivity of the machine (e.g. the detection limit at 
unit mass resolution is about 5 X 10'° cm' 3 ) enables good product detection. 
A Lambda Physics excimer laser was used to study the photolysis of tert- 
butyl nitrite at 248 nm where the photon energy is about 120 mJ pulse' 1 . 

2. Results and discussion 

tert-butyl nitrite (TBN) was chosen as a trial molecule to test the system 
because of itc large cross section [3] at 248 nm and the relative stability of 
its photolysis products: 

CH 3 CH 3 

CH 3 -C-CH 3 NO + CH 3 — C— CH 3 

ONO O 

O 

* CH 3 -C-CH 3 + CH 3 
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Fig. 1. Schematic diagram of VLP4> apparatu*. 


The cell that was used has been characterized and described in the literature 
(1]. Prior to u *», it was thoroughly cleaned and coated with Teflon (DuPont 
?TFE 852-200) according to specifications. The TBN (ICN Chemicals) was 
introduced into the cell at a concentration of about 2 X 10 13 cm” 3 with an 
average residence time of 0.7.5 s. The sample was irradiated at 248 nm (120 
mJ pulse' 1 ; 50 pulses ‘) A mass spectrum was recorded before and during 
photolysis; an example is given in Fig. 2. Care was taken to avoid hitting any 
surfaces other than the windows with the laser beam and to perform a blank 
photolysis run on the empty cell each day. TFr*< did not have a parent peak 
(m/e * 103), so the cracking peak at m/e ■ 88 served to monitor the concen- 
tration of unreacted TBN. Because an integrated mass spectrum was recorded, 
the 50 pulse s' 1 laser repetition frequency was used to maintain a steady 
state concentration of photolysis products in the cell. Our experimentally 
measured absorption cross section for TBN compared favorably with the 
literature value (31 under these conditions. Mass balances were performed, 
and in the absence of any gas other than TBN the acetone yield (m/e - 58) 
was 60% - 80% of the photodissociated TBN. Unfortunately, the mass spec- 
trometer was equally sensitive to m/e * 30 from the NO parent and the TBN 
cracking peak; thus no product information was ascertained at this mass. 

A product peak at m/c * 28 appeared, and if it was interpreted as CO 
formation its yield was 10% - 20% of the dissociated TBN. Some peak inten- 
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F'g. 2. Massspectrum of TBN : , unphotolyzec! TBN; , TBN undergoing photol- 

ysis. 

sity at m/e - 15 was assigned to CH 3 although no CH 4 (m/e = 16) was ob- 
served. At higher concentrations of TBN ( 10 14 cm -3 or more) the percentage 
photodissociated TBN decreased as a result of recombination competing with 
acetone formation from the fert-butoxy radical [4]. Also, ferf-butanol (m/e = 
58) was formed as a product (about 10% yield) at the higher TBN pressures 
which was consistent with secondary reactions becoming important relative 
to acetone formation and/or secondary photolysis. However, the acetone 
yield was not improved over that from the lower pressure runs and no CH 4 
was measured. 

In order to quantify the production of CH 3 , HI was added to the cell 
in varying amounts and the concentration of CH 4 was monitored. Several 
effects occurred during the photolysis of TBN in the presence of HI (2 X 
10 14 - 2 X 10 1S cm -3 ): (1) the acetone yield increased to about 100%, within 
experimental error, although the amount of photodissociated TBN did not 
change; (2) the peak at m/e = 28 (CO) was no longer observed as a product; 
(3) the CH 4 yield varied between 50% and 100%, giving about 100% yield 
only at high concentrations of HI (10 ,s cm -3 or more). Many runs were per- 
formed to verify the reproducibility of these observations. Lower concentra- 
tions of HI that should have been more than adequate to trap CH 3 [5] re- 
sulted in neither CH 4 production nor increased acetone yield. Only at those 
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concentrations where gas-gas collisions became equal to or greater than gas- 
wail collisions were the above effects important. This suggested that perhaps 
the HI was acting primarily as a quencher rather than as a titrant. Under our 
conditions, the ferf-butoxy radical could have an excess energy of as much as 
60 kcal mol -1 which may have been sufficient to access the first excited elec- 
tronic state of the radical. Two different experiments were performed, one 
in which 0 2 was substituted for HI and the other in which CH 4 was used 
instead of HI; in both cases the acetone yield resembled that of the pure 
TBN experiments although the m/e = 28 (CO) production was not significant. 
These data may indicate that an excited electronic state is not involved as 
the addition of relatively efficient electronic quenchers such as 0 2 and CH 4 
should have increased the acetone yield. Perhaps HI is a much more efficient 
quencher. 

To check for complications from absorption of products by the Teflon- 
coated surface (which could have been reduced by monolayers of HI), the 
Teflon was completely removed and the photolysis of TBN was repeated in 
the Pyrex cell. The acetone yield still remained at about 70%. Although this 
test did not absolutely rule out product loss on the walls, it is rather unlikely 
to have been a significant process and remain unchanged by the absence of 
the Teflon surface. 

Reaction of the methyl radical with the excited terf-butoxy radical 
seemed inconsistent with the facts that no CH 4 was observed except with HI 
addition and that several lower TBN pressure runs (about 10 12 cm' 3 ) did not 
increase the acetone yield. Secondary photolysis could certainly have been 
a complicating process in this study, given that the m/e = 28 peak, if inter- 
preted as CO, nearly compensated for the incomplete acetone yield and that 
this peak disappeared when the acetone yield increased to about 100% with 
HI present Several runs were made in which acetone was introduced into the 
Teflon-coated cell (about 10 13 cm" 3 ) and was photolyzed under the same 
conditions as TBN; there were no changes in the mass spectra nor was there 
any apparent dissociation (m/e = 58). The laser repetition frequency was 
reduced to 10 Hz and the TBN spectra were again recorded. Clearly a single- 
shot experiment with real-time product analysis would have been the ideal 
test for secondary photolysis, but the above experiment permitted the same 
steady state product analysis with a reduced photon flux in the same experi- 
mental configuration. The negative results do not rule out the importance of 
secondary photolysis since several laser shots may have been adequate to 
produce our observations. It does seem unusual that HI addition [6] should 
have reduced this process, especially in view of the many wall collisions that 
occurred in the pure TBN experiments. However, the conclusion that sec- 
ondary photolysis took place is consistent with much of the experimental 
data. 


This work was supported in part by the National Aeronautics and Space 
Administration under Contract NAS-7-100/JPL-954815. 
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SUMMARY OF PRELIMINARY DATA FOR THE PHOTOLYSIS OF ClONOj AND N 2 0 5 AT 285 nm 

P. L. Trevor, M. J. Rossi, and D. M. Golden, Department of Chemical Kinetics, 
SRI International, Menlo Park, CA 94025 


A series of preliminary experiments were conducted to determine the 

primary photolysis pathways of N 2 O 5 and CIONO 2 at 285 nm. The study was 

performed under low-pressure conditions, using the SRI Quantei laser system 

aud an apparatus (see Figure 1) specifically built for such studies. 

The low-pressure photolysis (VLP$ ) technique^ - has been described 

previously. Basically, the sample gas flowed through the Knudsen cell 

1 / *1 

(equipped with white cell optics ~ 50 passes) at low pressures (£ 10 i4 cm'*' 5 ) 
where the average residence time was controlled by the rate of escape out of 
an aperture . 2 The composition of the effluent gas was analyzed by the modu- 
lated molecular beam mass spectrometer. The high sensitivity of the machine 

Q _ *1 

(e.g., detection limit at unit mass resolution: 5 x 10* cm for CIONO 2 at 

m/e » 46) enabled good product detection. 

The Knudsen cell, as shown in Figure 2, can be characterized by the 
following: k e min aperture - 8.83 x 10 " 2 (T/M) 1 / 2 , vol - 300 cm 3 , 

length mirror _ Tn £ rror - 29 cm. The stopcock attached directly to the cell 
metered in varying concentrations of 0 ( 3 P) from a steady-state microwave 
source (2450 MHz). The cell was also appropriately equipped for resonance 
fluorescence detection 3 of both 0 ( 3 P) and Cl( 2 P 3 / 2 ) on a real-time, single- 
laser-ahot basis. The CaF 2 windows in front of the resonance lamp and photo- 
multiplier tube (EMR 542G-08-13) conveniently discriminated against Lyman-a 
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radiation. The region between the lamp and the cell was flushed with a mix- 
ture of 1% O 2 in He In order to filter out 135 nm-175 nm light; the area 
between the photomultiplier tube and the cell was flushed with dry nitrogen. 

A CaF 2 lens (fl - 5 cm) waa placed Immediately In front of the resonance lamp 
to focus the light Into the center of the cell. Single photon pulses from the 
photomultiplier went to an amplif ler/dlscrimlnator (PRA Model 1763) whose out- 
put waa stored In a multichannel analyzer (Nlcolet Model 1072). 

The white cell optics consisted of two mirrors (CVI Laser Corporation), 
2-inch diameter, 0.3-m radius of curvature. The front optic was made from an 
AR-coated optosll base which was coated for high laser power at 290 nm * 20 nm 
over the entire surface except for the 0.5-cm diameter laser entrance area at 
the very edge. The back optic was completely coated with the same hlgh-power 
dielectric coating. The mirrors were carefully glued into aluminum mounts 
which were sealed to the cell with vlton o-rlngs. 

The laser light was aligned by eye into the cell and onto the optics in 
such a way as to produce the brightest, most uniform, eliptical pattern of 
approximately 50 multiple passes.^ The laser output from the frequency- 
doubled Quantel dye laser consisted of 10 nsec, 8 mJ pulses, generated at 10 
hertz. The beam was steered Into the cell with four mirrors (CVI Laser Cor- 
poration) coated for high laser power at 290 nm ± 20 nm; the coated mirrors 
also served to filter out the undoubled laser light (570 nm). A power meter 
was placed In the beam immediately after the second steering optic, and the 
total energy did not change by any measurable amount when the undoubled light 
was blocked with a cutoff filter (Corning 7-54). 


The reagent and citrant gases were Introduced from an associated vacuum 
line through the capillary inlets, and the flow rates were measured using a 
calibrated 0.1 psi Validyne transducer. Just prior to data collection, the 
surface of the pyrex cell was silanized. 

Photolysis of C10N0 ? 

The possible product channels involved in CIONO 2 photolysis at ~ 300 nm 
are given In Table 1. Several studies 3- ® at various wavelengths have found 
that path (b) is the dominant channel, while others^ »*° have concluded that 
path (d) is the important one. Thus, a reliable, quantitative measure of any 

O-atom production from the photolysis of CIONO 2 was an essential part of this 

/ 

experiment. 


Table 1 ClONOj PHOTOLYSIS PRODUCT CHANNELS 


Product 

Channel 

C10N0 2 * CIO + NOj 
C10N0 2 * 0( 3 P) + C10N0 
C10N0 2 ♦ 0(*D) + CIONO 
ciono 2 ♦ Cl + N0 3 


Thermodynamic 

Threshold ** 

1108 nm 

509 na 

282 nm 

722 nm 


The resonance fluorescence detection system was first checked by monitoring 
the signal from 0( 3 P) generated on c steady-state basis. A mixture of 10Z N 2 
in He was dissociated in a u-wave discharge (2450 MHz), and a fraction of this 
flow was metered directly into the cell. Excess NO was added to the cell and 
the N-atoms were converted to 0( 3 P) (N + NO ♦ N 2 + 0) vary quickly.* 3 The 
Instantaneous concentration of 0-atom was determined by the Joss of NO obser- 
ved at m/e - 30 for which the mass spectrometer was previously calibrated. 


For a period of several hours, Che plot of Che resonance fluorescence signal 
Intensity remained reproducibly linearly proportional Co Che concentration of 
Che reacted NO; unfortunately , after longer times, the slope of such a plot 
changed. In order to eliminate the poaslblllty that the Instantaneous measure 
(A [NO]) of [0] was not representative of the true steady-state [0] In the call 
because of complicating wall reactions, a mixture of 10Z Oj in He, Instead of 
the N 2 /He mixture, was discharged and metered Into the cell. 

The [0] was measured by titrating with excess NO 2 and monitoring the cal- 
ibrated signal at m/e - 46 at the mass spectrometer. The decrease In the 

1 o 

[NO 2 ] gave the steady-state [0] by NO 2 + 0 ♦ NO + 02* Although there was no 
major difference In fluorescence Intensity, the same problem with long-term 
reproducibility of the resonance fluorescence signal Intensity at a set [0] 
persisted. Not without precedent, the resonance lamp UV-output slowly varied 
during the day, regardless of the precautionary measures taken such as cold 
trapping for Impurities In the makeup gases. 

In order to avoid this problem in addition to calibrating the system 
under conditions as similar as possible to the CIONO 2 photolysis runs, a known 
concentration of NO 2 was photolyzed at 285 nm In the multipass cell 
(NO 2 — * NO + 0), and the time-dependent resonance fluorescence signal was 
averaged over many laser shots ( ~ 10,000), and recorded. An example of such 
a signal is shown in Figure 3. 

At NO 2 concentrations of 3-10 x 10 cm , the decay rate of the 0( P) 
signal was consistent (within experimental error) with the reaction rate for 
NO 2 + 0 ♦ NO + 02»^ The signal was too small to be significant, without more 
time-averaging, when [NOj] <~ 3 x lO 1 ^ cm - ^. The Instantaneous intensity (I Q ) 
of the 0(^P) signal was determined by extrapolation of the exponential decay 
to time ■ 0; It was then noted along with the NO 2 concentration. I Q versus 



[NOj] was noc measured over a wide concentration range due to the time llmlta 
tlon on these preliminary experiments. However, there was no evidence that 
the resonance fluorescence signal was not well behaved. Although both the 
Irradiated volume and the actual laser power In the cell could not be accur- 
ately measured, reasonable estimates were made and an upper limit to the 

Instantaneous 0 ( P) concentration was determined using the reported value for 
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the NO 2 absorption cross section at 285 nrn and a quantum yield of one. No 

more than 1Z of the NO 2 within the Irradiated area was dissociated per pulse, 
thereby yielding 0( 3 P) £ 5 x 10 11 cm -3 at typical NOj concentrations 
(l.a., 5 x 10 13 cm" 3 ). At those 0( 3 P) concentrations, no significant devi- 
ation from proportional linearity was expected"* for I 0 versus [0]. Also, the 
steady— state dissociation of NO 2 (~ 2 Z) as observed at the mass spectrometer 
(m/e ■ 46) provided a good consistency check on the above estimates. 

NO 2 was photolyzed both before and after a C.IONO 2 photolysis run was 
made. Thus the system was phenomenologically calibrated for determining the 
quantum yield of 0( 3 P) product for any other molecule such as ClONOj whose 
absorption cross section 13 at 285 nm is known. Errors resulting from drift in 
the lamp output were minimized with this arrangement. 

The CIONO 2 sample was prepared according to an established method. 1 ^ 
Before purification, It contained substantial levels of CI 2 (m/e ■ 70, 72, 

74), CIO 2 (m/e - 67, 69), and C1 2 0 (m/e - 86 , 88 , 90). After line passivation 
and partial purification, ClONOj (~ 5-8 x 10 13 cm -3 ) was introduced into the 
cell and photolyzed at 285 nm. A CIONO 2 mass spectrum (steady— state) before 
and during photolysis is shown in Figure 4. CIONO 2 had no measurable parent 
peak (m/e • 97,99). Several results are particularly noteworthy, even with 
this contaminated sample: (1) both HNO 3 (m/e ■ 63) and NO 3 (m/e * *>2) were 

produced; ( 2 ) after normalizing relative to NO 2 , the quantum yield for 
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O(^P) production from C10N0 2 appeared to be greater than one; (3) the percen- 
tage of CIONO 2 (®/« “ ^ 6 ) photolyzed was ~ 15Z, lnatead of the anticipated 
2-3Z. Theae reaulta certainly confirmed the ability to detect product NO 3 
(m/e - 62, with time-averaging on the peak), aa Indicated In Figure 3. How- 
ever, the photolysla of the lmpurltlea and the aecondary reactlona of thoae 
producta with C10N0 2 (see Table 2) were obvloualy obecurlng any primary 
product Information. 


Table 2 RATE CONSTANTS FOR SOME C10N0 2 AND I » UR ITT REACTIONS 

Abaorption 

Reaction Croaa Section (T -g, 298 K) Ref . 

C10N0 2 + hv (285 nm) ♦ Cl + NO 3 8.3 x 10 “ 20 cm 2 12 

Cl 2 + hv (285 nm) * 2C1 5 x 10 " 20 cm 2 14 

NOj + hv (285 nm) ♦ NO + 0 6.99 x 10 " 20 cm 2 12 

C10 2 + hv (351 nm) * CIO + 0 ~ 1.2 x 10 " 18 cm 2 15 

C1 2 0 + hv (285 nm) ♦ CIO + Cl 1.1 x 10 " 18 cm 2 16 

Reaction Rate Constant (T^. 298 K) Ref . 


0 + C10N0 2 * Products 2.0 x 10~ 1 ^ molec -1 s“ l 12 

Cl + ClONOj * Products 2.2 x 10 "' 8 cm^ molec -1 s -1 12 


The C10N0 2 sample was further purified by distillation at 175 K to remove 


Che Cl 2 and by pumping at 195 K to remove the N0 2 , C10 2 , and C1 2 0. The pump- 
ing at 195 K was continued until the parent peaks of the impurities were very 






minimal and the m/e ratio of 46/51, 53, no longer changed. The ClONOj 
reservoir was left at 195 K (vapor pressure ~ 1 torr) for the photolysis 
experiments. A much cleaner sample of CIONO 2 (~ 5 x 10^ cm~^) was photol- 
yred, as given in Pigure 6 , which yielded the following results: (1) the 

quantum yield for 0( 3 P) product, normalized to N0 2 , was very 3 mall (~ 0.06); 

(2) Cl 2 (wall recombination of Cl) and a slight amount of HQ were produced; 

(3) NO 3 and HNO 3 were formed. 

13 

The mass spectrometer was locked onto NO 3 (m/e • 62); NO (~ 2 x 10 
cm" 3 ) was added to the cell, and the NO 3 was titrated away via the reaction 
NO + NO 3 ♦ 2N0 2 .^ As indicated in Figure 7, this confirmed the fact that 
m/e ■ 62 was indeed due to NO 3 . Prior to these experiments, the mass spectro- 
meter had been calibrated with HNO 3 , which gave an m/e ratio of 63/62 ^ 100/1, 
where the same unit mass resolution was maintained. HNO 3 was clearly not con- 
tributing to m/e ■ 62, although the HNO 3 formation during photolysis (sue 
Figure 5) may have been the result of NO 3 finding H-atoms on the cell surface, 
analogous to the formation of HC1 when Cl is present. 

Since the CIONO 2 purity cannot b*' determined to better than ~ 90Z by the 
method described above, the 0.06 0( 3 P) quantum yield is a conservative upper 
limit to path (b) and most likely is the result of the photolysis of impur- 
ities (NO 2 and/or C 102 )* 

In order to eliminate the possibility that the NO 3 was being formed by 
the reaction Cl + CIONO 2 ♦ Cl 2 + NO 3 , ethane was added to the cell ([C 2 Pg] 

- 7 x 10« cm -3 ) to preferentially react away Cl according to the reaction 
Cl + C 2 H 6 ♦ HC1 + C 2 H 5 »^ Although the ethane reaction was not a "clean" 
titration reaction in that the C 2 H 5 could further react in the cell, even with 



(NO pump away) 
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FIGURE 7 EXPANDED VIEW OF N0 3 <m/e - 62) FORMATION AND CHEMICAL 
TITRATION TESTS 



ClONOj (k e ■ 0.155 s" 1 ), Figure 7 shows that the [NO3] remained the same, 
which strongly suggests that the NO3 came directly from the photolysis of 

ciono 2 . 

In summary, this preliminary study of the photolysis of CIONO2 at 285 nm 
indicates that Cl + NO3 Is the major, if nor. the only, product channel. 

Photolysis of 

N2O5 was photolyzed at 285 nm in exactly the same manner as C10N0 2 . The 
0(^P) detection sensitivity was calibrated with NO2, as described above. The 
two most likely product channels for N2O3 photolysis at 285 nm and their 
thermodynamic thresholds are given in Table 3. 


Table 3 


Prouuct Channel 
N 2 0 5 + hv ♦ N0 2 + NO3 
N2O5 + hv ^ 2NO2 + 0 


Thermodynamic 
Threshold^ * 
1288 nm 
393 nm 


One study^ has concluded that pathway (f) is the product channel at 254 nm. 
Therefore, serious efforts were again made to monitor any 0(^P) product formed 
during photolysis. 

IQ 

NjOg was prepared by reacting N0 2 with excess O3 and trapping the 9olid 
N2O5 into a pyrex vessel at 195 R. Before use, the was trap-to-trap dis- 

tilled and the NO2 was pumped off at 250 R. The ^2^5 reservoir was maintain- 
ed at 250 R (vapor pressure ■ 5 torr) during the experiments. appeared 

to be very reactive on the pyrex surfaces of the vacuum line and inlets, even 
after hours of passivation, as previously observed by others. To hasten tht 


passivation process and to prevent N 2 O 5 from substantially converting to HNO 3 , 
the reservoir and needle valve were moved as close as possible to the cell, 
with a bypass to the vacuum line only for meassurlng pressure and flow rates. 

This procedure reduced the HNO 3 Impurity level to - 20Z; lower levels were not 

1 2 

achieved. However, due to the low absorption cross section at 285 nm for 
HNO 3 ( 8.8 x 10 ~ 21 cm 2 ), no complications from its photolysis products were 
anticipated. 

The NOj impurity proved to be a more serious problem, since it was 

18 

continually being formed from N 2 O 5 decomposition. After each run, the gas 
above the reservoir before the ueedle valve, and the NO 2 , were pumped away at 
250 K. 

A mass spectrum of (~ 5 x lO 1 ^ cm~^) is given in Figure 8 . N 2 O 5 

exhibited no parent peak (m/e ■ 108). A very small peak was meausred at m/e • 

20 21 

62 (m/e ratio 46/62 ~ 900/1), which had not been previously observed. * 
jj 0 (~ 3 x cm “3) was added to the cell to confirm that it was a fragmenta- 
tion peak from ^Oj and not the parent peak of the NO 3 decomposition pro- 
duct. At no time was there any measurable change in the intensity at m/e ■ 62 

in the presence of NO. 

22 —20 2 
The absorption cross section for N 2 O 5 at 285 nm (8.4 x 10 cm ) is 

very similar to that for C 10 N 02 * Normalized to NO 2 , again only 2-3Z photol- 
ysis was expected. Given both the weak signal intensity at m/e ■ 62 from N 2 O 5 
and the high probability f or similar mass spectral sensitivities for m/e ■ 62 
from N 0 O 5 and NO 3 , any NO 3 photolysis product would probably have not been 
identifiably measurable mass spectrometrically , even with long time averag- 
ing. Nevertheless, the peak at m/e » 62 was carefully monitored during 
several photolysis runs; no significant change in the signal Intensity was 



noted. A 2-3Z change In the peak Intensity would not have been measured with 
any certainty, given the low slgnal/nolse. Thus, this observation was consis- 
tent with either product channel. 

There was very little change In any part of the mass spectrum during 

photolysis, which was expected because of product contributions at the major 
peaks. The peak at m/e ■ 46 was recorded (with long time-averaging), which 
Indicated a very slight decrease (< 1Z) during photolysis. From prior cal- 
ibrations, the mass spectrometer was found to be ~ 20Z more sensitive to 
than to NO 2 at m/e ■ 46. Sensitivities relative to NO 3 at m/e ■ 46 could not 
be determined or well estimated. Unfortunately, these data (or lack thereof) 
at m/e • 46 did not provide any confirmatory evidence as to the dominant 
product channel and, again, were consistent with either pathway. 

Our final experimental efforts were directed at determining the quantum 
yield for 0("P) production SO 2 resonance fluorescence calibration runs were 
conducted both before and after each N 2 O 5 experiment, just as In the CIONO 2 
study. In addition, after several of the NjOj ruas » steady-state [0] 

(~ 1 x 10 ^ cm*^) was added to the cell from the u-wave Inlet In order to 
titrate any NO 2 Impurity without reacting away any of the N 2 O 5 (k g ■ 0.147 
g -l).20 j n a j_j_ ^ Jut Qne titration, ~ 10Z of the Intensity at m/e 46 reacted 
away, indicating a significant NO 2 impurity level. In one of the titrations, 
as much as 20Z of m/e ■ 46 disappeared. (NO 2 was added afterwards to verify 
that excess [ 0 ] was still present, even in the case with high NO 2 impurity). 
Under these experimental conditions, an 0( P) quantun yield £ 0.10 would most 
likely be due to NO 2 photolysis rather than N 2 O 5 photolysis. 


After normalization, the resonance fluorescence decays, an example of 
which Is given In Figure 9, yielded an 0( 3 P) quantum yield ~ 0.10 (within 
experimental error). No single N2O5 run (for which the data were collected) 
Indicated a higher quantum yield than the NO2 Impurity level; however, the 
purity levels were not checked on all the resonance fluorescence experiments. 

The NjOj concentration was not systematically varied In our data. None- 
theless, at [N2O5] ; 5 x 10 13 cm” 3 only minimal, If any, quenching of the pho- 
tolysis an proposed by Connell^ would have been taking place. Clearly, this 
is an area that requires further Investigation. 

In summary, since no significant 0( 3 P) product can confidently be 
assigned to N2O5 photolysis, and, because there la no evidence to the 
contrary, the product channel giving NO2 + NO3 Is Indicated as the primary 
pathway for the photolysis of N2O5 at 285 nm. 
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FIGURE 9 0( 3 P- RESONANCE FLUORESCENCE DECAY 
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INTRODUCTION 

Nitrogen oxides of natural and of anthropogenic origins can interact 
In the stratosphere to form dinitrogen pentoxlde. Once formed, NjOj acts 
as a short-lived reservoir species, Its major fate being photolysis to 

regenerate odd oxygen species.* It Is known from laboratory studies that 

2 

N 2 Oj hydrolyzes rapidly on surfaces and very slowly in the gas phase, and 
it has been recognized that hydrolysis of NjOg on stratospheric aerosol 
particles may produce nitric acid, which can rain out of the atmosphere. 
Because so much of the atmospheric NO x spends short periods of time tied up 
in N^j, a realtively minor reaction pathway to produce HNO 3 can have an 
important quantitative effect on stratospheric composition and the ozone 
distribution. 

The stratospheric aerosols have been determined to consist primarily 

of sulfuric acid solutions in water (typically, 70-90 weight percent 

^SO^).^ The 'Mxlinum aerosol number density occurs at about 20 ka with 

4 

aerosol diameters ranging from 0.01 to 1.0 tim. The structure of the 
aerosol particles is not well characterized, but they oav consist of frozen 
super-cooled liquid cores surrounded by a thin film of liquid.^ 


1 


It Is desirable co determine Che efficiency with which N2O5 reacts 
with f^SO^/l^O mixutres in order to assess whether such reactions can play 
significant roles in stratospheric chemistry. In a previous survey^ by 
this laboratory of heterogeneous reaction involving 1002 H2SO4 surfaces, 
the reaction probability (y) of ^2^5 was measured to be 3.3 x 10“^, and 
there were indications that a small quantity of NO2 was released as a gas- 
phase product. In the same study, the expected hydrolysis product, HNO3, 
was found to react with y • 2.4 x 10~ 4 . These reaction efficiencies are 
too low to be important in stratospheric chemistry, but the question 
remained of the effect of adding water vapor to the system. This question 
is addressed in the present study. 


EXPERIMENTAL 


The experimental apparatus and technique have been described in detail 
elsewhere.^* ^ Briefly, the apparatus consists of a two-chamber Knudsen 
cell reactor (Figure 1) coupled with a modulated molecular beam mass spec- 
trometer (MMBMS), which has been described in the literature. A regulated 
flow of reactant gas enters the reactor through one of the two capillary 
inlets and a flow of water vapor can be independently introduced through 
the other. Absolute flow rates are measured by timing pressure Increases 
in calibrated volumes. The total pressure in the reactor is maintained 
below about 10 mtorr so that the mean-free path of the gas molecules is 
greater than the diameter of the escape aperture, ensuring that Knudsen 


flow through the aperture controls the residence time (x) in the reactor. 
The volume of the reactor is 427 cnr and the escape rate constant (1/x) is 
k e ■ 1.21/±M^2at room temperature for species of molecular weight M. The 
composition of the effluent gas is monitored with the MMBMS , and the 
results are recorded on a strip-chart recorder or are stored in the memory 
of a signal averager. 

The main chamber of the two-chamber reactor is isolated from the sul- 
furic acid surface by an externally operated valve. When the valve is 
opened, the H 2 SG 4 surface is exposed, and the reactant gas may undergo 
reaction, altering the composition of the effluent gas monitored by the 
mass spectrometer. The mass spectrometer used in the series of experiments 
reported here was a Finnigan model 400 quadrupole. 

Dinitrogen pentoxide was synthesized by mixing ozone generated in a 
stream of pure oxygen with a small amount of NO 2 and condensing the N 2 O 5 at 
195 K. Ozone was always present in excess and the resulting ^Og was ini- 
tially free of NO 2 , although small amounts of HNO 3 were present, presumably 
due to heterogeneous reaction of with traces of H 2 O. The pure anhyd- 

rous nitric acid used for mass spectrometer calibration was synthesized by 
mixing NaNOj with concentrated H 2 SO 4 and condensing the product gas at 195 
K. Nitrogen dioxide was mix ' ernight with a great excess of pure oxygen 
to convert any nitric oxide to and was then thoroughly degassed by 
bulb-to-bulb distillation. 

The sulfuric acid used as the heterogeneous surface was placed in a 
pyrex watch-glass to form a pool of circular diameter 4.8 cm. The 
concentrated acid was pumped at < 10 ”^ torr over-night prior to use to 
insure that it was thoroughly degassed and as dehydrated as fecsible. 


RESULTS AND DISCUSSION 


Typically, the freshly synthesized N 2 O 5 showed little evidence of NO 2 
contamination, but after 30 minutes, the mass spectrum had changed signifi- 
cantly. Three mass; peaks were routinely monitored: m/e ■ 30, 46, 62/63. 

The sensitivity of the mass spectrometer was too low to permit unit mass 
resolution at m/e » 62/63 and total intensity of both peaks was measured. 
Since no mass 62 peak was observed fci pure N 2 O 5 , but there was a measur- 
able contribution from pure HNO 3 , signal at this mass was indicative of the 
presence of HNO 2 . However, the mass spectrometer sensitivity was so poor 
and the signal was so noisy that quantitative measures of HNO 3 could not be 
obtained. 

In the absence of added water vapor, NjOj reacted with the sulfuric 
acid surface in a consistent manner, such that 40 ± 10% of the N 2 O 5 was 
lost to the surface. There was no evidence that gas was released as a 
product, contrary to the survey results retorted earlier. 

For this experiment, 



I 



where I and I Q are the triss spectrometer signal intensities without and 
with -he sulfuric acid exposed, respectively. i^ the < fective area of 


For 


the escape aperture and A g is the area of the sulfuric acid pool.^'^ 
the present series of experiments, A^/A g ■ 4.5 x 10~ 4 . Thus, in the 
absence of 1^0 vapor, y - 3 x lO - ^, about an order of magnitude larger than 
reported earlier and almost identical to the result obtained earlier for 
HNO-j. 6«7 This result was found to be independent of N 2 O 5 concentration, 
and it was reproduced using several different batches of The 

uncertainty on y is no larger than ± 50%. 

When water vapor was added to the system, a contribution was observed 
at m/e ■ 62/63, Indicative of the formation of HNO^. Unfortunately, the 
signal was too noisy to be of quantitative usefulness. Analysis of the 
results indicated that as much as 50% of the NjOj was converted to HNO 3 
when large amounts of H 2 O vapor were introduced. Despite this complica- 
tion, the loss of N 2 O 5 on the H 2 SO 4 surface seemed to be substantially 
unaffected by the presence or absence of water vapor. It was also conclud- 
ed that HNO 2 is lost on the H 2 SO 4 surface w about the same y as measured 
for N 0 O 5 (consistent with the earlier results). 

The reason for the discrepancy between the present results for ^05 
and those reported earlier is not known, but it may be related to the way 
concentrations were monitored in the earlier work.^i^ Pre’ iously, only m/e 
* 46 was monitored even though the system can contain N 2 O 5 , hC> 2 , and HNO 3 , 
all of which give contributions at that mass. In the present experiments, 
the situation is a little better, but not completely satisfactory due to 
the poor sensitivity for m/e = 62/63. In future experiments on this reac- 
tion, a tetter mass spectrometer will be employed. 

If we accept that for ^05 y ■ (3 ± 1.5) x 10“^, we can estimate the 
effect of heterogeneous reactions on loss of N 2 O 5 * n the stratosphere. For 


typical concentrations and sizes of stratospheric aerosols, the collision 

frequency of an N 2 O 5 molecule with the aerosol particles has been estimated 

to be ~ IO - 5 s - 1 . 6 * 7 Thus, the total first-order rate of loss of N 2 O 5 due 

—9 —1 

to heterogeneous reaction is"' 3 x 10 s . This value can be compared 
with the photolysis loss rate for NjOg at 20 km, ~ 5.3 x lO - -’ s - ^,® whi.h 
is about four order s-of-magnitude faster. Unimolecular decomposition of 
N 2 O 5 tc form NO 2 and NO 3 has a rate constant of ~ 2 x 10 - ® s - ^ at 20 km,® 
an order of magnitude faster than the estimated heterogeneous loss rate. 
Thus, we conclude that heterogeneous reaction of N 2 O 5 plays only a small 
role at night and is insignificant when photolysis is possible. 
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